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Notations 



This is a table of symbols and acronyms used in the thesis. Sign conventions and notations 
are as in the classic book Gravitation by C.W. Misner, K.S. Thorne and J. A. Wheeler if 
otherwise not stated. In particular Greek indices may take values — 3 and Latin indices 
may take values 1 — 3. 



SYMBOL NAME 

KE/NKE kinematical/non-kinematical evolution (of the FHSs) 

FOs/FHSs fundamental observers/fundamental hypersurfaces 

NKR/NKCE non-kinematical redshift/non-kinematical coupling equation 

GTCS global time coordinate system 

HOGS hypersurface-orthogonal coordinate system 

WEP/EEP weak/Einstein equivalence principle 

GWEP/SEP gravitational weak/strong equivalence principle 

LLI/LPI local Lorentz/positional invariance 

SR/GR special/general relativity 

SC Schiff's conjecture 

t — t{x^) global time function, conventionally we define t = x°/c when 

using a GTGS (where is a global time coordinate) 

Ai/Ai metric space-time manifolds obtained by holding t constant 

J\f/J\f quasi-metric space-time manifolds, foliated into the set of 
FHSs 

or g(^t)fiu one-parameter family of space-time metric tensors on Af 

gt or g(t)ij.u one-parameter family of space-time metric tensors on J\f, so- 
lution of field equations 

gt or g(t)nu one-parameter family of auxiliary space-time metric tensors 
on AT 

lit or n^*^-) unit vector fields normal to the FHSs in the metrics g^ 

ht or n^^^j unit vector fields normal to the FHSs in the metrics gj 

h( or h(^t)ij spatial metric tensors of the FHSs obtained from gt 

ht or h(^t)ij spatial metric tensors of the FHSs obtained from gt 

N{x^) lapse function field in the metrics gj 

Nt(x^, t) lapse function field family in the metrics gt 
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N^{x^) components of the shift vector field in the metrics gtj, 

N^{x^) components of the shift vector field in the metrics g^g 

r"j)^^ components of the metric connections Vt compatible with gt 

r"j)^^ components of the metric connections Vt compatible with gj 

r"^^^^ components of the metric connections Vt compatible with gj 

r^, components of the degenerate 5-dimensional connection V 

associated with the family gj 

* * ± 
r^, components of the degenerate 5-dimensional connection V 

associated with the family gj 

or a^j-j family of metric 4- accelerations (any observer) 

coordinate expression for degenerate covariant derivation 

compatible with gt 

()jq, coordinate expression for degenerate covariant derivation 

compatible with gt 

coordinate expression for spatial covariant derivation com- 
patible with ht or as appropriate 

coordinate expression for metric covariant derivation com- 
patible with a member of gt or gt as appropriate 
±/l projection symbols (projection on the normal direction to 

the FHSs) 

£x Lie derivative with respect to X in A4/ 

* 

i?x Lie derivative with respect to X in H /J\ 

£x projected Lie derivative with respect to X in /TW 

vCx projected Lie derivative with respect to X in M /'JJ 

a^/ajF or a^j:/a}j: four-acceleration of the FOs ui M./ M. 

hj: or 6^ "local distance vector from the center of gravity" generalized 

from the spherically symmetric case 
e^/e^ or e\/e\ unit 3-vector/covector field in the bjr-dircction 
Wt or w^j^ coordinate 3- velocity of test particle with respect to the FOs 

Ut or u^^-^ 4-velocity of test particle in M 

Uf or ■u|'^j 4-velocity of test particle in AT 

Vt or 3-vector fields giving the motion of FOs in JJ compared to 

FOs in N 
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a or a/* degenerate 4-acceleration in J\f (any observer) 

G{t)iJ.v/G{t)iJ.v Einstein tensors m. M./M. 

R{t)iiv / R{t)iJ.v Ricci tensors in Al/ 

H(^t)ij/H(^t)ij spatial Einstein tensors of the FHSs in A^/Al 

P(t)ij / P{t)ij spatial Ricci tensors of the FHSs m. M/M. 

Pt/Pt spatial Ricci scalars of the FHSs m. M./M. 

or K{t)ij family of extrinsic curvature tensors of the FHSs in M. 

or K(t)i3 family of extrinsic curvature tensors of the FHSs in Jv[ 

Tit the space-time tensor family projected into the FHSs 

Ht/yt local measures of the NKE 

Xt local measure of the KE 

drt/drt proper time interval (any observer) in M 

drjrjdTjr proper time interval for a FO in A/'/AT 

Tt or T^^ stress-energy tensor as an active source of gravitation 

(in AT) 

Tt or T^^ passive stress-energy tensor (in M) 

Gjn Newton/Einstein gravitational constant 

mt/Mt active masses (scalar fields measured dynamically) 

Ft=Ntct local scale factors of the FHSs in AT 

^j=Fj~^ a scalar field describing how atomic time units vary in 

space-time 

—U Newtonian potential 

Pm passive mass density in local rest frame of the source 

p passive pressure 

Pm active mass density in local rest frame of the source 

p active pressure 

Pm the amount of active mass per coordinate volume 

p pressure corresponding to pm 



3 



4 



Prologue 



Observationally the Hubble law in its most famous form may be found by measuring 
spectral shifts of "nearby" galaxies and thereby inferring their motion; it relates the 
average recessional velocity v<^c of such galaxies to their distance d via the equation 

V = Hd, 

where H is the Hubble parameter. We see that the Hubble law does not depend on 
direction. This is merely a consequence of the fact that it is an empirical law; had the 
observations suggested the existence of anisotropic recessional velocities the Hubble law 
could still be formulated but in an anisotropic version. 

The question now is if this apparent lack of anisotropy may follow from some hitherto 
undiscovered physical principle rather than being a consequence of some rather special 
kinematical initial conditions. Once one suspects this it is natural to assume that the 
potential new fundamental law is local in nature. That is, by inspecting the behaviour 
of the Hubble law for small distances d one may get a hint of how it may be understood 
as a special case of a local law. 

We see that the particular form of the Hubble law implies that recessional velocities 
may be made to vanish to arbitrary accuracy by considering sufficiently small distances. 
This observation may be significant since there in fact exists a local geometrical quantity 
the effects of which have just the property that they may be made to vanish to arbitrary 
accuracy in a sufficiently small region. This geometrical quantity is curvature. Hence one 
may get the idea that a hypothetical local version of the Hubble law has something to do 
with curvature. Furthermore the natural way of relating the Hubble law to curvature is 
via the Hubble parameter. Since this is a scalar and not a tensor (as it would have to be 
in an anisotropic version of the Hubble law) one may suspect that the Hubble parameter 
in fact may be expressed as some function of curvature scalars. 

In standard cosmology the Hubble law applies only to a particular set of observers 
associated with the smeared-out motion of the galaxies. This suggests that a hypothet- 
ical local version of the Hubble law should also be associated with a privileged class of 
observers. At least parts of the interrelationship between nearby privileged observers 
should consist of an "expansion" . For reasons explained in the thesis this kind of expan- 
sion is called "non- kinematical" . Furthermore one may guess that the Hubble parameter 
should depend on curvature scalars calculated from the 3-dimensional space defined by 
the privileged observers taken at constant values of some privileged time coordinate. This 
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privileged time coordinate then acts as a global notion of "simultaneity" and should be 
a basic element of the space-time geometry. 

It should be clear from the above that any attempt to construct a local version of the 
Hubble law and implementing this into a general geometrical structure in fact necessitates 
the construction of a new framework of space and time. We show in this thesis that 
this new framework disposes of the space-time metric as a global field. This means 
that the mathematical structure of the new framework is not based on semi-Riemannian 
geometry. Thus not all aspects of metric theory will hold in the new framework. On 
the other hand, important physical principles, previously thought to hold only in metric 
theory, also hold in the new framework. For this reason we call the new framework 
"quasi-metric". Since the quasi-metric framework must be relativistic one may always 
construct a local space-time metric on the tangent space of each event by identifying the 
local inertial frames with local Lorcntz frames. But rather than demanding that the set 
of local metrics constitutes a global space-time metric field we require the less stringent 
condition that the set of local metrics constitutes a semi-global metric field. That is, the 
domain of validity for the semi-global metric field is not the entire space-time manifold 
but rather a 3-dimensional submanifold defined by some constant value of the privileged 
time coordinate. The quasi-metric space-time manifold may be thought of as consisting 
of a family of such submanifolds each of being equipped with a semi-global space-time 
metric field. The crucial fact is that the set of such semi-global space-time metric fields 
does not necessarily constitute a single global space-time (non-degenerate) metric field. 
The reason for this is that the affine connection compatible with the set of semi-global 
metric fields depends directly on the existence of a privileged time coordinate. That is, 
unlike the Levi-Civita connection the affine connection compatible with the set of semi- 
global metrics cannot in general be derived from any single space-time (non-degenerate) 
metric field. 

In this thesis it is shown that the absence of a global Lorentzian metrical field and the 
existence of a privileged time coordinate does not make it necessary to give up important 
physical principles such as e.g. the various versions of the principle of equivalence (but 
the strong principle of equivalence does not hold in its most stringent interpretation). 
On the other hand the generalization of non-gravitational physical laws to curved space- 
time is more complicated in the quasi-metric framework than in the metric framework. 
This should not be surprising since non-gravitational fields may possibly couple to fields 
characterizing the quasi-metric geometrical structure in a way not possible in metric 
geometry, but such that terms representing such couplings vanish in the local inertial 
frames. (That is, in the local inertial frames the non-gravitational physics take their 
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standard special relativistic form just as in metric theory.) Thus, even if there is a 
preferred notion of space and time, this does not mean that there must exist a preferred 
coordinate frame in the sense that the outcomes of local experiments (gravitational and 
non-gravitational) depend on the velocity with respect to it. And this is easily justified 
since the geometric structure of quasi-metric theory does not depend on the existence 
of any particular coordinate frame. On the other hand there exists a set of "preferred" 
coordinate systems especially well adapted to the geometric structure of quasi-metric 
space-time. This is a consequence of the fact that the preferred time coordinate represents 
an "absolute" geometrical element. However, the existence of non-dynamical fields is 
possible in metric theory also. Thus there should not be any a priori objections to 
constructing theories of gravity based on the quasi-metric framework. Consequently, in 
this thesis it is shown how to construct such a theory. This theory corresponds with 
General Relativity in particular situations and may possibly be viable. 

When it comes to comparison of testable models based on the quasi-metric theory to 
experiment, it still remains to develop a suitable weak field expansion analogous to the 
parameterized post-Newtonian expansion valid for the metric framework. Obviously this 
is a subject for further work. However, even if a suitable weak field expansion is missing 
one may still try to construct specific models for idealized situations. In the thesis 
this is done for certain spherically symmetric systems and it is shown that the classical 
solar system tests come in just as for General Relativity. But one important difference 
is that the quasi-metric theory predicts that the gravitational field of the solar system 
expands with the Universe whereas General Relativity predicts no such thing. Moreover, 
gravitationally bound bodies made of ideal gas are predicted to expand in this manner 
also. As discussed in chapter 3, chapter 4 and chapter 5 of the thesis, observational 
evidence for expanding gravitational fields so far seems to favour the quasi-metric theory 
rather than General Relativity. Thus the observational evidence seems to indicate that 
the recessional velocities of galaxies are non-kinematical in nature and moreover that 
space is not flat. If the observational evidence holds up one must have in mind that the 
quasi-metric framework in general and the predictions of expanding gravitational flelds 
and certain gravitationally bound bodies in particular, are based on a reinterpretation of 
the Hubble law. 

It may be surprising that a reinterpretation of the simple empirical Hubble law leads 
to nothing less than the construction of a new geometrical framework as the basis for 
relativistic physics. However, from a philosophical point of view one should prefer a 
theoretical framework having the property that general empirical laws follow from first 
principles. Thus seen the existence of general empirical laws which do not follow from first 
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principles may be interpreted as a sign of incompleteness for any theoretical framework 
and may potentially lead to the demise of that framework. 

Now that the reader has had a foretaste of what this thesis is all about, it is useful 
to get a overview of the mathematical concepts used in the thesis before starting on its 
main parts. 
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Some Basic Mathematical Concepts 

by 

Dag 0stvang 

Institutt for Fysikk, Norges teknisk-naturvitenskapelige universitet, NTNU 

N-7491 Trondheim, Norway 

Abstract 

The canonical viewpoint of space-time is taken as fundamental by requiring the 
existence of a global time function t corresponding to a global foliation of space- 
time into a set of spatial hypersurfaces. Moreover it is required that space-time can 
be foliated into a set of timelike curves corresponding to a family of fundamental 
observers and that these two foliations are orthogonal to each other. It is shown 
that this leads to a "quasi-metric" space-time geometry if t is taken to represent 
one extra degenerate time dimension. The quasi-metric framework then consists 
of a differcntiable manifold equipped with a one-parameter family of Lorentzian 4- 
metrics parameterized by the global time function t, in addition to a non-metric 
connection compatible with the family g^. 

1 Introduction 

In this thesis we construct a new relativistic theory of gravity. This theory is compatible 
with a geometric framework that differs from the usual metric framework. While there 
is nothing remarkable about the mathematical description of the new "quasi-metric" 
framework (it is just ordinary differential geometry), the mathematics and its applications 
may seem unfamiliar to the traditional relativist who is used to working only with metric 
geometry. For this reason we give an intuitive introduction to the basic mathematical 
concepts of the quasi-metric geometry; the physical motivation for introducing it will 
have to wait until the main part of the thesis. 

2 From 3-geometry to 4-geometry 

The basic premise defining metric space-time geometry is that space-time can be de- 
scribed as a 4-dimensional semi-Riemannian manifold. In particular, in metric geometry 
the affine geometry follows uniquely from the space-time metric (barring torsion). Thus, 
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since the affine geometry does not depend on any preferred time coordinate, the view- 
point is taken as fundamental, that there is no need to specify any particular splitting of 
space-time into space and time except for reasons of convenience. And certainly no such 
splitting can have anything to do with fundamental physics. 

Yet one may wish to describe the geometry of space-time via the evolution with time 
of fields defined on a spatial hypersurface as in the ADM formahsm [1] (this is called a 
canonical description of the space-time geometry). However, in metric theory the choice 
of said spatial hypersurface is arbitrary, and its deformation along some arbitrary timelike 
vector field is totally dependent on the exact form of the space-time metric tensor field g. 
In particular this implies that in metric theory, one cannot impose general requirements 
that certain canonical structures should exist without simultaneously imposing restric- 
tions on the form of g. And from a physical point of view this is unacceptable since 
in metric theory, the afiinc structure follows directly from the space-time metric, so any 
restrictions on that will interfere with the dynamics of any potential theory of gravity. 

Suppose then, that for some reason one wishes to take the canonical viewpoint of 
space-time as fundamental but such that space-time needs not be metric. Then there is 
no reason to believe that the requirement that certain canonical structures exist should 
restrict the afiine geometry unduly. In particular, one may require that a preferred time 
coordinate should exist. In this case the affinc geometry docs not necessarily depend 
critically on a space-time metric. We illustrate this in the next section. Meanwhile 
we describe the canonical structures that are required to exist within our quasi-metric 
framework. 

Mathematically, to justify the view that a canonical description of space-time is fun- 
damental, such a description must involve canonical structures that are taken as basic. 
Thus such structures must always exist, and they should in some way be "simpler" than 
canonical structures in metric theory. One should regard these criteria as fulfilled if one 
requires that the space-time manifold H always can be foliated into a particular set of 
spatial fundamental hypersurfaces (FHSs) St parameterized by the global time function 
t. To ensure that the parameterization in terms of t exists independent of any space-time 
metric we let t represent one extra degenerate time dimension. Besides we require that 
H can always be foliated into a family of timelike curves corresponding to a set of funda- 
mental observers (FOs), and such that these two fohations are orthogonal to each other, 
see figure 1. 

If this is to make sense, there must be a space-time metric tensor available at each 
tangent space such that scalar products can be taken. To be sure that such a space-time 
tensor does exist we require that on St, there must exist a family of scalar fields Nt and 
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t=constant 



Figure 1: Basic canonical structure: space-time accepting two mutually orthogonal foliations. 

a family of intrinsic spatial metrics h^. Moreover there must exist a set of preferred 
coordinate systems where the time coordinate x° can be identified with the global time 
function t. We call such a coordinate system a global time coordinate system (GTCS). 
The coordinate motion of the FOs with respect to a GTCS defines a family of spatial 
covector fields tangential to the FHSs. Now A^^^, and uniquely determine a family of 
Lorentzian space-time metrics gj. This metric family may be thought of as constructed 
from measurements available to the FOs if the functions A^^^ are identified with the lapse 
functions describing the proper time elapsed as measured by the FOs, and if the covectors 
are identified with the shift covector fields of the FOs in a GTCS. With the help of A''^ 
and Nt we may then construct the family of normal vector fields rij being tangent vector 
fields to the world lines of the FOs and such that ht(nt, ■) = 0. Thus by construction the 
foliation of M into the set of world lines of the FOs is orthogonal to the foliation of M 
into the set of FHSs, as asserted. To illustrate this we may write 

gt = -gt(nt,-)(g)gt(nt,-) +hj. (1) 

Note that there is no reason to believe that any affine connection compatible with the 
above construction must be metric. Thus the point is to illustrate that it may be possible 
to create a geometrical structure where a space-time metric tensor is available at every 
event but where the affine geometry does not necessarily follow directly from a space- 
time metric. The reason for this is that the existence of the global time function t may 
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Figure 2: The geometrical structure of quasi- metric space-time, visualized. 

be crucial when determining the affine connection; this possibility would not exist if the 
space-time geometry were required to be metric. 

3 Connection and curvature 

The existence of a Lorentzian space-time metric gGg* for each FHS S&St implies that 
each FHS can be viewed as a spatial submanifold of a metric space-time manifold Ai. 
We now perform a mathematical trick inasmuch as whenever we take this viewpoint 
we may define a global time coordinate x° on Ai such that = ct on each FHS; any 
coordinate system where this relationship holds is by definition a GTCS. We may then 
regard x° and t as representing separate time dimensions such that t =constant on Ai. 
It is thus possible to view A/" as a 4-dimensional submanifold of x R. The point with 
this is that on the FHSs we may take scalar products of fields defined on space-time in 
A4 whereas time derivatives of such fields depend on t as well as on in A/". Thus the 
existence of extra time derivatives with respect to t of gt means that we can find an affine 
connection that is not fully determined by any single gGgf. See figure 2 for visualization 
of quasi-metric space-time geometry. 

Since each gGgt is defined on a spatial submanifold only, there is in general no way 
to extend it to all of A/". However, the Levi-Civita connection corresponding to each 
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gegt can be calculated on each FHS by holding t constant. But this family of metric 
connections does not fully determine the connection compatible with the family g^. To 
find the full connection compatible with the family gt it is convenient to view g* as a single 
5-dimensional degenerate metric on AlxR. One may then construct a linear, symmetric 
and torsion- free connection V on AlxR by requiring that this connection is compatible 
with the non-degenerate part of the metric family g^ in addition to the condition that 
Va.nt vanishes. For obvious reasons this connection is called "degenerate", and it should 

dt 

come as no surprise that it is non-metric. One may then restrict the degenerate connection 
to jV by just considering the submanifold x'^ — ct in a GTCS. For the exact form of the 
connection coefficients determining V, see e.g. [2]. 

At this point a natural question would be if there are physical phenomena which are 
naturally modelled by the non-metric part of quasi-metric geometry. And as we shall see 
in the main part of this thesis this is indeed the case. 



4 Two simple examples 

We finish this mathematical introduction by giving two simple examples. The first ex- 
ample illustrates how a family of space-time metrics can be represented by a family of 
line elements, whereas the second example illustrates how affinely parameterized curves 
are represented. 

Example 1 

A family of Minkowski metrics can be represented by the family of line elements 

dsl = -{dx^'f + 52(i) (dx" + dy"" + dz^) , (2) 

where S(i) is a scale factor, x^ — ct and x, y, z are ordinary Cartesian coordinates. Note 
that the t-direction is treated as degenerate. It is easy to show that the free fall curves 
obtained from the equations of motion (derived in chapter 1, section 5 of this thesis) and 
applied to (2) are different than their counterparts obtained by setting dx^ — cdt in (2) 
and using the geodesic equation compatible with the resulting single metric. 

Example 2 

Let {t{X),x"{X)} represent an affinely parameterized curve, where x°' are space-time 
coordinates in a GTCS. Then the tangent vector field along the curve may be represented 
as 

d dt d dx^ d . . 

dX ^ dX'di^ ~dX^' ^ ' 
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Now the length of ^ may be calculated along the curve, e.g. for a timehke curve 
gt(^, ^) < along the curve. On the other hand the degenerate connection taken 
along ^ is 



* dt* dx'^* 

9^ dX St dX Bxf^ 

Thus we have illustrated that the degenerate part of ^ does not contribute to its length 
(since gt(^, ■)=0)- However, as can be seen from equation (4), the degenerate part of ^ 
does influence parallel-transport of space-time objects. 

Hopefully, anyone who has understood the basic mathematical concepts underlying 
the quasi-metric framework will not be so easily distracted by mathematical hurdles and 
will be in a better position to follow the derivation of the main results presented in this 
thesis. Accordingly, having read this mathematical introduction, the reader should now 
be ready for the main parts of the thesis. 

References 
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A Non-Metric Approach to Space, Time and Gravitation 



by 

Dag 0stvang 

Institutt for fysikk, Norges teknisk-naturvitenskapelige universitet, NTNU 

N-7491 Trondheim, Norway 

Abstract 

This document surveys the basic concepts underlying a new geometrical frame- 
work of space and time developed elsewhere [2]. Moreover it is shown how this 
framework is used to construct a new relativistic theory of gravity. 

1 Introduction 

A proper description of Nature's fundamental forces involves their formulation as dynam- 
ical laws subject to appropriate initial conditions. This defines the initial value problem 
for whatever dynamical system one wishes to study. 

The initial value problem for gravitation generally involves what one may loosely 
call "evolution of space with time". More precisely, given a one-parameter family of 
3-dimensional manifolds St each of which is identified as "space" at "time" t (such that 
t =constant on each SeSi), the family defines an evolution of the 3-manifolds St with t. 
However, every conceivable theory of gravity includes an overlying 4-dimcnsional geomet- 
rical structure called "space-time" with which the given family must be compatible. Thus 
it is not possible to analyze a given family St in isolation; the evolution with t of fields 
that determine the relationship between the family and the overlying space-time must be 
analyzed simultaneously. When this is done one recovers the space-time geometry. 

Naturally the structure of any space-time framework should be more general than 
any particular theory of gravity formulated within that framework. This implies that in 
general there should exist geometrical relations between the space-time framework and 
any given family St that do not depend on the detailed form of the dynamical laws one 
uses. A particular set of such relations are called kinematical relations. The kinematical 
relations of different space-time frameworks are recognized by their sharing of a common 
feature; that they are not sufficient by themselves to determine solutions of the initial 
value problem. That is, in addition to the initial conditions one must supply extra fields 
coming from the dynamics at each step of the evolution of St if one wants to recover the 
space-time geometry. Thus the kinematical relations explicitly exhibit the dependence on 
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the existence of some coupling between geometry and matter fields, even if the detailed 
nature of the coupling is not specified. This is what one should expect, since any relations 
defined as kinematical always should give room for appropriate dynamics. 

As an example we take the case of the stratified non-relativistic space-time framework 
of Newton-Cartan [4] . Here there exists a unique family St consisting of Euclidean spaces 
(and left intact by Gahlean transformations), where t is Newton's absolute time. The 
overlying geometrical structure belonging to this space-time is the space-time covariant 
derivative V represented by the only non-zero connection coefficients r^j=$,j (using a 
Cartesian coordinate system {a;*}). Here $(t, x*) is a potential function and the comma 
denotes partial derivation. The equation of motion for inertial test particles is the geodesic 
equation 



+ nt = 0, (1) 



so to calculate the world lines of such test particles one must know F^^ for every member 
of the family Sf. Now we notice that there is a general relationship between the only 
non-zero components of the Riemann tensor and the potential function given by [4] 

That this relation is kinematical we see from the fact that to determine the FJ^ for all 
times and thus all information relevant for the motion of inertial test particles, one must 
explicitly supply the tensor Ritjt for each time step. Furthermore, equation (2) is valid 
regardless of the detailed coupling between Rujt and matter fields. 

As a second example we take the more realistic case of a metric relativistic space-time. 
In this case the choice of family St is not unique, so the kinematical relations should be 
valid for any such choice. Furthermore the overlying geometrical structure belonging to 
space-time is the space-time metric field, and this field induces relations between different 
choices of St. Analogous to the first example the kinematical relations are given by (using 
a general spatial coordinate system and applying Einstein's summation convention) 

Ri±j± = N'^CNnKij + Ki^Kkj + c~'^ai\j + c~'^aiaj, (3) 

see e.g. [1] for a derivation. Here a is the 4-acceleration field of observers moving 
orthogonally to each SeSt, N is their associated lapse function field and '±' denotes 
projection onto the normal direction of St- Furthermore denotes the projected Lie 
derivative in the normal direction and '|' denotes spatial covariant derivation intrinsic 
to each SeSt- The extrinsic curvature tensor K describes the curvature of each SeSt 
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relative to the overlying space-time curvature. There exists a well-known relationship 
between K and the evolution in the normal direction of the intrinsic geometry of each 
SeSt, namely (see e.g. [1], [4]) 

Kij — — — £nhij, (4) 

where £„ denotes Lie derivation in the normal direction and h denotes the intrinsic 
metric of each S&St- Another, geometrically equivalent expression of the kinematical 
relations (3) is found (see e.g. [1]) by projecting the Einstein space-time tensor field G 
into respective normal to the spatial hypcrsurfaces St. The result is {K=K'^k) 

Gij = -N-'CN^{Kij-Khij)-2Ki''Kkj + 3KK,j-^{KkiK''' + K^)hij 

- c~^(aj|j - a}'\khij) - c''^{aiaj - a^akhij) Hij, (5) 

where H and P are respectively the Einstein tensor and the Ricci scalar intrinsic to the 
hypersurfaces. The equation (4) in combination with equation (3) (or equivalently the 
totally spatial projection Gij shown in (5), the other projections representing constraints 
on initial data) may be recognized as kinematical equations. This follows from the fact 
that in addition to initial values of hij and K^j at an initial hypersurface, one needs 
to specify the tensor field Gij or equivalently the tensor field Ri±j± at each step of 
the evolution of St if one wishes to recover the space-time metric field and thus all 
information about the geometry of space-time as predicted at subsequent hypersurfaces. 
Notice that this does not depend on Einstein's equations being vahd (see e.g. [1] for a 
further discussion). 

In light of the two above examples we may call the evolution represented by a family 
St a kinematical evolution (KE) if it depends on the supplement of extra fields coming 
from an overlying geometrical structure at each step of the evolution. The question now 
is if this type of evolution exhausts all possibilities. The reason for this concern is that 
one may possibly imagine a type of evolution of a spatial hypersurface that does not 
depend on the supplement of extra fields at each step of the evolution. It would be 
natural to call such a kind of evolution non-kinematical (NKE). Since by construction 
the NKE does not give room for explicit dynamics, any geometrical framework having 
the capacity of accommodating both KE and NKE must be different from the standard 
metric framework. 
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As illustrated by the above examples we see from equations (2) and (3) that the 
extra fields one must supply at each time step to get the KE working, take the form of 
curvature tensors. Thus we may expect that to get the NKE working, it too must have 
something to do with curvature. However this curvature must somehow be available 
without referring to an overlying geometrical structure. The only possibility of realizing 
this is if the wanted curvature is obtainable as an intrinsic property of each member S of 
a unique family St- Since we want our alternative space-time framework to be relativistic, 
t must represent a "preferred" time coordinate and any member of the unique family St 
must represent a "preferred" notion of space. Observers moving normally to any S&St 
get the status of preferred observers. However, the existence of the unique family St does 
not imply that there exists a preferred coordinate frame in the sense that the outcomes of 
local experiments depend on it [2]. On the other hand it will exist a class of "preferred" 
coordinate systems especially well adapted to the geometry of space-time. One may 
expect that equations take special forms in any such a coordinate system. 

Now to the questions; are there physical phenomena which may be suspected to have 
something to do with NKE and if so, do we get any hint of which form it should take? The 
fact that our Universe seems to be compatible with a preferred time coordinate at large 
scales makes us turn to cosmology in an attempt to answer these questions. With the 
possibility of a rather dramatic rcintcrprctation of the Hubble law in mind wc guess that 
the NKE should take the form of a local increase in scale with time as measured by the 
above defined preferred observers. Moreover, analogous to the couphng between the KE 
and curvature as shown in the above examples, there should exist a connection between 
the NKE and some aspects of the intrinsic geometry of each member S of the preferred 
family St- That is, we guess that the Hubble law has nothing to do with kinematics; 
rather the Hubble law should be the basic empirical consequence of the NKE associated 
with the global curvature of space. 

The above guess represents the key to finding a new framework of space and time which 
incorporates both KE and NKE in a natural way. We call this new framework "quasi- 
metric" since it is not based on semi-Riemannian geometry and yet the new framework is 
compatible with important physical principles previously thought to hold in metric theory 
only. In this thesis we explore some of the possibilities of the quasi-metric framework; 
among other things we derive a new theory of gravity and show how this theory helps to 
define the exact relationship between curvature and the NKE. The thesis is organized as 
follows: In chapter 1 we define the quasi-metric framework and list the important results 
without going through all the gory details while chapters 2, 3, 4 and 5 are self-contained 
articles containing the detailed calculations underlying the results listed in chapter 1. 
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2 An alternative framework of space and time 



Our first problem is to arrive at a general framework describing space and time in a way 
compatible with both KE and NKE, as discussed in the introduction. Since at least parts 
of the KE of any member S of the preferred family St should be described by a Lorentzian 
space-time metric g we may assume that the KE represented by g is discernible at each 
time step. That is, at least in a (infinitesimally) small time interval centered at the 
hypersurface t =constant it should be possible to evolve this hypersurface purely kinc- 
matically in the hypersurface-orthogonal direction with the help of the 4-metric g without 
explicitly taking the NKE into account. But since this should be valid for every SESt, to 
be compatible with the NKE different SESt should be associated with different metrics 
gGgt- Thus the preferred family St of 3-dimensional spatial manifolds should define a 
one-parameter family of 4-dimensional Lorentzian metrics gj, the domain of validity of 
each being exactly the hypersurface t ^constant, but such that one may extrapolate each 
metric at least in a (infinitesimally) small interval centered at this hypersurface. 

We may think of the parameter t as representing an extra time dimension, since 
fixing t is equivalent to fixing a space-time metric gGgt on (at least part of) a space- 
time 4-manifold. This indicates that the mathematically precise definition of the quasi- 
metric framework should involve a 5-dimensional product manifold AlxR where M. is 
a Lorentzian space-time manifold and R is the real line. We may then interpret t as 
a global coordinate on R and the preferred family St as representing a foliation into 
spatial 3- manifolds of a 4-dimensional submanifold of AlxR. jV is thus by definition 
a 4-dimensional space-time manifold equipped with a one-parameter family of Lorentzian 
metrics gt in terms of the global time function t. This is the mathematical definition of 
the quasi-metric space-time framework we want to use as the basis for constructing a 
new relativistic theory of gravity. 

Henceforth we refer to the members of St as the fundamental hypersurfaces (FHSs). 
Observers always moving orthogonally to the FHSs are called fundamental observers 
(FOs). Note that the FOs are in general non-inertial. When doing calculations it is often 
necessary to define useful coordinate systems. A particular class of coordinate systems 
especially well adapted to the existence of the FHSs is the set of global time coordinate 
systems (GTCSs). A coordinate system {x^} is a GTCS if and only if the relationship 
in A/" between the time coordinate and the global time function t is explicitly given 
by the equation = ct. Since the local hypersurface-orthogonal direction in Ai by 
definition is physically equivalent to the direction along the world line of the local FO in 
A/", the one-parameter family of metrics is most conveniently expressed in a GTCS. That 
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is, in a GTCS the spatial line elements of the FHSs evolve along the world lines of the 
FOs in M] there is no need to worry about the explicit evolution of the spatial geometry 
along alternative world lines. Note that there exist infinitely many GTCSs due to the 
particular structure of quasi-metric space-time but that in general it is not possible to 
find one where the FOs are at rest. However, in particular cases further simplification 
is possible if we can find a comoving coordinate system. A comoving coordinate system 
is by definition a coordinate system where the spatial coordinates are constants along 
the world lines of the FOs. A comoving coordinate system which is also a GTCS we 
conventionally call a hypersurf ace- orthogonal coordinate system (HOGS). Expressed in 
a HOGS the 3- vector piece goi of each 4-metric vanishes on its domain FHS. Thus the 
proper time along the world lines of the FOs in J\f may be found directly from the scalar 
piece goo of each 4-metric in Ai if each 4-mctric is expressed in a HOGS. But in general 
it is not possible to find a HOGS; that is possible for particular cases only. 

The condition that the global time function should exist everywhere implies that in 
Ai, each FHS must be a Gauchy hypersurface for every allowable metric gEgf Further- 
more each FHS must be a member of a family of hypersurfaces foliating Ai equipped 
with any allowable metric gEgt when this metric is extrapolated off the FHS. Note that 
the hypothetical observers moving orthogonal to this family of hypersurfaces may be 
thought of as FOs with the explicit NKE "turned off' and that the normal curves of 
these hypersurfaces are in general not geodesies of gGg*. 

Actually it may not be possible to ehminate all the effects of the NKE just by holding 
t constant. The reason for this is the possibility that some of the NKE is not "realized" . 
That is, it is possible that the FOs move in a way such that some of the "expansion" 
representing the NKE is exactly cancelled out. In other words, the FOs may be thought 
of as "falling" with a 3-velocity — Vt and simultaneously expanding with a 3-velocity 
due to some of the NKE, the net result being "stationariness" . But this type of sta- 
tionariness is not equivalent to being stationary with no NKE; for example one typically 
gets time dilation factors associated with the extra motion. Thus some non-kinematical 
features may be disguised as kinematical features belonging to the metric structure of 
each member of the family g^. This represents a problem since the family gj should be 
found from appropriate field equations. But due to the implicit effects of the NKE, gj 
cannot represent a solution of these field equations. To overcome this problem we intro- 
duce another metric family g^ on a corresponding product manifold AlxR and consider 
a 4- dimensional submanif old jV" just as for the family g^. The crucial difference is that the 
family g^ by definition represents a solution of appropriate field equations. Thus when 
we construct such field equations we will study the general properties of gt rather than 
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those of gf Of course it must be possible to construct from gt and the vector field 
family v^, see [2] for more details. 

The family of hypersurfaces foliating Al (or A^) equipped with any allowable metric 
gegt (or gegt) defines a one-parameter family of 3-manifolds in terms of a suitably 
defined time coordinate function A sufficient condition to ensure the existence and 
uniqueness of the global time coordinate function on M. (or M.) such that the gradient 
field of this coordinate is everywhere timelike, is to demand that each FHS is compact 
(without boundaries) [3]. Furthermore, since the NKE of the FHSs by hypothesis should 
depend on aspects of their intrinsic curvature, it seems reasonable that the requirement of 
compact FHSs implies some sort of prior 3-geometry on the FHSs. That is, we guess that 
any global curvature associated with the compactness of the FHSs should not depend on 
the dynamics of gf Thus, in order to construct gravitational field equations we must 
investigate which restrictions should be laid on gt to ensure that any prior 3-gcometry 
associated with the compact FHSs does not interfere with the dynamics of gt. It turns 
out that a sufficient condition for this is that the dynamics of gt is represented by one 
scalar field. In fact, this condition means that in a GTCS we may write the general form 
of the family g^ on the form [2] 

__ _ 4. _ i2_ 

dst = \n,{x^, t)N'{x^) - N^ixi", t)] {dxy + 2-Nidx'dx° + -^N^Sikix'')dx'dx\ (6) 

where Sikdx'^dx^ is the metric of the 3-sphere (with radius cto) and to is an arbitrary 
epoch setting the scale of the spatial coordinates. Furthermore Nt is the lapse function 
family of the FOs in g^ and ^Nidx"^ is the family of shift covector fields of the FOs in g^. 

Compared to the metric framework, the restriction that each metric must be of the 
particular form compatible with equation (6) represents a reduction in the possible ex- 
plicit dynamical degrees of freedom for any theory of gravity formulated within the quasi- 
metric framework. Besides, within the metric framework any such restriction automati- 
cally implies the existence of "prior geometry", i.e., any dynamically relevant aspect of 
the geometry that cannot be changed by changing gravitating source terms. However in 
the quasi-metric framework the form (6) of gt merely indicates that an explicit dynamical 
degree of freedom has been abandoned in order to accommodate the NKE and that this 
is valid for all theories of gravity consistent with the quasi-metric framework. This is 
reahzed by the fact that in quasi-metric theory there exists a "basic" 3-geometry of 77 in 
the sense that the geometry hehj of each FHS is conformal to this basic 3-geometry but 
such that this does not interfere with the dynamically relevant aspects of the space-time 
geometry. Furthermore we may expect that the abandoned dynamical degree of freedom 
will be put back implicitly via the NKE when constructing gt from gt. It turns out that 
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this will work and that we get a consistent theory if we choose the basic 3-geometry to 
be S^, as we have done in equation (6). 

Next an obvious question is which sort of affine structure is induced on H by the 
family of metrics gt (and analogously which affine structure is induced on U by gj). To 
answer this question it is convenient to perceive gt as a single degenerate 5-dimensional 
metric on AlxR and construct a torsion- free connection V (almost) compatible with 
this degenerate metric. But in addition to metric-compatibility we should also require 
that the unit normal vector field family of the FHSs does not change when it is parallel- 
transported in the t-direction. That is, we should require that V^rij vanishes. When a 

at 

connection satisfying these requirements is found we can restrict it to M and we have the 
wanted affine structure. 

In order to construct V on x R, the metric-preserving condition 
d * * 

-Ki^tiYt, zt) = gt( Va Yt, zt) + gt(yt,Va zt), (7) 

ot St at 

involving arbitrary families of vector fields and in may be used to find a candidate 
connection where the connection coefficients are determined from g^ alone. It turns out 
that such a candidate connection is unique and and that its difference from the usual Levi- 
Civita connection is determined from the only non-zero connection coefficients containing 
t. These are given by f^^ and are equal to ^9ft)'§i9it)iicr- (Here we have introduced the 
coordinate notation g(t)ij.a for the family g^.) But since the explicit ^-dependence of gj 
is equal to that of gt shown in equation (6) (by construction N does not depend on t), 
this candidate connection does in general not satisfy the criterion Va.nt = 0. But for the 

dt 

particular case when the shift vector field vanishes it does. We may thus choose the 

connection on the grounds that it should reproduce the results of the above-mentioned 
candidate connection for the case when N'' vanishes. We then find the unique connection 
coefficients 

= ]sr6i, f° ^ (8) 

valid in a GTCS. All the other connection coefficients than (8) are found by requiring 
that they should be identical to those found from the family of Levi-Civita connections 
explicitly defined by the g^. It may be readily shown that equation (8), together with an 
explicit t-dependence of gt similar to that of gt shown in equation (6) yield 

Vagt = 0, Vant = 0, Vaht = 0, (9) 

at dt dt 

* * r 

thus the connection V has the desired properties. The restriction of V to is trivial 
inasmuch as the same formulae are valid in jV as in xR, the only difference being the 
restriction x° = in a GTCS. 
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Now we want to use the above defined affine structure on H to find equations of motion 
for test particles. Let A be an affine parameter along the world line of an arbitrary test 
particle. (In addition to the affine parameter A, i is also a (non-affine) parameter along 
any non-spacelike curve in J\f.) If we define coordinate vector fields + ^gfs" 

is the coordinate representation of the tangent vector field ^ along the curve. We then 
define the covariant derivative along the curve as 

* dt* dx"* 

V^=— Va + ^V^. 10 

9^ dX 8t dX 8^°' 

A particularly important family of vector fields is the 4- velocity tangent vector field family 
Ut of a curve. That is, by definition we have 

_ a ^ _dx" d 
"*=''W9^=rf^9^' ^^^^ 
where is the proper time as measured along the curve. 

The equations of motion are found by calculating the covariant derivative of 4- velocity 
tangent vectors along themselves using the connection in J\f. According to the above, 
this is equivalent to calculating Vut along Using the coordinate representation of 
^ in an arbitrary coordinate system we may thus define the vector field a by 
* * / dt * dx" ^ \ dt * 

a=V_9_Ut = — + ^— V_§_ )ut=—\/a_Ut + a*. (12) 

9-rt ^drt 8t drt 9x0'/ st 

Wc call this vector field the "degenerate" 4-acceleration. It may be shown that a is 
orthogonal to u^. 

The coordinate expression for a shows that this yields equations of motion, namely 

d^x"" , ff^ dt dxf^-.dx'' _ (dTf.2^a 

'd)^ + y^'^-dx + ^w'^'^ dA ^ 'dx-^dx)'' ■ ^^^^ 

To get these equations to work the degenerate 4-acceleration must be specified indepen- 
dently for each time step; thus they are kinematical equations. But a cannot be chosen 
freely. That is, given an initial tangent vector of the world line of an arbitrary test 
particle, as we can see from equation (12) a is uniquely determined by a^ and -^V ±'^t- 

However, in a subsequent section we will see that V^u^ vanishes so we in fact have a=at. 

dt 

Now a^ has the usual interpretation as an expression for the inertial forces experienced 
by the test particle. The standard coordinate expression for a^ is 

d^x" ,-r.a dxPdx'' _(dTtVi 

However, to find the curves of test particles equations (13) will be solved rather than 
equations (14); the latter are insufficient for this task since they do not take into account 
the dependence of Ut on i.e., they hold only in M.. 
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3 The role of gravity 



We are now able to set up the postulates which must be satisfied for every theory of 
gravity compatible with the quasi-metric framework. These are: (postulates for the 
metric framework [5] are stated in parenthesis for comparison) 

• Space-time is equipped with a metric family gt as described in the previous section, 
(metric theory: space-time is equipped with one single Lorentzian metric), 

• Inertial test particles follow curves for which = calculated from (13), (metric 
theory: the curves are geodesies calculated from (14)), 

• In the local inertial frames of gt the non-gravitational physics is as in special relativ- 
ity. (Metric theory: in the local Lorentz frames of the metric the non-gravitational 
physics is as in special relativity.) 

The third postulate means that we may always find a local coordinate system such that 
the connection coefficients in equation (13) vanish. Later wc will sec that in fact a=at, 
meaning that inertial test particles follow geodesies of V, see equations (48), (49) and 
the subsequent discussion for verification. Thus the local inertial frames of can be 
identified with local Lorentz frames, just as in metric theory. This means that there 
can be no local consequences of curved space-time on non-gravitational physics. But it 
is important to note that this applies strictly locally. That is, at any finite scale the 
geodesies of V do not coincide with the geodesies of any single space-time metric as long 
as there is a dependence on t. This means that the non- gravitational physics at any 
finite scale may vary in a way incompatible with any metric geometry. Despite this any 
allowable theory does possess local Lorentz invariance (LLI), and furthermore it follows 
that we do have local positional invariance (LPI), but that deviations from metric theory 
should be detectable at any finite scale. (LPI is by definition that the outcome of any local 
nongravitational test experiment is independent of where or when it is performed.) Note 
that quasi-local deviations from LLI and LPI typically occur explicitly in differential laws 
such as conservation laws. That is, such laws may be different from their counterparts in 
metric theory since non-gravitational fields may couple to fields characterizing the FHSs 
in such a way that these couplings vanish in the local inertial frames. 

The second postulate ensures that the WEP is vahd (i.e., that the free fall curves 
traced out by test particles are independent of the particles' internal structure). The 
sum of LLI, LPI and the WEP is called the Einstein equivalence principle (EEP) [5] and 
is satisfied by any metric theory. Note that the EEP was previously thought to hold for 
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metric theories only. Nevertheless, despite the differences in mathematical structure, the 
EEP should also hold for the class of quasi-metric theories according to the arguments 
made above. 

Now we want to construct a theory of gravity satisfying the above postulates. This 
particular theory should contain no extra independent gravitational fields other than 
and also no prior 4-geometry. Since the quasi-metric framework fulfills LLI and LPI we 
might hope to construct a theory which embodies the strong equivalence principle (SEP) 
(that the EEP is valid for self-gravitating bodies and local gravitational test experiments 
as well as for test particles). But as we shall see, the quasi-metric framework predicts 
scale variations in space-time between non-gravitational and gravitational systems, and 
this is not compatible with the SEP in its most stringent interpretation. So we cannot 
ask for more than that the gravitational weak equivalence principle (GWEP) is satisfied 
(that the WEP is valid for self-gravitating bodies as well as for test particles). However, if 
the GWEP is not violated, this is sufficient to eliminate the most serious effects resulting 
from the violation of the SEP. 

One distinct feature of our framework is the existence of a unique set of FHSs. This 
means that we can split up the metrics into space and time in a preferred way; i.e., we 
have a preferred family of unit normal vector fields n^. In addition we have the 3-metrics 
hi intrinsic to the FHSs. We may then define a family of 4-acceleration fields aj^ by 



where ';' denotes covariant derivation defined by the metric connection V compatible 
with any single metric gGgf. 

It is useful to denote a scalar product of any space-time object with —fit by the symbol 
'X'; this operation defines a normal projection. Analogously any projection into a FHS 
may be performed by taking scalar products with hj. To avoid confusion about whether 
we are dealing with a space-time object or a space object intrinsic to the FHSs we will 
label space objects with a "hat" if necessary. That is, the space objects At arc the result 
of a maximal number of spatial projections of the space-time object At with hf 

By repeated projections with — or any space-time object eventually reduces to 
a space object. For example the space-time vector field a^^ has the normal projection 
ajrj_= — a'^n(t)a and the components of the horizontal projection are just a^. In particular 
one may easily show that 




(15) 



d 



(16) 



a^l = 0, 



c aj:-j 
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Furthermore we define another useful quantity 

c x, + c y*=-^ + ^, (17) 

which takes into consideration how Nt changes in the direction normal to the FHSs. The 
reason why the left hand side of (17) consists of two terms is that the first term determines 
the effect on Nt through the KE of the FHSs and the second term determines the effect 
on Nt through the NKE of the FHSs. We also note that the family of extrinsic curvature 
tensors corresponding to the metric family of the type given in equation (6) has the 
form 

Km = - jc-Wk^)hm + ^(^^b- + ^^k)- (18) 

Now we want to postulate what the nature of gravity should be in our theory. To arrive at 
that, we notice that the lapse function Nt in (6) also plays the role of a spatial scale factor. 
This means that we may think of N^ as a conformal factor in (6). We may then interpret 
Nt as resulting from a position-dependent units transformation. This indicates that the 
scale factor Ft=Ntct of the FHSs may be thought of as representing a difference in scale 
between atomic and gravitational systems. That is, Ft may be thought of as representing 
a gravitational scale measured in atomic units. An equivalent way of interpreting Ft 
is to postulate that fixed operationally defined "atomic" units are considered formally 
variable throughout space-time. Since the intercomparison of operationally defined units 
only has meaning locally for objects at rest with respect to each other, any non-local 
intercomparison is purely a matter of definition [7]. 

The variation in scale between gravitational and atomic systems means that grav- 
itational quantities get an extra formal variation when measured in atomic units (and 
vice versa); in particular the dimensional constants of nature may depend on this formal 
variation. That is, the "constants" of nature may in principle depend on whether they 
are expressed in atomic or gravitational units. We notice that when performing a units 
transformation we are free to define as constants physical constants that cannot be re- 
lated to each other by dimensionless numbers [6]. We exploit this fact by defining c and 
Planck's constant h not to be formally variable. This means that measured in atomic 
units, the variation of gravitational lengths and times are identical and inverse to that of 
gravitational energy (or mass). 

The requirement that c and h are not formally variable yields no restrictions on any 
possible relationship between the variation of units and physical fields. On the other hand 
we may combine the unit charge e with h and c to form the dimensionless number /he. 
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Thus by defining e not to be formally variable we deny any possible direct connection 
between gravitation and electromagnetism. 

Formally variable atomic units represent changes of scale between non-gravitational 
and gravitational systems, and such changes of scale should be described by the metric 
family g^. Moreover should be found from field equations. We now define the scalar 
field "^t'l by definition tells how atomic time units vary in space-time. Any gravitational 
quantity is associated with to some power describing how the quantity changes due 
to the difference in scale between gravitational and atomic systems. In particular we 
find from dimensional analysis that Newton's gravitational coupling parameter Gt has 
the dimension of time units squared. Since Gt is a gravitational quantity, measured in 
atomic units it will then vary as a factor But since Gt usually occurs in combination 
with masses it is convenient to define Gt to take a constant value G and rather separate 
between active mass rrit which is measured dynamically and passive mass m (i.e., passive 
gravitational mass or inertial mass). This means that any active mass rrit is not a constant 
but a scalar field. That is, unlike c, h and e, rrit gets a formal variation when measured in 
atomic units. Since the formal variability of G should be included into rrit we find that the 
formal variation of active mass goes as Another example of a gravitational quantity 
is the active stress-energy tensor which plays the role as a source of gravitation. One 
may show that gets an extra formal variability as \1'^. The extra variability must be 
taken care of in relevant differential laws, e.g. in local conservation laws. 

The difference between active mass and passive mass represents a violation of the 
SEP. To see if this is ruled out from experiments testing the constancy of C, one may 
consider laboratory experiments (e.g. Cavendish experiments) and space experiments in 
the solar system (e.g. lunar laser ranging). In the latter class of experiments it is assumed 
that any variation in G can be found simply by replacing the gravitational constant with 
an independently variable G in Newton's equations of motion [5]. But this cannot be 
done within the quasi-metric framework since it would be inconsistent (see equation (19) 
below). Moreover, the global NKE of gravitationally bound systems cannot be neglected 
when testing the temporal constancy of G. (For an exphcit example of this, see [10].) 
Similarly, in laboratory experiments a variable G should not simply be put into Newton's 
force law without taking into account the global NKE. Due to the global effects of the 
NKE one must also separate between bodies bound by atomic forces and gravitationally 
bound bodies when calculating the time variation of G. In particular, when testing 
whether or not G changes with time, according to quasi-metric theory the result should 
depend on the specific experiment performed. For example in Cavendish experiments, 
where the experimental system is bound by atomic forces, one should expect a different 
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result than by trying to measure the variation with time in the acceleration of gravity 
on the surface of the Earth, say. (The former experiment is expected to show that G 
increases with time whereas the latter experiment is expected to show that G decreases 
with time.) Now laboratory experiments testing the constancy of G are nowhere near 
the precision level attained from the corresponding space experiments. But since the 
interpretation of the space experiments is explicitly theory-dependent, they are not very 
useful as direct tests on the variability of G as predicted by a quasi-metric theory. And 
laboratory experiments are not yet sufficiently precise to rule out the predicted variations 



We conclude that any difference between active and passive mass-energy is a violation 
of the SEP since there are locally measurable consequences of this. But the fact that the 
locally measured G is variable only via the variation of (atomic) units (in contrast to any 
"physical" variation in G independent of the variation of units) is important regarding 
the interpretation of experiments. So far it seems that no experiment has been capable 
of ruling out the formal variation of G. 

The following postulate is consistent with the above considerations: 

• Gravity is described by the metric family (constructible from another metric 
family g^) measured in atomic units. Measured in atomic units, gravitational quan- 
tities get an extra formal variability depending on dimension. By construction c, 
h and e are not formally variable. But to be able to define G as constant it is 
necessary to separate between active and passive mass-energy. In particular the 
formal variability of active mass-energy is included into the active stress-energy 
tensor considered as a source of gravity. The formal variation of atomic units in 
space-time is by definition identical to the variation of the scalar field How '^t 
varies in space-time depends on the fields Oj^j, ^^j^-, and 

To get any further we must know the formal dependence of on aj^j, and |. 

Prom equation (6) we note that the scale factor Ntct of the PHSs is measured in atomic 
units. This suggests that any change in this scale factor may be explained solely by 
means of the formal variation of units. We thus define ^t=^-j: so we have 



Note that the Newtonian approximation of equation (6) (which is obtained by neglecting 
the shift covector, setting the spatial scale factor equal to unity and replacing Sij with 



of G. 




mt,j = c "^arj-mt. 



(19) 
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an Euclidean 3- metric), is inconsistent with any variation in space-time of active mass. 
That is, in quasi-metric relativity G must be a constant in the Newtonian limit. 

Then there is the question if the FHSs play any explicit role in the couphng of matter 
to geometry. And indeed they do, but only via Kj and derivatives of Ft as we shall 
see below. The curvature intrinsic to the FHSs should not couple directly to matter 
fields. Before we try to find gravitational field equations we first find how the quasi- 
metric local conservation laws differ from their counterparts in metric theory. That is, 
the degenerate covariant divergence V -T^ should not be expected to vanish in general 
due to the postulates determining the role of gravity within the quasi-metric framework. 

It may seem natural to believe that local conservation of mass-energy and momentum 
is necessarily violated since V -Tj may be non-zero. However, given the possibilities of 
coupling Tt to the fields a^rj, ^^^^j ^jf^ and ^ such that LLI and LPI hold, it is in fact 
possible that local conservation laws within the quasi-metric framework do indeed take 
the form V -T^t^O, but such that their special relativistic form still holds in the local 
inertial frames. Note that the formal variation of the units is taken care of both in 
derivatives of nit and in derivatives such as V-T^. 

Proceeding with the task of finding local conservation laws; to separate out most 
effects coming from the conformal factor in (6) we define the metric family g(=iVi~^gi 
and its associated degenerate covariant derivative V • In coordinate notation we denote 
the operations of degenerate covariant derivation V and V with the symbols and 

respectively. Using this notation, a short calculation yields the following expression, 
which turns out to be helpful: 

{t)tJ,*u=(^ ^ {t)n*t + ^ {t)iJ.;u — ^ {t)iJ,*u + ~/V^ Wi" 

We now project this expression into respective normal to the FHSs. It is possible to 
figure out expressions for the projections of T^t)iJ,*v guessing how they should behave 
for a null fluid. These calculations can be found in [2] , where it is shown that the local 
conservation laws in 77 read (note that these laws take into account the effects of the 
formally variable units) 

T^t)L.. = 2c-^a^4)i - 2(^ + {cNt)-'[\ + ^])T^t)ll. (21) 

T(t,r*^ = 2c--arrflt,, - 2(^ + (ciV,)-^[^ + ^])T(,i, (22) 

Note that since c'^aj^j = -^t^jl^ ^ = --/Vtfjj and | + ^ = ^1^^ in a GTCS, these 
local conservation laws take their special relativistic form in the local inertial frames, as 
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asserted previously. Besides, by applying equations (21) and (22) to a source consisting 
of a perfect fluid with no pressure (i.e., dust) and projecting the resulting equations with 
the quantity gj + c~^Ut(8)Ut, one may show that the dust particles move on geodesies of 
V in AT. 

Next, the projections of (21) and (22) yield the conditions [2] (with Kt=K^^f.-^-) 

^iitT{t)ll = {Kt - 2-^)r(t)ii + K{t)ikflt) - (23) 



t iV* 

+c-Wifl,^j + c-^ar,j^T^t)U - Tltw (24) 



Equations (23) and (24) determine the evolution in the hypersurface-orthogonal direction 
of the matter density and the momentum density m M.. To find these evolutions in AT 
we add the dependence on t. Using the projected Lie derivative Cat valid in AT we find 

KTit)li = - {cN,)-'[\ + ^nt)lr (26) 

Note that one expects that the global NKE implies local non-conservation of mass-energy 
and momentum in special situations; in these cases equations (25) and (26) will be 
violated. That is, whenever the global NKE is so strong that non-gravitational quantum 
fields cannot be treated as localized to sufficient accuracy, one should get spontaneous 
pair production from excitations of vacuum fiuctuations of such quantum fields, violating 
equations (25) and (26). This will be important only in the early universe and possibly 
in extremely deep gravitational wells, though. On the other hand, equations (23) and 
(24) are required to hold always. 

Since Tj is the active stress-energy tensor it is not directly measurable locally. That 
is, one must know how it relates to the passive stress-energy tensor % in H. Now the 
relationship between and % depends in principle on the general nature of the matter 
source. To illustrate this, consider and % for a perfect fluid: 

Tt = {Pm + C-2p)Ut0Ut + pgt, %={pm + C-2p)UtOUt + pgt, (27) 

where p^, p are the active mass-energy density and the pressure in the local rest frame of 
the source. Here we have used the deflnitions pjn\/ht=PmVht and pJht='p\/ht, where pm 
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and p are the passive counterparts to pm and p, respectively. Also ht and ht are defined 
to be the determinants of and h^, respectively. Furthermore the relationship between 
Pm and pm is given by 

{j-NtPm for a fluid of material point particles, 
I - (28) 
jiNfpm for a null fluid, 

and a similar relationship exists between p and p. The reason why the relationship 
between and pm is different for a null fluid than for other perfect fluids is that spectral 
shifts of null particles influence their passive mass-energy but not their active mass-energy. 
Besides, note that more "physically correct" local conservation laws can be found by 
calculating V-Tt when gt is known. But these local conservation laws for passive mass- 
energy do not take any predetermined form. 

To flnd fleld equations, the requirement that we must have metric correspondence 
with Newtonian theory for weak stationary flelds leads to the postulate {k='&t^G/c^) 

2i?(t)ll = 2(c-2a!,|, + c-^-ar^a]p - Km^) + ^^^Kt) = + (29) 

where Rt is the Ricci tensor family corresponding to the metric family (6) and is the 
extrinsic curvature tensor family of the FHSs in A4. Equation (29) represents one fleld 
equation with the particular property that it does not depend explicitly on the intrinsic 
curvature of the FHSs. We postulate a second set of fleld equations having this property, 
namely 

Note that metric counterparts to equations (29) and (30) hold in general relativity (GR) 
for projections on an arbitrary spatial hypersurface. 

In particular cases it is possible to choose a convenient coordinate condition which 
makes the field equations more easy to handle. That is, the fleld equations become less 
complicated if 

N%^3c-'ariN^ => ^t = ^~. (31) 

Note that the geometrical interpretation of the condition (31) is that the extrinsic cur- 
vature family Kt (calculated from the metric family gj) is traceless. In other words the 
condition (31) imphes that the FHSs are maximally embedded in g^. 

Several comments apply to the fleld equations (29), (30). Firstly we notice that they 
do not represent a full coupling of Tj to space-time geometry but only a partial one. 
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Secondly, the field equations superficially look like just a subset of the Einstein field 
equations in ordinary GR. But the crucial difference is that the equations (29), (30) are 
valid only for projections into the FHSs; for projections into any other hypersurfaces they 
do not hold. Thirdly we observe that the FHSs enter the field equations geometrically 
only via their extrinsic curvature; the intrinsic curvature of the FHSs does not couple 
directly to matter. Fourthly equation (29) is a dynamical equation whereas equations (30) 
are constraints. We thus have one dynamical equation and three constraint equations. 
One may compare this to GR where from (5) one can see that the Einstein equations 
split up into six dynamical equations and four constraint equations. Fifthly, the local 
conservation laws do not follow from the field equations as they do in GR. Rather the local 
conservation laws are coupled to the field equations and represent independent restrictions 
on the nature of the matter source. Finally we observe that our field equations contain 
only one explicit dynamical field, namely Kt- Since equations (18) and (31) indicate 
that Kt = S^^^ at least for some cases, this means that the true dynamical field is the 
lapse function Nf. That is, our field equations have only one explicit dynamical degree of 
freedom whereas GR has two [4]. However, a second implicit dynamical degree of freedom 
is added to our theory via the NKE when constructing gt from gj [2]. 

4 Kinematics and non-kinematics of the FHSs 

In this section we take a closer look at the evolution with time of the geometry of the 
FHSs in J\f and in J\f. To begin with (6) yields the evolution equation [2] 

= 2((ciV,) [- + _]-_ + -c a^,^)/.,,,--(_ + ^). (32) 

On the other hand we may define the split-up into respectively kinematical and non- 
kinematical terms of the time evolution of the spatial scale factor in (6). We may thus 
define 

c-%^c-^HMcN,r{l^^)-% (33) 

where Hf is a quantity associated with the NKE of the FHSs. (We define in equation 
(34) below.) 

In the introduction we defined the NKE as a local increase of scale with time of the 
FHSs. To quantify this definition we introduce the local scale factor Ft=Ntct of the 
tangent spaces of any FHS in AT. By definition the local change of scale with time as 
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defined by nearby FOs is entirely expressible as a change in Ff. Thus we may define the 
non-kinematical scale rate Ht by 

Ht=cF-%^Ft - c-^xt. (34) 

Note that we have subtracted the kinematical contribution c~^Xt to the scale change 
in equation (34). Since Ht is a purely non-kinematical quantity, it should be available 
directly from the geometry intrinsic to the FHSs. We argue in [2] that the correct form 
of this so-called non-kinematical "coupling" equation is given by 

= i + c-ya^a^j. (35) 

Besides we see from (17) that yt represents the non-kinematical part of the evolution of 
Nt with time. Combining equations (17), (33) and (35) we find 

J/t = -^t - ^ = c~^\[^r^i, (36) 



,-,, = _^+(eiV.)-'(i + ^)-c-'i?,. (37) 

For cases where there is no coupling to non-gravitational fields via the equations of 
motion (13) in jV' we now have a well-defined initial- value system in jV": on an initial 
FHS we choose initial data TV*, -/V*,!, -/V*,*, -/V*, (-/V*),x and the various projections of 
which satisfy the constraint equations (30); then we evolve through jV" using equations 
(29), (30), (18), (19), (23), (24), (25) and (26). For each step of the evolution we must 
supply the quantity 2^(^^xi be able to construct and for the subsequent steps. 
Then the non-kinematical quantities follow automatically via (35) and (36). Note that 
the local conservation laws (23), (24), (25) and (26) are coupled to the field equations 
since they represent non-trivial evolution equations of T^^^y^-^ and But to know 

the time evolution of the full stress-energy tensor we need to know the equation of state 
characterizing the matter source; i.e., we need to know how the passive mass-energy is 
related to the passive pressure. 

Note that for the initial- value system defined by the equations (4) and (5) it is possible 
to choose the lapse function N freely throughout space-time since the evolution of iV is 
not determined by the constraint preservation [8]. On the other hand, for our initial- value 
system Nt is the lapse function for a particular set of observers; the evolution of the FHSs 
tells how Nt evolves with time and there is no arbitrariness of choice involved. 
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But in general the evolution scheme in JJ is coupled to the evolution in scheme H via 
non-gravitational fields. This means that these evolution schemes are intertwined, and 
that we need to construct gj for each step of the evolution. This can be done by finding 
the family of 3-vector fields — Vj which tells how much the FOs in "move" compared 
to the FOs in J\f. Expressed in a GTCS said vector field may be defined as [2] 



^lt)=yt^T^ v=ytyh(t)ikbrb'r, (38) 

where b^F is a 3-vector field representing "the local distance from the center of grav- 
ity" generalized from the spherically symmetric case. Furthermore bjF is defined by the 
equation 

o^lfe + c"^ajrfea^] b^r - + c"^a^fea^] 6^ - 2a^ = 0, (39) 

which has the unique solution bjF = r-^ for the spherically symmetric case (using a 
spherical GTCS). This means that we intuitively expect that b:^ should vanish when ajF 
does and that a^F and hj^ should have the same direction. 

From equation (38) we see that the construction of from g^ involves that part of 
the NKE which is not "reahzed" in the evolution equation (34) . That is, from equation 
(36) we see that yt — c~^^aj:ia)p, and this quantity is not included in the explicit time 
evolution of the spatial scale factor. Note that v should not depend explicitly on t. It 
is also worth noticing that there exists only one propagating degree of freedom in jV". 
However the quantity v represents a second degree of freedom so one would expect that 
the number of propagating degrees of freedom in H is two, the same as in GR. 

In [2] it is shown how to construct g* (or, equivalently h^, N and N^) algebraically 
from gt (or, equivalently h^, Nt and N^) and the vector field family Vj for each step of 
the evolution; thus we also have a well-defined initial- value system in H . 

Actually, transformations similar to gt— >gt apply to any tensor field which norm is 
required to remain unchanged when g*— >gt. As an example we list the formulae valid for 
the transformation Z^— ^Zj where Z^ is a rank 1 tensor field family. These transformation 
formulae involve the unit vector field e.i,=^e\-^ along the bjp-direction in TV" and the 
corresponding covector field e^=j^e^dx\ The transformation formulae read 



(1-^)^)0, = (40) 



2- _ -7 2- 



= + Y±^S^Me], 4) = - ^(e1Z(,))el. (41) 
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Note that it is also possible to find analogous formulae (involving the unit vector field 
66=^6^^ in the bj^-direction and its associated covector field e^=j^e\dx^ in M) for the 
inverse transformation Z^— >Zt. 

One example of a (co)vector field obeying the rules (40)-(41) is a general 4-velocity 
(co) vector family Uj. Another important example is any 4- acceleration field family 
determined from a non-gravitational force (e.g., if is due to the Lorentz force acting 
on charged matter) . For such cases the norm of a^ must be invariant under the transfor- 
mation. In particular may vanish, and this means that geodesic motion in 77 implies 
geodesic motion in M . Note that ajc- does not in general transform according to the rules 
(40)-(41) since S-jr is determined from the requirement that the FOs move normal to the 
FHSs rather from some non-gravitational force. 

For completeness we also list the transformation formulae valid for the transformation 
Qt^Qt, where is a rank 2 tensor field family. We have 

Q(t)oo = (l - <9(t)oo, <3°°) = (l - (42) 




2- 

Q{t)i3 = Q{t)ii + Jvy e^{^\Q{t)kj + Q{t)ik^^j), 

Q%) = Ql)-^^vyAi^A'^) + QtA)- (44) 

The transformation formulae shown in equations (42)-(44) may in particular be used to 
find gt from g^; moreover they may easily be generalized to tensor field families of higher 
rank. 



5 Test particle motion revisited 

Once we know g^ we may calculate paths of arbitrary test particles using the equations 
of motion (13). But before we are able to do such calculations we need to have an 
independent expression for the degenerate acceleration field a. In this section we will see 
that we in fact have a^a^. 
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In the following it is convenient to introduce the difference in coordinate 3-velocity 
between an arbitrary test particle and the local FO. We call this quantity Wt, and we 
write its square as w'^. Thus we define (using a GTCS) 



where Tt is the proper time as measured along the curve. Note that Wj is an object 
intrinsic to the FHSs and that we may separate out the t-dependence of Wj, namely [2] 

One may readily show that the tangent 4- velocity can be split up into parts respectively 
normal and tangential to the FHSs, i.e. 

Ut = 7(cnt + wt). (47) 

Since any 4-velocity vector in Af by construction has constant norm we may locally 
calculate the transformation Uj— >Ut from equations (40) and (41). We may then show 
that the norm w of the coordinate 3-velocity Wt (and thus 7) is invariant under the 
transformation, i.e., w — w. 

Using equations (46) and (47) it is now possible to calculate the quantity VjiUt and 

dt 

thereby the degenerate acceleration a from equation (12). That is, using equation (9) it 
is straightforward to show that 

Va.wt = 0, VAUt = 0, a = at. (48) 

dt at 

The coordinate expression (13) for the equations of motion then take the form (this form 
is valid in a general coordinate system, but equation (8) is valid only in a GTCS) 

^ + y^'^'dx + ^(^^^^Ix) -dx-^^)' ^^^^ 

For timelike curves one may parameterize the curve by proper time by choosing X — cTf. 
For null curves A must be chosen as some other affine parameter such that ^ = 0. 

The equations of motions (49) are the geodesic equations obtained from V. Moreover, 
it can be readily shown that (49) does not reduce to the geodesic equation (14) whenever 
the dependence on t of the metric components cannot be neglected. That is. if the global 
part of the NKE (see equation (35)) is not negligible inertial test particles do not 
move as if they were following geodesies of any single space-time metric. 
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6 Discussion 



As we have seen in this survey, the reahzation that physical space should not necessarily 
be treated as a purely kinematical field may have important theoretical consequences; 
this idea leads ultimately not merely to a new theory of gravity but to a new geometrical 
framework as the basis for relativistic physics. One of the reasons that makes the pos- 
sible replacement of the currently prevailing metric framework with our alternative even 
remotely probable, is the fact that this can be done without sacrificing the EEP. 

Even though our theory seems well-posed and self consistent, the question is still 
open if it is viable according to experimental criteria. We show in [2] that it is at 
least possible to find certain limits where the predictions of our theory match those of 
GR within experimental precision. Particular cases here are static systems with spherical 
symmetry; we show that by assuming that the effects of a global scale factor are negligible, 
it is possible to find a one-parameter family of identical, conformally flat metrics as 
a vacuum solution of our fleld equations. We then include the local effects of the NKE 
into the family gt which we construct from the conformally flat metrics. That is, when 
the global effects of the NKE can be neglected the metric famihes gt and gt can be 
represented by single metrics g and g since there is no explicit dependence on t in this case. 
Moreover, for the static, spherically symmetric case it turns out that g agrees with the 
Schwarzschild metric to sufficient order in the quantity r^/r, where is the Schwarzschild 
radius. However, there is disagreement for higher orders since neither g nor g possesses 
an event horizon. This illustrates the fact that the geometrical structure of quasi-metric 
space-time does not allow the formation of black holes. Furthermore we easily flnd a toy 
cosmological solution of our fleld equations; a family of metrics with the geometry S^xR 
constitutes a vacuum solution of the fleld equations. (In this case the families gt and gt 
are identical.) This indicates that our theory does not allow any physical singularities 
since even at the cosmological singularity no physical quantities diverge. (But we need a 
matter creation mechanism for a Universe with an empty beginning, this may perhaps be 
constructed along the lines described in [9].) Moreover no coimterparts to the Friedmann- 
Robertson- Walker models exist in quasi-metric theory; besides the empty toy model the 
only allowed isotropic cosmological solutions of the field equations represent toy universes 
filled with a null fiuid. This means that the usual cosmological parameters are not useful 
in quasi-metric theory. 

Another prediction of our theory is that the use of Newtonian dynamics for weak 
field systems has its limitations; in addition to the weak field criterion systems should be 
"small" so that the global eSects of the NKE (i.e., the effects due to a "real" expansion) can 
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be neglected. (That is, even if these global effects are small locally, over large distances 
or long time spans they may not be small.) We explore some of the global effects of the 
NKE in [10], where we discuss the possible relevance of the results for explaining the 
observed asymptotically non-Keplerian rotation curves of spiral galaxies. Furthermore, 
even in the solar system there may be observable consequences of the NKE. In [11] we 
discuss one such effect and compare to data the predictions of the model considered in 
[10]. In addition we show in [12] that the global effects of the NKE are predicted to 
apply inside gravitationally bound perfect fluid bodies as well. In particular quasi-metric 
theory predicts that the radii of bodies made of ideal gas should not be constant in time 
but rather change similarly to the cosmological scale factor. This model can be compared 
to data on expansion of the Earth. Also possible is comparison of the model to data on 
spin-down of the Earth's rotation rate. 

A lot of work remains to be done before one can decide whether or not our theory 
is viable. For example it is necessary to work out if it correctly models binary pulsars, 
quasi-stellar objects, cosmological observations etc. Should our theory turn out to be 
viable, it is believable that its possible acceptance hinges on if it is able to explain the 
workings of the Universe from first principles better than competing theories. Only then 
will the NKE get the status as a fundamental principle realized by nature. 
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Abstract 

We introduce the "quasi-metric" framework consisting of a 4-dimensional dif- 
ferential space-time manifold equipped with a one-parameter family of Lorentzian 
4-metrics parameterized by a global time function. Within this framework, we de- 
fine the concept "non-kinematical evolution of the 3-geometry" and show how this 
is compatible with a new relativistic theory of gravity. It is not yet clear whether 
or not this theory is viable. However, it does agree with the classic solar system 
tests. 

1 Introduction 

The class of space-time theories defined as metric theories of gravity satisfies the following 
postulates [1]: (1) space-time is equipped with a Lorentzian non-degenerate metric, (2) 
inertial test particles follow non-spacelike geodesies of this metric, and (3) in the local 
Lorentz frames of the metric, the non-gravitational physics is as in special relativity. 
In light of the highly successful record of the metric framework, and with the lack of 
any obvious experimental motivation, one might wonder why anyone should bother to 
construct a non-metric alternative. The answer to this question is that even though it is 
possible to model a wide range of phenomena within the metric framework, things do not 
always follow from first principles. If one within a non-metric framework can do at least 
as well as with any metric theory, and if additionally one gets a better connection to first 
principles, the non-metric framework should be considered as a serious alternative. The 
key questions are, of course, firstly; which physical/theoretical principles should guide one 
towards constructing a non-metric framework, and secondly; which phenomena should be 
considered being sufficiently fundamental to motivate the formulation of such principles. 

In the search for a suitable non-metric framework part of our phenomenological mo- 
tivation consists of the fact that the observable large-scale Universe seems to take part 
in a global increase of scale. Even though this observation has a natural interpretation 
within the metric framework, it does not follow from first principles. That is, within the 
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metric framework one may model many types of universe besides the expanding ones, e.g. 
collapsing universes, static universes or universes whose evolution at large scales cannot 
be modelled by a global scale factor. Contrary to this we want to construct a non-metric 
framework which explains the above-mentioned observation from first principles. 

The above considerations suggest that a central feature of our non-metric framework 
should be the postulated existence of a specific set of (hypothetical) observers whose 
interrelationship consists at least partly of an increase of local scale with time. For 
reasons explained in section 2 this increase of local scale with time is called the "non- 
kinematical evolution" (NKE) of the 3-geometry spanned by these observers. In this 
paper we construct a particular type of non-metric framework which accommodates this 
concept as part of the basic geometrical structure. 

The introduction of the NKE can be done without sacrificing one of the basic principles 
defining metric theory; namely the Einstein equivalence principle (EEP). However, the 
mathematical structure of the new framework is not based on semi-Riemannian geometry, 
so postulate (1) above does not apply. But despite that the new framework is so similar 
to the metric framework that we may call it "quasi- metric". 

This paper is organized as follows: In section 2 we define the quasi-metric framework in 
general terms, in section 3 we construct a theory of gravity consistent with this framework, 
and in section 4 we find special solutions and discuss the question of testing. A list of 
acronyms and mathematical symbols can be found in an appendix. 

2 The quasi-metric framework 

We start by briefly describing the well-known concept of "kinematical evolution" (KE) 
of the 3-geometry, also called geometrokinematics, within the metric framework. 

By slicing an arbitrary 3-dimensional (Cauchy) hypersurface S out of the space-time 
manifold M. such that its normal vector field is everywhere timehke with respect to the 
metric g, one splits space-time into space and time. Suppose one furthermore specifies 
an arbitrary "deformation" timelike vector field on S. Then, from the intrinsic metric 
of S, the extrinsic curvature tensor field and its rate of change along the deformation 
vector field, one may deduce the metric on the deformed spatial hypersurface. This 
is the essence of geometrokinematics, and its detailed relations hold irrespective of any 
dynamical laws. The reason is that these relations describe how a single space-time tensor 
field, sufficiently represented by the Einstein tensor field G, is projected into tangential 
and normal directions to S, and this severely restricts how the projections behave under 
deformations of S, see e.g. [2]. 
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The kinematical relations of Ai have one characteristic defining property, namely that 
they are not sufficient by themselves to determine the space-time geometry. That is, in 
addition to the initial data on S one must supply extra fields coming from the overlying 
space-time geometry at each step of the deformation of S if one wants to predict the 
metric on S for each time step of the deformation process. This is as it should be, since 
the geometry of space-time should be dependent on the coupling to geometry of matter 
fields; only when appropriate dynamics is specified should it become possible to calculate 
the space-time geometry. 

Any equations recognized as kinematical must leave room for dynamics without spec- 
ifying the exact nature of the coupling between matter and geometry. This is realized 
geometrically by the fact that the extra fields that must be supplied take the form of 
curvature tensors. As noted above it is sufficient to supply the projection of G into S 
for each step of the deformation; in itself this does not depend on how G is coupled to 
matter. See [2] for a detailed discussion. 

The kinematical relations hold equally well for the KE of any spatial 3-hyper surf ace. 
By contrast, one might think of restricting the KE of some particular spatial hypersurface 
to take some prescribed form along some particular deformation vector field. This would 
be equivalent to just laying restrictions on the form of the space-time metric, and the KE 
of the special hypersurface up the prescribed foliation would just represent the kinematical 
relations applied to this particular foliation. 

However, one might imagine that there are aspects of 3-geometry evolution which do 
not depend on the continuous supplement of curvature fields coming from an overlying 
space-time metric. Rather, it would seem conceivable to try to construct a type of 
evolution which may be defined in terms of the intrinsic curvature of S. Such aspects 
of 3-geometry evolution one may call "non-kinematical" since they by construction are 
self- advancing, i.e., they leave no room for dynamics. It is clearly not possible to make 
any sense of the NKE concept within the metric framework due to the global restrictions 
of the metrical field. However, one may try to get around this by postulating more than 
one space-time metric field; thus one may define some sort of evolution of the space-time 
metric field. Then there is no reason to expect that the NKE should have anything to 
do with the kinematical relations defined by any of the single metric fields. Of course 
it might be possible to construct a single metric field containing a particular foliation 
such that the evolution of S up this foliation is indistinguishable from the evolution of S 
resulting from an evolution of space-time metric fields. But the evolution of S should be 
interpreted as physically different for these two cases; in principle it should be possible 
to perform test experiments that are capable of showing whether or not the evolution of 
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S is purely kinematical. 

We want to create a non-metric framework which incorporates the NKE in a natural 
way. Since this framework should include more than one space-time metric field and some 
way to evolve one of those into another, there must exist some evolution parameter by 
means of which this can be defined. This parameter must have an existence independent 
of the details of each single metric field. The straightforward way to define such a 
parameter is to equip the space-time manifold with a global time function; this quantity 
should be a natural feature of the general geometrical structure on which we base our 
non-metric framework. We proceed by defining such a basic geometrical structure. 

Consider a 5-dimensional differential manifold with the global structure AixH, where 

is a Lorentzian space-time manifold and R is the real line. We may then introduce 
local coordinates {x"^} on Ai (e taking the values 0,1,2,3 where is the time coordinate), 
and a global coordinate t on R. Suppose next that we slice a 4-dimensional hypersurface 
A/" specified by the equation t = t(x'^) out of TVf xR. On this slice there is required to exist 
a unique invertible functional relationship between t and throughout the coordinate 
patch covered by {x*^}- Now we let a point with coordinate to in R represent a space-time 
metric gtu on Ai. Then A/" is defined to be a 4-dimensional space-time manifold equipped 
with a one-parameter family of Lorentzian 4-metrics gj parameterized by the global time 
function t. This general geometric structure is the basis of the quasi-metric framework, 
and it has the capacity to accommodate both KE and NKE. We come back to questions 
of uniqueness and existence of such a geometrical structure in the next section. 

The global time function t singles out of a set of fundamental 3-dimensional spatial 
submanifolds St- Each of these hypersurfaces is "preferred" in some sense since each 
of these hypersurfaces constitutes the domain of validity for one single member of the 
metric family; there corresponds exactly one metric field gGg* to each SeSf. We call 
the hypersurfaces St the fundamental hypersurfaces (FHSs). Observers always moving 
orthogonally to the FHSs are called fundamental observers (FOs). It is important to have 
in mind that each 4- metric applies exactly only on its domain FHS. This means that 
each member of the metric family is relevant to the evolution of the FHSs only at the 
corresponding FHS. Thus it is sufficient to extrapolate each 4-metric to an infinitesimally 
small neighbourhood of the FHSs whereas any further extrapolation may be of dubious 
value physically. 

Now the FHSs are subject to two types of evolution with proper time Tjf (depending 
directly on t) as measured by the FOs when traversing the one-parameter family of 
4-metrics. Firstly, it is the KE of the FHSs as defined from its dependence on the 
supplement of curvature information from J\f; this information should be available from 
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the 4-geometry defined by the metric family. Secondly it is the NKE of the FHSs; this 
type of evolution is defined by being explicitly dependent on some aspects of the intrinsic 
curvature of the FHSs and not on the space-time curvature calculated from the metric 
family. Prom this it follows that the concept of NKE is intimately connected with the 
existence of "preferred" hypersurfaces from which the relevant intrinsic curvature can be 
calculated. On the other hand, within the metric framework all spatial hypersurfaces are 
equivalent as far as the kinematical relations are concerned. This is true even though 
different spatial hypersurfaces in general will have very different intrinsic curvatures; that 
is why any sort of NKE is impossible in metric theory. 

If one considers a metric space-time manifold as given a priori, there are in principle 
two different ways of splitting up space-time into space and time. The first one is the 
foliation of space-time into a family of space-like hypersurfaces. This corresponds to 
choosing a time coordinate representing "surfaces of simultaneity" . Note that the gradient 
field of this time coordinate is in general not orthogonal to the hypersurfaces. The second 
way of splitting up space-time is by foliating space-time into a congruence of time-like 
curves representing a family of observers. This corresponds to assigning a set of constant 
spatial coordinates to each observer, i.e., each curve of the congruence is characterized by 
a set of three numbers. There is, however, no natural choice of simultaneity connecting 
the local rest spaces of the different observers representing the family of curves. 

In quasi-metric relativity each metric of the metric family is defined from the local 
measurements of the FOs with the global time function serving as the (unique) choice of 
simultaneity. The construction of quasi-metric space-time as consisting of two mutually 
orthogonal foliations means that it must be possible to choose a global time coordinate 
for each metric of the metric family. Consequently, it should exist a set of "preferred" 
coordinate systems which adapt to the FHSs in a natural way. A set of such coordinate 
systems is found by first noting that it exists, by definition of the global time function, 
a unique invertible functional relationship between and t throughout the coordinate 
patch covered by {x^}- We may then define a particular type of coordinate system where 
this relationship is given explicitly by x^ — ct. Such a coordinate system where x° and 
ct are identified and thus physically equivalent in J\f we call a global time coordinate 
system (GTCS). It is usually most convenient to use a GTCS since such a coordinate 
system is better adapted to the quasi-metric space-time geometry than any alternative 
coordinate system. In particular cases further simplification is possible if one can find a 
comoving coordinate system, i.e., a coordinate system where the spatial coordinates {x*} 
assigned to any FO do not vary along its world line. That is, the spatial coordinates 
along such a world line do not depend on t in J\f. A GTCS which is also a comoving 
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coordinate system we call a hypersurf ace- orthogonal coordinate system (HOCS). Notice 
that in principle infinitely many GTCSs exist, but that a HOCS can be found only in 
particular cases. 

If we define a GTCS {x^}, each single metric of the one-parameter family can be 
expressed by its components g^v{x^). The fact that we have a one-parameter family of 
metrics rather than a single metric on J\f we express by an implicit dependence on the 
coordinates via the global time function t. In component notation we write this formally 
as g{t)fj,u{x'^) ■ where the dependence on t is put in a parenthesis not to confuse it with 
the coordinate indices. By treating i as a constant, we get the components of each single 
metric. Note that the off-diagonal components g(t)oi=g{t)io (where the spatial indices i 
take the values 1,2,3) vanish in a HOCS. 

By postulating the existence of the family we automatically define a family of Levi- 
Civita connections Vf, where each connection is found from the corresponding metric by 
holding t constant. But we are more interested in the affine structure induced on J\f 
by the metric family than in the Levi-Civita connection of any single metric. Thus we 
want to define a unique covariant derivative in A/", given the functions g(t)n.v The natural 
way to do this is by thinking of the metric family as one single degenerate 5-dimensional 
metric on TWxR, where the degeneracy manifests itself via the condition gi(^, ■)=0. 
We may then introduce a torsion-free, metric-compatible 5-dimensional connection V 
on A^xR and consider the restriction of V to A/"; in this way we obtain the wanted 
affine structure on J\ . We name the connection V the "degenerate" connection since 
it is (almost) compatible with which may be perceived as a degenerate 5-dimensional 
metric. The requirement that V must be metric-preserving reads (a comma denotes 
partial derivation, and we use Einstein's summation convention throughout) 

Moreover we define the degenerate connection coefficients F^^, and F^^ by 

* d * d ** * d * d * 

— _ pa pa _ pa vy — pa pa —pa (Q Q\ 



where F"^^^^ are the connection coefficients obtained from the Levi-Civita connections 
Vt. Note that equation (2.2) implies that V does not have torsion since it is symmetric. 
This means that any counterpart to equation (2.1) obtained by e.g. interchanging the 
coordinate vector fields ^ and ^ does not hold in general. So V does not preserve the 
degenerate part of degenerate 5-dimensional metric on AlxR. Furthermore notice that 
other degenerate connection coefficients than those given by (2.2) vanish by definition. 
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This implies that the gradient of the global time function is covariantly constant, i.e., 
thectV a_dt —V a dt — 0, see [4] for an analogy with Newton-Cartan-theory. 

We must now find expressions for the degenerate connection coefficients f^^ defined 
in equation (2.2). If the metric- preserving condition (2.1) were the only restriction one 
wanted to lay on the degenerate connection, a natural choice of degenerate connection 
coefficients would be ones totally determined from partial derivatives with respect to t of 
the functions g{t)ij.v, i-e., from terms containing expressions like ■§^g(t)ij.v If can be shown 
that the requirement that the degenerate connection coefficients take this form yields 
that must be equal to \g"t)9{t)an)t- But in addition to the condition (2.1) we also 
have the condition that the geometrical structure of quasi-metric space-time should not 
be affected by parallel-transport in the t-direction. That is, we should require that the 
degenerate connection ensures that the unit normal vector field family xit (i.e., a vector 
field family having the property gf(nf, nj) = —1) of the FHSs will be parallel-transported 
in the t-direction. We thus require that 

V^nt = 0. (2.3) 

dt 

Now it turns out that the degenerate connection defined from the candidate degenerate 
connection coefficients shown above does not satisfy equation (2.3) in general. However, 
in particular cases it does, and this should be significant. 

Furthermore, it turns out that a degenerate connection which fulfills both (2.1) and 
(2.3) must necessarily be defined from the form the degenerate connection coefficients 
take in a GTCS. We choose such a degenerate connection on the grounds that it should 
reproduce the results of the above-mentioned candidate degenerate connection for those 
particular cases where it satisfies (2.3). It may be shown that this criterion yields the 
unique choice (vahd in a GTCS) 

f^t = (2-4) 

Note that the spatial indices occurring in equation (2.4) make the degenerate connection 
explicitly dependent on the structure of quasi-metric space-time. 

On AlxR, the Levi-Civita connections Vt define covariant derivatives along vectors 
tangent to M.. We extend the Vt by applying these covariant derivatives to the vector 
^ also. This must be done in a way that meshes with the degenerate connection V such 
that the torsion vanishes. By definition we then have (to avoid confusion we drop the 
f-label of Vt in expressions such as Vx) 

1 . d * d * d 

=r" = C = V a — =V a — = Va — , (2.5) 

t^t nt ^ I g^Q^ rndxt" ^ ^ 
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and this makes the torsion vanish by definition, since coordinate vector fields commute. 

Any vector field z on AlxR can be split up into pieces, where one piece z^-^ is 
tangential to A4 and the other piece is normal to At. To take the degenerate 

connection V along a given direction in Al xR it is sufficient to know V a and Va since 

J, dx'^ at 

the degenerate covariant derivative Vz is defined by using the appropriate Levi-Civita 

connection along the piece tangent to M. and using Va. along the remaining piece. In 

at 

particular, we find by direct calculation that 

The expressions in equation (2.6) vanish provided that g(t)oo ^-nd n^^^ do not depend on t, 
g(t)ojOct, n^'^^^(xt~^ and g{t)ij(xt^ (expressed in a GTCS). This dependence on t ensures that 
members of the metric family and the unit normal vector family are parallel-transported 
along ^. Thus the degenerate connection is compatible with the metric family. This 
guarantees that different metrics at different events can be compared unambiguously. 

Let A be some affine parameter describing the world line of an arbitrary test particle. 
The tangent vector ffeld ^ along the curve then has the coordinate representation ||^ + 
^Tsf"' 'w'^^re the are coordinate vector ffelds. Using the coordinate representation 
of we may deffne Vs =-^V9 + a ; this is the dcffnition of taking the covariant 

derivative along a curve in A^xR. The restriction to A/" is easily obtained by only 
considering curves in A/", i.e., curves where x° = ct in a GTCS. 

As an example we ffnd the degenerate covariant derivative of a family of vector ffelds 
Zf along an affinely parameterized curve x^{X) in A/". Expressed in (arbitrary) coordinates 
this is 

Whether or not the ffelds Zt are parallel-transported with respect to the degenerate 
connection V along the above curve depends on whether (2.7) vanishes or not. 

An important special case is when a family of vector ffelds is identiffed with the 
4- velocity tangent vector ffeld of the curve, that is 

where Tt is the proper time measured along the curve. Then we may calculate the 
degenerate covariant derivative of Ut along the tangent vector ffeld this yields the 
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degenerate acceleration vector field a defined by 

* * dt* dx"'^ dt* ,^ ^, 

a=V_8_Ut = — V8.Ut + — — V a Uf=—Va Ut + Bif (2.9) 

dTt ^* drt dTt ^* 

Using coordinate notation we see that (2.9) yields equations of motion 

d'^x^^ /*,, dt dx^\dx'^ rdTt\2*„ 

^ + (C^ + ri-„„,-)— (2.10) 

where ^ is constant along the curve since A and Tt are both affine parameters. We may 
compare (2.10) to the standard coordinate expression for the family of 4-acceleration 
fields at=VutUt compatible with the Levi-Civita connections Vt, namely 

However the equations (2.11) cannot be used as equations of motion in J\f since they do 
not include the dependence of on t. That is, the equations (2.11) only hold in Ai. 

A geodesic for the degenerate connection V (or equivalently for the metric family) in 
A/" is defined as a curve along which a vanishes. Clearly we need to specify a independently 
to determine which curves of Af are geodesies for V and which are not. We will show later 
(section 3.4) that in fact VjiUt = 0, so we have a=at from equation (2.9). Thus, similar 
to the Levi-Civita connection the geodesies for the degenerate connection are identified 
with world lines of inertial test particles. However, as we can see from equations (2.10) 
and (2.11), in quasi-metric theory the paths of inertial test particles are different from 
their counterparts in metric theory. 

Notice that since Ujr = cn^ is the 4-velocity of the FOs, the FOs always move orthog- 
onally to the FHSs. We may also show that we in general have 

* dx'^ 

«'^^^(tW = 0, (2.12) 

for any curve in A/", thus the degenerate acceleration field is normal to the tangent vector 
field (see equation (3.73) for direct verification). This is necessary to ensure that the 
length of Ut is constant along the curve. 

A simple illustrative example of a one-parameter family of 4-metrics is defined by the 
line elements (given in spherical GTCS (a;°, r, 9, (f) j) 

dsl = -{dx'f + (t^ + r'dn^) , (2.13) 

where dQ^=d9'^ + sm^Odcf)^, and to and Sq are constants. We see that this is a one- 
parameter family of static metrics, each one having the geometry of S^xR, but with 
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different scale factors. When doing kinematical calculations using (2.13) one must of 
course apply the equations of motion (2.10). Notice the difference between this situation 
and that one gets by applying the equations (2.11) to the single non-static metric one 
obtains by substituting cdt — dx^ into (2.13); the geodesies of this single metric are not 
equivalent to the geodesies calculated from the metric family (2.13). And although the 
null paths are the same, the role of i as a parameter along these paths is different for the 
two cases. (See section 4.2 for some calculations.) 

One immediate question is how to calculate the metric family gt, the obvious answer 
being that it should be found from appropriate field equations. However, this does not 
necessarily mean that should represent a solution of said field equations. Rather, we 
should take the more general view that it must be possible to construct the family g^ 
from another family gt which does represent a solution of said field equations, together 
with fields obtainable from the family gj. Why we should take this view is explained 
below. 

As mentioned earlier the NKE should be coupled to curvature of some sort. At this 
point, we may guess the qualitative form this "coupling" equation should take. To do that 
we remind of two points. Firstly, in contrast to the kinematical relations, this equation 
should not apply to an arbitrary spatial hypersurface, but only to the FHSs (where the 4- 
metrics of the family apply exactly). Secondly, the coupling equation should not depend 
on the surrounding 4-metrics via extrinsic curvature etc., but only via the evolution 
parameter. Thus this equation should essentially be scalar and it should look something 
like 

Ht=cFi%,Ft - c-^xt = /*. (2.14) 

Here Ft are the local scale factors of the tangent spaces of the FHSs in g^, dn denotes 
the Lie derivative in the direction normal to the FHSs and c^^Xf represents a possible 
kinematical contribution to the change with time of Ft. Furthermore ft are scalars 
constructed from the geometry intrinsic to the FHSs and the only constant of nature 
we have implicitly used this far, namely the velocity of light. 

We want the right hand side of (2.14) to be positive everywhere. This means that the 
NKE should take the form of a local "expansion" (defined by the FOs). In particular, 
doing cosmology within the quasi-metric framework, the Hubble law is interpreted as 
a special case of (2.14), and the cosmic expansion thus follows from first principles, in 
contrast to the situation within the metric framework. We will give the exact form of 
the non-kinematical coupling equation (NKCE) in section 3.3. 

From the above we expect that Af contains the non-kinematical "expansion" as part 
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of the geometrical structure. However, it may not be possible to eliminate all the effects 
of this "expansion" just by holding t constant. This means that in general each metric 
of the family gt should contain some non-kinematical features integrated into the metric 
structure. That is, for each metric of the family some non-kinematical effects are treated 
as kinematical effects. But this means that the family should in general not represent 
a solution of appropriate field equations, since such solutions should not contain any 
features coming from the NKE of the FHSs. Rather, as mentioned previously, the family 
gt should be constructed from a solution family gt and fields obtainable from gt. In 
section 2.3 we show that this construction involves a family of vector fields. 

2.1 Uniqueness and existence of the global time function 

In this section we investigate which restrictions should be put on each member of the fam- 
ily of metrics gt to ensure that they are acceptable according to the quasi-metric frame- 
work. The reason that we put the restrictions on gt and not directly on the "physical" 
family gt is that the former family by hypothesis represents a solution of field equations; 
if we know any general properties of solutions this may be useful when constructing such 
equations. 

Obviously we want the postulated global time function to exist always. That is, 
as observed earlier, there always should exist a unique invertible functional relationship 
between the coordinate time and t throughout the coordinate patch covered by an 
arbitrary local coordinate system {x^}. The requirement that the global time function 
exists leads to restrictions on the metric family gt. In particular we show that the 
existence of the global time function implies that it is possible to describe all the time- 
dependent aspects of the KE of gt with the help of a scalar function plus a spatial vector 
field. Besides, extrapolating gt off the domain FHSs, we require that each FHS is a 
member of a unique foliation of Ai (equipped with any allowable metric) into a family 
of spatial 3-hypersurfaces if we hold t constant and propagate the FHS out into the 
surrounding 4-geometry of A4. (We write A4 instead of A4 and J\f instead of J\f when 
we deal with the family gt.) 

To ensure the existence of a global time function, firstly it is important to ensure 
that the gradient field of the global time function in AT, i.e., dt = ^^rfx", is indeed 
timelike everywhere. (If it were not, a unit time-like vector field family fif would not 
exist everywhere on the FHSs.) Besides the world lines of different FOs should never 
intersect; this condition concerns the geometry of each 4-metric of gt in the sense that 
they cannot contain event horizons. By hypothesis the family gt should represent a 



55 



solution field equations. This means that if solutions of candidate field equations do 
indeed contain event horizons, such candidates must be discarded as inconsistent with 
the quasi-metric framework. 

Secondly, even if the FHSs are spatial Cauchy hypersurfaces, one must ensure that 
each FHS can be seen as a member of a family of spatial 3-hypersurfaces fohating A4 
equipped with any allowable 4-geometry gGgt. A sufficient condition to ensure this is 
to demand that the FHSs are compact (without boundaries) [3]. Thus, to ensure that 
the global time function exists, we will require that the FHSs are Cauchy and compact. 
Also this will be sufficient to ensure that the global time function is unique in the sense 
that it splits space-time into a unique family of FHSs. (The above result carries over to 
the asymptotically Minkowski case when we replace the condition of compactness with 
the boundary condition that the FHSs asymptotically can be identified with the fiat 
3-spaces of one of the global Lorentz frames of Minkowski space-time, namely the one 
where the global time coordinate =constant. We feel free to use asymptotically fiat 
FHSs in particular cases. However, doing this one may leave out crucial physical effects; 
it is important to be aware of that asymptotically flat FHSs may be used only as an 
approximation. Thus in principle the FHSs should be compact.) 

As observed earlier, the existence of a global time function singles out a preferred 
set of coordinate systems in A/", namely those where the identification = at is made. 
In the following we specify such a GTCS {x'^,x'^}, where the are spatial coordinates 
for i — 1, 2, 3. Using the notation of the ADM-formalism [4], an arbitrary metric family 
gt can be expressed in a GTCS via the family of fine elements (we use the signature 
convention ( — I- ++) throughout): 



N(t)sNl,) - N^ix", t)\ {dxy + 2N^t)i{x^, t)dxV + Qit^x", t)dx'dx\ (2.15) 



where Nt is the family of lapse functions for the FOs, the iV(t)j are the components of the 
family of shift covectors for the FOs and g{t)ij are the spatial components of the metric 
family (these are identical to the components of the metric family intrinsic to the FHSs) . 
Note that the form (2.15) is too general. That is, one should require the spatial part of 
(2.15) to be conformal to some "basic" 3-geometry to ensure that there will be no conflict 
between the NKE and the KE of the FHSs. Besides one needs to specify the role of the 
global time function t in (2.15) such that V^gt and VaUj vanish. 

dt dt 

As mentioned above, to ensure that the global time function exists we have restricted 
the 4-manifold AT in a topological sense by requiring the existence of compact Cauchy 
hypersurfaces (without boundaries). Besides there must exist a "basic" 3-geometry in 
the sense that the geometries g(^t)ij of the FHSs are restricted to be conformal to this 
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basic 3-geometry. A natural choice of such a basic 3-geometry is the geometry of the 
3-sphere since it is compact and since its positive intrinsic curvature fits well with an 
equation of type (2.14). We may take this as a definition, thus any metric family (2.15) 
should contain "prior" 3-geometry that cannot be changed by changing the dynamics of 
gravitating sources. 

Now any dependence on t in (2.15) should reflect non-kinematical aspects of the 
theory whereas any dependence on x° should reflect kinematical aspects. Since the non- 
kinematical evolution should take the form of an "expansion" , it should affect the spatial 
fields only (if local conservation of energy and momentum hold), so we guess that the 
family of lapse functions should in general not depend on t in AT. However, even if 
this holds, it is still possible that Nt may be a function of the quotient ^. Besides, in 
particular cases (in the early Universe and in extremely strong gravitational fields) one 
expects that local conservation laws may be violated. That is, when the global NKE is so 
strong that non-gravitational quantum fields cannot be treated as localized to sufficient 
accuracy, one should get spontaneous particle pair creation due to excitations of vacuum 
fiuctuations of such quantum fields. But besides these possibilities, we guess that the 
family of lapse functions is equivalent to the ordinary lapse function. We also guess that 
the explicit dependence on t of the spatial metrics g{t)ij enters (2.15) as a conformal 
factor. As justified in section 4, this factor is proportional to t^. But consistency then 
requires that the explicit dependence on t of the family of shift covectors must be linear 
and that the dependence of the family of shift vectors on t goes as |. Thus we may write 

H) = jNi,,)ix^)=jN^ix''), = l^x^, t)^9^t)ij%, (2.16) 

so that the family of shift vector fields is adequately represented by one ordinary shift 
vector field. The ^-dependence in equation (2.15) is now determined such that Ve.gt and 
Vd_nt vanish. But since may possibly depend on t (see above), the affine structure on 

dt 

AT will be slightly more general than that on Af. That is, the connection coefficients in 
77 read (compare to equation (2.5)) 

-0 ^t,t 1 Nt,t i - -a -a -a 

Tot — ~/Y^' ^jt ~ ("^ "/Vj"'^'^''' ^t=^tij.j ^i'ij.=^ {t)uix- (2-17) 

To avoid any restrictions on the dynamically relevant parts of gt, the dynamical content 
of the family gj should be adequately described by the lapse function family A^^. Thus, 
to ensure that the existence of a prior 3-geometry should not interfere with the dynamics 
of gt, it is sufficient that (2.15) takes the form 

t - - 
{dx^ + 2-Nidx'dx^ + -^N^Sij{x^')dx'dx\ (2.18) 

^0 ^0 



dSf 



Ns{x^,t)N'{x'')-N^{x^,t) 
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where Sijdx^dx^ is the metric of the 3-sphere (with constant radius equal to cto, see 
section 4). Note that the components of the inverse metric family g^"^ are given by [4] 

where Ji^^-^ = ^N^'^S^^ are the components of the inverse spatial metric family h^^. 

The consequences of (2.18) are rather far-reaching, as it means that the dynamical 
content of the family should adequately be described by the scalar field Nf for any 
theory of gravity formulated within the quasi-metric framework. This means that if one 
specifies initial conditions on an initial FHS it is sufficient to know the time evolution 
of one scalar field to predict the geometry of all of the hypcrsurfaces, i.e., to recover the 
space-time geometry. (This should be seen in contrast to e.g. General Relativity (GR), 
where analogously one must know the time evolution of the full extrinsic curvature tensor 
field since it is not restricted a priori.) 

We close this section by noticing that the proper time fjr of the FOs in JJ" is given by 
dfjT=Nt{x'^{t),t)dt, which may in principle be integrated along the world line of any FO. 
We thus see that (1) quasi-metric space-time is asymmetric; once time-orientation of the 
FHSs is chosen t must increase when fjr increases, (2) it is convenient to fix the scale of 
Nt by setting it equal to unity for an isotropic, globally empty toy universe (section 4), 
(3) it is possible to scale t (and thus 7>-) with a constant (this corresponds to an ordinary 
units transformation), (4) it is natural to choose i = for 7>- = as the beginning of 
time. Once this is done and the units and the scale of Nf are fixed, t is unique. 



2.2 The question of general covariance 

As mentioned above, the existence of a global time function implies that space-time 
can be foliated into a set of "preferred" spatial hypcrsurfaces, that is the FHSs. This 
geometric structure naturally defines a class of "preferred" coordinate systems, namely 
the GTCSs. The FOs by definition always move orthogonally to the FHSs, but the FOs 
usually have non-zero 3-velocities in any GTCS. That is, the FOs are in general not 
stationary with respect to any particular coordinate frame. This indicates that there 
is no preferred coordinate frame naturally associated with the quasi-metric geometrical 
framework. (Averaged over cosmic distances however, the FOs are at rest with respect 
to a particular set of GTCSs defining the cosmic rest frame. But this is not a preferred 
frame in the sense that the outcomes of local experiments depend on it.) 

It may be shown [15] that the local inertial frames of gt can be identified with local 
Lorentz frames. This means that in quasi-metric theory, there can be no local conse- 
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quences of curved space-time for non-gravitational physics and that local Lorentz in- 
variance (LLI) holds for non-gravitational test experiments. But since there should be 
no preferred-frame effects associated with quasi-metric geometry, there is no reason not 
to search for a quasi-metric theory of gravity where LLI holds for all local experiments 
including local gravitational experiments (such as Cavendish experiments). Thus there 
should be no a priori reasons to prefer metric over quasi-metric geometry as the basis for 
a new theory of gravity. 

But another question arises due to the natural existence of GTCSs within quasi-metric 
geometry. The GTCSs represent a set of preferred coordinate systems in quasi-metric 
relativity, so what happens to general covariance? It was recognized long ago that if one 
with "general covariance" understands "diffeomorphism-invariance" , it is physically vac- 
uous. In this sense the quasi-metric framework is, of course, generally covariant. Nor does 
"general covariance" in the sense of invariance of form under general coordinate transfor- 
mations contain any general physical principle. (We do not expect our framework to be 
generally covariant in this sense, since (2.18) is form invariant only for transformations 
from one GTCS to another.) Rather, the modern classification of theories depends on 
whether or not there exists "prior geometry" [1], [4]. This concept refers to any dynam- 
ically relevant aspect of the geometry that cannot be changed by changing the matter 
distribution. 

By construction the quasi-metric framework guarantees the existence of prior geom- 
etry. After all, we have restricted the 4-metrics of the family gt in such a way that they 
are of the specific type shown in (2.18). But as mentioned in the previous section, the 
existence of a prior 3-geometry does not necessarily interfere with the dynamics of the 
metric family g^. And even if this metric family has only one explicit dynamical degree 
of freedom, as we shall see later a second, implicit dynamical degree of freedom coming 
from the NKE is present in the "physical" metric family g^. So even if the structure 
of quasi-metric geometry implies the existence of prior geometry there is no reason to 
expect that this prior geometry lays any extra restrictions on the dynamics. 

In the next section we try to find for which general cases there exists some sort of 
correspondence with the metric framework. 

2.3 The metric approximation 

In the following we describe how to construct the family gt from the family gt and a 
family of vector fields. We do this by first studying a particular example which we will 
call a "metric approximation". A metric approximation is defined as a metric family 
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gt (or gt) which does not depend exphcitly on the global time function t; thus in effect 
we have only one metric. This metric necessarily represents an approximation since one 
neglects the effects of a changing global scale factor. Note that metric approximations 
make physically sense for isolated systems only. 

Metric approximations cannot work if the FHSs are taken to be compact since the 
NKE associated with the global curvature of space should result in global scale changes 
of the FHSs. And in quasi-metric theory global scale changes cannot be described by any 
single metric since the global curvature of space has nothing to do with gravity. (This 
is easily seen from the spatial prior geometry present in equation (2.18).) However, in 
particular cases it may be a good approximation to relax the requirement of compactness, 
replacing it with the requirement of asymptotic flatness. If the FHSs are considered as 
asymptotically fiat and embedded in asymptotically Minkowski space-time the problem 
of global scale changes is avoided, since wc expect the NKE to vanish for fiat FHSs. But 
we still want to identify the FHSs with the hypersurfaces x° =constant, where is the 
global time coordinate we started out with. That is, when making metric approximations 
we keep the unique global time coordinate wc had to begin with and define the FHSs 
from the condition x^ =constant. Moreover, wc may define a GTCS as a global time 
coordinate system where this global time coordinate is used. But note that after having 
made a metric approximation, it is possible to choose other global time coordinates than 
the one we started out with. 

Metric approximations can be constructed from any one-parameter family gt of the 
form (2.18) by setting ^ to unity and replacing Sijdx^dx^ by an Euclidean 3-metric. 
Thus any metric approximation may be written formally as gt^g (or gi^g); that is, 
metric families are approximated by single metrics. Now we will study the particular case 
where a GTCS exists in which the family gt does not depend on x° and the shift field 
Nj vanishes; this is called the metrically static case. Besides, for reasons explained later, 
we will require spherical symmetry. In this case the metric approximation of the family 
gt is a single static, spherically symmetric metric g. That the metric g is independent 
of time means that the FHSs are static, i.e., there is no net change of scale anywhere 
on the FHSs. But this does not mean that there are no effects of the NKE. Rather, to 
construct a new metric g from the metric g we need to put the NKE of the FHSs into 
the metric structure of g. This means that if the FOs are considered as stationary for the 
metric g, then the stationariness of the FOs for the metric g should be seen as a result 
of both the non-kinematical expansion and the motion of the FOs taken together. That 
is, the FOs in AT may be perceived as moving with a velocity — v say, with respect to the 
FOs in A^. Of course this motion is cancelled out when one adds another vector field v 
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constructed from the NKE. Thus one may think of the FOs as "expanding" and "faUing" 
simultaneously such that the combination results in stationariness. However, this way of 
recovering stationariness is not equivalent to being stationary with no NKE in the metric 
g. We may find the difference between g and g by quantifying the two different versions 
of stationariness in terms of the vector field v. As we shall see the effect of v is that 
one must introduce correction factors to the components of g; these corrections yield the 
wanted metric g. 

Before proceeding we should choose an appropriate coordinate system. When there is 
spherical symmetry with respect to one distinguished point we naturally define a spherical 
GTCS (a;°,r, ^,0) with the distinguished point at the spatial origin. The distinguished 
point plays the role as a natural "center of gravity". Now, according to (2.14) the NKE 
takes the form of a (in general different) uniform increase in scale for each tangent space 
of the FHSs. Due to the spherical symmetry this increase in scale can be expressed 
as a purely radial, outwards-pointing 3- vector field v = v'^{r)^. We thus define v by 
applying the "Hubble law" in the tangent space of each point of the FHSs. Besides, the 
distance £ to the spatial origin in the tangent space of each point on the FHSs is given by 
£ = r^/g^. The radial 3-vcctor field v and its speed v=\/v^ are thus defined by (where 
H is the metric approximation of Ht given in (2.14)) 

v'^(r) = H(r)r, v(r) = H(r)e{r) = H{r)^grr{r)r. (2.20) 

We now calculate the new metric g when g is static and spherically symmetric, as assumed 
above. In this case the spherical GTCS defined above is in fact a HOGS and the metric 
takes the standard form 

ds" = -B{r) {dxy + A{r)dr^ + r^dQ^, (2.21) 

where A{r) and B{r) are smooth functions without zeros and dVL^=d9'^ + siv?9d(j)'^. 

We now try to find a new metric g from the metric g of the type (2.21) and the vector 
field V. Prom the NKE the FOs in AT get a net 3- velocity y — v'^^in the positive radial 
direction of the spherical radial coordinate system according to equation (2.20). To get 
rid of this the FOs in JJ must move with a velocity —v. That is, the static condition of 
the FOs in jV" may be seen as the result of adding two opposite 3- velocities. Moreover 
all observers in H are also static, but in a way that does not involve any "motion" by 
definition. We now want to find the time and length intervals as expressed by the static 
FOs in M. We do this for each tangent space, and taken together the result will be a 
new metric g derived from the metric g given in equation (2.21). And just as g is the 
metric approximation of a spherically symmetric, metrically static family gt, the new 
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metric g may be thought of as the metric approximation of a new spherically symmetric, 
metrically static family g^. 

We proceed in two steps, each step giving the same correction factors to the metric 
(2.21). Step one takes into account how the length and time intervals of "static" FOs are 
seen by observers moving with 3- velocity — v and by definition "at rest" , before taking 
the NKE into account. That is, the "static" FOs are not "at rest" but move outwards 
with velocity v compared to the observers "at rest" when the NKE is not taken into 
account. In each tangent space this extra motion yields a correction to the radial interval 
due to the radial Doppler effect, the correction factor being [^izfj ■ There will also be 
an inverse time dilation correction factor (1 — ^)^/^ to the time interval. There are no 
correction factors for the angular intervals. 

Step two transforms the effects of the NKE in AT into correction factors to the metric 
coefficients. In 77 the NKE is described as a velocity v in the positive radial direction 
for each tangent space, and the effects of v as seen by static FOs, now by definition "at 
rest" , are taken care of by using a pair of correction factors identical to those in step one. 

So the new metric g constructed from the metric g given in equation (2.21) and v 
reads 

ds' = -(l - "^fBirXdxr + (i±^)'A(r)dr^ + rW, (2.22) 

c 

and it is clear that we must have v{r) < c everywhere if the above expression is to make 
sense. Note that any metric approximation works because equations (2.10) reduce to the 
geodesic equations (2.11) for metric approximations. 

It must be stressed that any metric approximation is adequate only if one can sensibly 
neglect the NKE coming from a global scale factor. This observation leads to the question 
which vector fields one should use when constructing g^ from gt for the general case. The 
answer clearly depends on the evolution of the FHSs. That is, the FHSs are not stationary 
as is in the above case but will in general expand. However the FHSs may in principle 
even contract in some places due to gravitation. This means that the net change of 
scale for each tangent space at points of each FHS may be positive or negative. But to 
construct the radial vector fields for the general case we are interested in that part of 
the evolution of the FHSs which is not "reahzed", i.e., not refiected directly in g^. When 
a metric approximation is sufficient this non-realized part of the evolution is represented 
by H. In the general case, however, the non-realized part of the evolution of the FHSs is 
not represented by H but some other quantity yt, say. 

This means that for the general case, we should replace Ht with yt. But besides 
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this there is a more subtle problem. Namely, that there is no obvious choice of "the 
center of gravity" in the absence of spherical symmetry (this is why we imposed spherical 
symmetry when deriving equations (2.20) and (2.22)). That is, in the general case one 
needs to find a 3- vector field family bjr representing "a local distance vector field" from 
some fictitious center of gravity, and such that bjF reduces to the radial coordinate vector 
field for the spherically symmetric case. We will define hjr in section 3.3. Meanwhile 
we show how to construct gt from and when it is assumed that b^F is known. 

To begin with, using the quantities yt and bjr, we may define the vector field family 
just as we defined equation (2.20); this 3- vector field family represents the generalization 
of the vector field v used for spherically symmetric metric approximations. Instead of 
(2.20) we define (using an arbitrary GTCS) 

'^{t) = yt^-F: V = yt\l9(t)ik^T^r- (2-23) 

Note that the explicit ^-dependence should cancel out in the second formula and that 
we can eliminate any possible implicit ^-dependence via Nt by making the substitution 
t — x°/ c. Also note that the NKE in the tangent space of each point of the FHSs still 
arises from Hf. But the part of the NKE coming from a global scale factor must be 
subtracted from Ht when constructing v^; we thus get (2.23) instead of (2.20). We will 
call the effects of the NKE associated with the local effects of the NKE. 

In the following we will write down the formulae for the components of gt which must 
hold for the general case. The main difference from the spherically symmetric case is 
that the "extra motion" of the FOs in AT is locally along the bjc-direction rather than in 
the radial direction. But besides that, the principles behind the construction of gt are 
similar. To begin with we notice that the transformation formula should be identical to 
the spherically symmetric case when no spatial components are involved. That is, as in 
equation (2.22) we define 

g{t)oo = (l - ■^)^^(t)oo- (2.24) 

To guess reasonable transformation formulae when spatial components are involved it 
is convenient to define the unit vector field eb=^el^ and the corresponding covector 
field e^=j^e\dx'^ along bjF- But it is in general not possible or practical to construct a 
GTCS where e;, is parallel to one of the coordinate vector fields. For this reason one 
expects that the general transformation formulae for spatial components should involve 
components of e*. By trial and error (involving the spherically symmetric case using 
non-spherical coordinates) we find the following transformation formula for the purely 
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spatial components of g*: 



9{t)ij = 9{t)ij + 




(2.25) 



The transformation formula (2.25) has the correct correspondence with the result found 
in equation (2.22) and it may be shown that this correspondence holds in any GTCS. 
Moreover, from equations (2.24) and (2.25) and consistency requirements we may find 
the remaining components of gt, namely 



Given equations (2.24), (2.25) and (2.26) it is now straightforward to find the transfor- 
mation formulae valid for g^^. 



Our goal now is to construct a theory of gravity consistent with the quasi-metric frame- 
work laid out in section 2. The class of such theories must satisfy the postulates: (1) 
space-time is constructed as a 4-dimensional submanifold A/" of a 5-dimensional product 
manifold A4xR as described in section 2. This implies that space-time is equipped with 
a global time function, an one-parameter family of Lorentzian 4-metrics gt (constructible 
from another family gt where each metric is restricted by the constraints stated in section 
2.1), and a degenerate connection V compatible with the metric family, (2) inertial test 
particles follow curves calculated from equation (2.10) with the condition = and 
where the relation between a and a^ must be specified, (3) in the local inertial frames of 
the metric family g^, the non-gravitational physics is as in special relativity (SR). 

Postulate (3) means that it should be possible to identify the local inertial frames 
of gt with local Lorentz frames, see [15] for verification. This implies that the quasi- 
metric framework by itself induces no local consequences of curved space-time on the 
non-gravitational physics. Thus both LLI and local positional invariance (LPI) should 
hold within the quasi-metric framework. But for any finite scale the free fall curves of 
gt as obtained from from (2.10) do not coincide with the geodesies any single space- 
time metric, see [15] for a further discussion. Furthermore this deviation depends on 
position and on the velocity with respect to the FOs. Such dependence implies that 
non-gravitational laws may appear different from their counterparts in any metric theory 
for any finite-sized observer. 




(2.26) 



3 Equations of evolution 
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On the other hand, postulate (2) ensures that the weak equivalence principle (WEP) 
is fulfilled (i.e., that test particles' trajectories do not depend on the particles' internal 
composition). The sum of the WEP, LLI and LPI is called the Einstein equivalence 
principle (EEP); every metric theory of gravity satisfies the EEP [1]. The proposition 
that every consistent theory which satisfies the WEP also satisfies the EEP is called 
Schiff's conjecture (SC) [1]. Thus SC continues to hold for any theory consistent with 
the above postulates. 

From equation (2.18) we see that how the scale of the FHSs varies from event to event 
is given by the scale factor Ft=Ntct. One may interpret this scale factor as resulting from 
a position-dependent transformation of length units. Because any dependence on scale of 
physical laws should be experimentally determined only by the variation with position of 
dimensionless numbers, we are free to define as constants physical constants that cannot 
be related to other such constants by dimensionless numbers (sec e.g. [5]). Thus, by 
definition the velocity of light c in vacuum and Planck's constant h are constants, i.e., 
independent of Ff. But when writing down line elements describing the geometry of 
space-time, some set of units are implicitly assumed used. Usually one docs not separate 
between units defined from purely non-gravitational systems and other types of units. 
However, this must be done if one wants to describe position-dependent scale changes 
between non-gravitational and gravitational systems. 

We are thus led to the notion that any fixed "atomic" units, i.e., units operationally 
defined by means of non-gravitational forces only, may be considered as formally variable 
throughout space-time; this applies both to M. and to JJ". That c and h are not formally 
variable by definition implies that the variation of length and time units are similar 
and inverse to that of energy (or mass) units. One may expect that the formally variable 
units should represent scale changes between non-gravitational and gravitational systems, 
where said scale changes should be found from a theory of gravity. 

This means that all physical laws are formulated in terms of "variable" atomic units. 
That is, gravitational quantities are considered as variable measured in such units. But 
which quantities do we define as gravitational? Not unreasonably we define the family 
gt as gravitational since gt is by definition a solution of gravitational field equations. 
Moreover the Newtonian gravitational "constant" Gt and the stress-energy tensor 
considered as a source of gravity should be defined as gravitational quantities. 

It is convenient to introduce a scalar field which defines how atomic time units 
vary in space-time. Any gravitational quantity should then be associated with \E'( to 
some power; this factor should tell how the field changes due to the formally variable 
units. To illustrate this point we note that Gt has a formal variation similar to that 
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of time units squared. Since Gt is a gravitational quantity this means that the extra 
function associated with Gt is a factor ^7^- That is, measured in atomic units Gt will 
not be a constant. Now, one obvious consequence of this is to postulate that Gt varies 
as But since Gt always occurs in combination with masses it is more convenient 

to define Gt to take a constant value G in atomic units and rather differentiate between 
active mass and passive mass. That is, active mass plays the role as a source of gravity 
and varies measured in atomic units. On the other hand passive mass (inertial mass or 
passive gravitational mass) does not vary and should be used in atomic interactions. In 
the following we will denote any active mass rrit (corresponding to a passive mass m) 
with an index. 

Mass considered as a gravitational quantity scales via the variation of mass units, 
which vary as Thus seen active mass is associated with a factor \E'f. But since 

we have defined the gravitational coupling parameter to be constant wc must explicitly 
include the variation of the gravitational coupling into the active mass. This means that 
active mass rrit varies as "^t^ measured in atomic units. Furthermore, since matter fields 
should be coupled to geometry via the active stress-energy tensor Tt which dimension 
is (active) mass density, considered as a source of gravity is associated with an extra 
factor This means that experiences extra changes due to the formally variable 
units, and this must be taken care of in local conservation laws. 

The question now is how (or equivalcntly rrit) changes from event to event in AT. 
Since line elements of the type (2.18) are measured in atomic units we propose that any 
change in the scale factor Ntct can be explained solely by means of the formally variable 
units. This will work if 

- *. = ^*.. ".. = ^"0. (3.1) 

where \E'o, miQ etc. refer to values taken at some arbitrary reference event. By convention 
we may choose iVo = 1, \l'o = ^ which means that the chosen (hypothetical) reference 
situation is an empty universe at epoch to (see section 4). Note that the formal variability 
of gravitational quantities should be unaffected by the construction of the family from 
the family since this construction only involves the formally invariant vector fields Vj. 

The strong equivalence principle (SEP) is essentially a generalization of the EEP. 
That is, for a theory to satisfy the SEP, the EEP must be vahd for self-gravitating 
bodies and gravitational test experiments as well as for non-gravitating bodies and non- 
gravitational test experiments, see e.g. [1]. Since theories consistent with postulate (3) 
satisfy the EEP, one might hope to find a quasi-metric theory that satisfies the SEP as 
well. However, since active mass is not equal to passive mass the SEP in its most stringent 
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interpretation cannot be satisfied within the quasi-metric framework; all we can hope for 
is that the gravitational weak principle of equivalence (GWEP) (that the WEP is vahd 
for self-gravitating bodies as well) holds. 

But the result of any laboratory experiment testing the constancy of G should reveal 
only the variation of G via the formally variable units and not any extra, independent 
kind of variability. Thus we should be able to investigate the dynamical properties of 
without having to worry about a "physically" variable G, independent of the formally 
variable units. Besides the formal variability of G should depend only on quantities 
constructed from pieces of the metric family gt and not on any extra gravitational fields. 
But this means that we may be able to construct a theory of gravity, consistent with the 
quasi-metric framework, which satisfies the GWEP. In fact the GWEP holds since it may 
be shown that active mass dust particles really move on geodesies of gt and of g^ [15]. 

3.1 Projections 

If we define a GTCS {x'^} = this can be used to find the coordinate repre- 

sentations of the orthonormal vector field basis {iit, Cj} defining the moving frame field 
spanned by the FOs, where ej=^. The vector field family fit is identified with the 
unit-length vector field family = —1) tangent to the world fines of the FOs, 

but is in general not parallel to the coordinate vector field For the coordinate 
representation of fit in a GTCS, see below. We can use iit and its corresponding covector 
field together with the metric tensor family intrinsic to the FHSs to project space-time 
objects into space objects intrinsic to the FHSs. That is, repeated scalar products with 
eventually reduce space-time objects to space scalars, whereas repeated scalar products 
with ht eventually reduce space-time objects to space objects of the same type and rank. 
Combinations of these operations are possible. To avoid confusion about whether we 
are dealing with a space-time object or its spatial projection, we will label space objects 
with a "hat" if necessary. That is, the space objects are the spatial projections of the 
space-time objects A^. 

Component-wise, we define a label ± to indicate a projection by —fit, whereas spatial 
covariant indices indicate spatial projections. Note that contravariant spatial indices do 
not indicate spatial projections; here one must use the "hat" label to differentiate between 
the components of a space-time object and its spatial projection. As an example we take 
the space vectors Q''(^i^j = — h^(t)(3Q^{t)a^''(t) which arc mixed projections of the the space-time 
tensors Qft)^- The corresponding space covectors are Q^^^jj_= — h^^^,^-Q{t)i3a'^'{ty 
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The decomposition of and its inverse with respect to each domain FHS is given by 

where h(t)^u are the components of h^. We see from (3.2) that h(^t)fiu represents the re- 
striction to the domain FHS of each space-time metric g(t)fiu- In the GTCS, hj considered 
as a space-time tensor field, has the potential non-zero mixed or time components 

kt)oo = A^fe7V^ h^tyo = km = l^Nj, hi^^, = ^jN^ = ^jN-'S^'N,, (3.3) 

whereas has the coordinate representation 

Using (3.2), (3.3) and (3.4) we may verify that h[i)n(t)v = h(t)i^^^lt) = 0' 9{t)ll = -1> 
= 9(t)jl = 0, = kt)ij and ^f^^ = /ij^^ - 
Any metric gGgt and gGg* (where gt=iVt~^g() define symmetric linear metric con- 
nection fields V and V respectively. In component notation, we denote the operation of 
V on space-time objects by the symbol To illustrate this by an example, if are 
vector fields, the components of the tensor fields VY^ are y(t)-p- The link between the 
Christoffel symbols corresponding to V and V is 

i (t)/37 - ^ ~ ~ ~N^^ -^9(t)9{t)pj- (3.5) 

By projecting the connections into the FHSs, we can define hypersurface-intrinsic con- 
nection fields At and At compatible with the hypersurface-intrinsic metrics ht and ht 
respectively. In component notation we denote the operation of At on space objects by 
the symbol '|'. (Since we never mix fields constructed from the different families gt and gt 

we let the symbols ';' and '|' denote space-time covariant derivation respectively spatial 
covariant derivation also in A4. Doing otherwise would imply a serious risk of running 
out of appropriate symbols.) 

Also, the manifold structure of Ai (and of A4) naturally defines a Lie derivative £x 
which operates on space-time objects. We may also define a spatial Lie derivative £x 
which operates only on space objects, where the vector field X is still treated as a space- 
time vector field. See ref. [7] for the detailed projection formulae relating £x and £x- 
Analogously the manifold structure of J\f (and of J\f) defines the Lie derivative £x and 
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its spatial counterpart >Cx- The only differences between the Lie derivatives in A4 and 
Af arise from the i-dependency of the metric family. 

Furthermore the metric families gt and gt a-re compatible with the degenerate con- 
nections V and V, respectively. It may be easily shown that the link between the corre- 
sponding connection coefficients is 

Vpt = (>p. r^T = r^T - - -^^p + -^9{t)m'yp- (3-6) 

In component notation wc use the symbols and to denote the operations on space- 
time fields of V and V? respectively. Just as for metric space-time geometry one may 
construct higher rank tensor fields by using these operators. As an example of how this 
is is done, take a vector field family Yt and construct VY^ defined by its components in a 

GTCS; ^(")ii=^(");i, Y(t)*o=^(ty,o + ^ ^^(t)*t = ^(t);o + ^ ^^(t)^t +c ^f"^^(t)- Since the extra 
terms represent the components of a 4-vector, VYt is a family of rank 2 tensor fields. 
Similarly the covariant divergence V-Y^, or in component notation ^(t)ja=c^^(t)*t + 
Yi^y,a = {cNt)-'iY^t)l + is a scalar field. 

Next we define the 4-acceleration field ajF corresponding to each normal unit vector 
field n(t) by c~^a^=n|^j-).^n^j-). This is a purely spatial vector field since a^x = —ajr^n^_^^ = 
—aj^n^t^i^ — 0. Besides, from (3.4) and the definition of a.jr we find 

= ^InTV,. (3.7) 

Furthermore we define a family of extrinsic curvature tensors describing how each 
FHS is embedded in the relevant g&gt- That is, we define 

K(t)fiu= - h1t)ijhlt)v^{t)a;p, (3-8) 

which is symmetric since there is no space-time torsion in Ai. The extrinsic curvature 
tensors of the FHSs in are defined analogously. Projection shows that each K^, 
is a purely hypersurface-intrinsic object, i.e., -^(t)Xl = ^{t)lj — ^(,t)Il — ^(t)Ij = 0- 1^ 
can also be readily shown [2], [4], [7] that K^, may be written as K^i-^^^, = —^£nth(t)nu, 
^{t)tiv — —\£^Tith{t)nvi respectively. After some calculations we find the spatial compo- 
nents 

= -2^nthit)ij = + ^jK) + ~ J^~^^:Fk^)h(t)ij, 

i^i - toN\ .Nt,j_ to -2- 
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(3.10) 



Note the difference between tfie spatial components of £^nth(t)tiv and Cnji{t)ij'i the latter 
Lie derivative acts on spatial objects only and this means that should be treated as 
a tensor on space when using Cn^ but as a tensor on space-time when using £ii^. The 
relationship 



(3.11) 



follows directly from the definitions of the Lie derivatives. Relationships similar to (3.11) 
can be found between the different Lie derivatives of any tensor field. 

Prom (2.18) we can find the first derivative of with respect to t. This is 



-Qlht)i3 = cNt{£iith^t)ij + '2K(^t)ij) = 2(- + ~)h(^t)ij, 



(3.12) 



^1%^ = cN,{hji%) - 2K^ = -2(1 + ^)^- . 



We now write down how ajrj varies in the direction orthogonal to the FHSs. We find 



(3.13) 



h - -r - ^Of- ^'i -2- 
and similar formulae for other covectors, and moreover 



t 



d 



(3.14) 



£r.taTj = \c 'a:fj,o -"-^[aTkN%j ), c = - c ^a:Fj^, (3.15) 



and still furthermore we find 



We also get 
where 



2 ,1 iV,„, . toN\, 



cNtH Nt 



(3.16) 



(3.17) 



(3.18) 



(3.19) 
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The projection of metric space-time covariant derivatives into spatial objects expressed 
in components may be subject to some confusion of notation. The important thing 
to remember is that when a component contains a metric space-time derivative, the 
projectors ht and hj are always used outside the derivative, whereas if the component 
contains a space derivative, the projectors are used inside the derivative. For example 
a mixed projection of the rank 2 space-time tensor fields VY^ yields the components 
y(j)j.^= — Y(t)i^.jn'^^^ and not — whereas the corresponding components of the 

spatial covector fields Ay^ are y{t)l\j= - 

Formulae expressed in component notation for various projections of metric space- 
time covariant derivatives of general families of space-time tensors may be found by 
finding the projection formulae for VZ^, where is a some family of space-time covector 
fields. One may then apply these formulae to every covariant index of a general family of 
totally covariant space-time tensors. Such formulae are derived in e.g. [2] (see also [7]), 
and are given by 

Z{t)i;j = Z(t)i\j - Z(^t^lK^t)ij, 

Z{t)j;l = -'^n(^(t)j + c"^ajpj(y-^^(t)i - ^(t)l) - (3.20) 

Similarly we may find the projection formulae for the the tensor fields VZj constructed 
by using the degenerate connection in 77. The result is 

Z{t)l*j = Z(t)l;j^ Z{t)iij = Z(t)i;j, = —^tZ(t)l — C '^a'jrZ(^t)i, (3.21) 

Z{t)j*l = -^tZ{t)j + c ^ajrj(-^— ^(t)i - 'Z'(t)i) 

+{cNtr\] + §f - ^mHty (3.22) 

Since we will explicitly need the projection formulae for the metric covariant derivative 
of tensor fields of rank 2, we fist the relevant ones. If Qt is a family of such fields, then 
we have [2] 

Q{t)li;j = Q{t)li\j — K{t)ijQ{t)l±_ — K(t)jkQ'lt)i: (3.23) 
Q{t)ij;k — Q{t)ij\k — -f^(t)ifeQ(t)lj — K{t)jkQ {t)li, (3.24) 

Q(t)ll;l — —^ntQ{t)ll — c ^a:Fj(5(t)i — c ^ajFi<5(t)l ) (3.25) 
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to iV* 

Q{t)jL;l = -^ntQ(t)jl + (C ^rj — ^ - Klt)j)Q( 



t N, 



(t}ij 



-c ^aj^iQu.j - c ^arjQ(t)ll, (3-26) 



-2- f-O^V -J. 2- 



Q{t)ij;l — -^ntQit)ij + {C I^(t)j)Qit)ik - C ajrj(5(j)jjL 



t 

+(c~^a^i j-^ - K^t^i)Q(^t)kj - c~'^ariQ{t)lj- (3-27) 



The analogous formulae for the degenerate covariant derivative of a family of tensor fields 
of rank 2 are 

Q{t)li*o = Q{t)li;jJ Qit)ij*k = Q{t)ij;k, (3.28) 

Q{t)II*I — —^tQ{t)II ~ c ^a;Fi(5(t)i — c ^a:Fi(3(t)i , (3.29) 



+ (ciV,)-i(^ + ^)Qit)jl - c-^ari%)^ - c-^a^,Qit)ll, (3-30) 
<5(t)ij*l = -^tQwii + - ^5)J<5(t)fcj + ^(- + 



^ - K^^)j)Q^t)ik - c~^{aj^jQ{t)il + a^iQ(t)lj)- (3-31) 
Note that we may show that h(^t)ij*l — from (3.31) by setting Q[t)ij — h(t)ij- 

3.2 The field equations 

In this section we use the formalism described in the previous section to find local conser- 
vation laws for the active stress-energy tensor and to express these laws in hypersurface 
language. But as we shall see the local conservation laws in quasi-metric theory are dif- 
ferent from their metric counterparts. We also postulate field equations which must hold 
for the metric family gf Notice that in quasi- metric theory, the stress-energy tensor is 
really a family of rank 2 tensor fields, this is what its t-label is meant to indicate. 

By construction the quasi-metric framework satisfies the EEP. But this does not 
imply that one should automatically use the usual rule "comma goes to semicolon" when 
generalizing the non-gravitational physics from flat to curved space-time. The reason 
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for this is that it is possible to couple non-gravitational fields to first derivatives of the 
spatial scale factor Ft such that the EEP still holds, see below. In particular we find 
that the covariant divergence V-TtT^O (expressed in terms of formally variable units) in 
the quasi-metric framework. It must be stressed that this does not necessarily imply any 
violation of local conservation laws in Af; rather it means that the conservation laws take 
unfamiliar forms in quasi-metric space-time. Any direct coupling to curvature does of 
course not occur in the conservation laws; rather there is coupling to other fields, namely 

- Nui Nt,t „„J _1_ 
U'Ti, , d-iii^ jvtrf- 

Within the metric framework it can be shown [1] that validity of the local conservation 
laws in the form V-T = follows from universal coupling (that gravity couples similarly 
to all non-gravitational fields, this follows directly from the EEP) and invariance of non- 
gravitational actions. Thus in any metric theory based on an invariant action principle 
the local conservation laws in their usual form V-T = are valid independently of the 
detailed form of the gravitational field equations. (Besides it can be shown [1] that the 
local conservation laws applied to the stress-energy tensor for a perfect fluid directly yield 
the Eulcrian equations of hydrodynamics in the limit of weak gravitational flelds and slow 
motions, establishing a Newtonian limit.) On the other hand, since c~^ajrj = —Ntt^j, 
^ = f^ip + ^IJc-^a^. = -Ntl\^ and i + ^ = |f in a GTCS, any coupfing to these 
flelds may be made to vanish in the local inertial frames so the EEP still holds. For this 
reason, in a GTCS the local conservation laws may be expected to take the unfamiliar 
form V-Tj7^0 for any theory compatible with the quasi-metric framework. (Note that 
one may straightforwardly show that the Newtonian limit of the local conservation laws 
applied to a perfect fluid yields the Eulerian equations of hydrodynamics even for the 
case V-TtT^O, see appendix B for some details.) 

By exploiting the form the local conservation laws take given the coupling dependence 
on OjFj, ^ffi and we flnd constraints which must be satisfled for any candidate 
geometrical fleld to which we want to couple T^. To begin with, using (3.6) we flnd that 
(expressed in atomic units) 

^ {t)lJ.*v — ^ {t)tji*v ^ Wm ^ (t)a 

= '^{t)ix*u + '^{c~'^o.jriTlt)i^ ^^(t)Ml) + - ^(t)i)> (3.32) 

where (3.32) compares the covariant divergences V-Tj and V"Tt, neither of which is 
expected to vanish. 

To get any further we need to know explicitly how the units vary; i.e. how ^'j (and 
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nit) in general depends on aj^i, ^j^, and |. This is easily found from (3.1). We get 

"^t,t = ^7"^lv~^"^*' '^t,!^ ^^m, mt.j^c aj^jUit. (3.33) 

It is also convenient to have expressions telling how and rrit vary in the direction 
normal to the FHSs. Such expressions are 

"'t, 



Note that we can use (3.34) to show that we may write 

mt*a = {mt,l + "^0^W°^ ^(*)°' mtio=c''^mt,t +mt,o , (3.35) 

and similarly for other scalar fields. Since Tt considered as a source of gravity has 
the dimension of (active) mass density, Tt has an extra variability like \E'^ measured in 
atomic units. For convenience we put all effects of this into the term T^t)iA.*u (3.32). We 
now want to project the relation (3.32) into the FHSs using the results of the previous 
section. We have 2 different projections; one yielding a spatial scalar field, the other a 
spatial vector field. These are T^t)l*v '^{t)j*i^^ respectively. Using the formulae (3.28), 
(3.29) and (3.30) we find 

~ -^*^(t)ll ~ ^{t)ikTlt), (3.36) 

'^{t)j*v — ~T(t)jiii + hft)T{t)ji*k = Ait^(t)jl + c ^ajc-jT(*,)^- + c ^ajFjT(t)ii 

^nm. + - '-'-'r/-fS;T,,,.. (3.37) 

Furthermore, using (3.32), (3.36) and (3.37) we get, after some short calculations 

Ait^(t)ll = Tl^t)l*u + -^t - ^{t)ikTlt) 

Nt,i — 

— + 2c"^a^iT(',)i - (3.38) 
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T, 



{t)j*u 



+ 



+ 3c ajTiT'l^^. 



mi 



(3.39) 



.to 

At 

To get any further with (3.38) and (3.39) , we must have an idea of how to find independent 
expressions for the terms T^^^-^j^-^ and Tl^t)j*u- do that we set up the general expressions 
(exphcitly including the effect of the formally variable units) 



T, 



{t)j*v 



+2 



— (- + 



(t)±± 



T, 



(3.40) 



■ {t)±*u 



(t)J 



Nt,: 



2 A^Nt,t. 



T, 



where ^i, ^2, ^^^d ^4 ^^e constants to be determined. As illustrated in appendix B, we 
must have ^1 = — 1 to ensure that the Newtonian limit of the resulting conservation laws 
comes out correctly. Furthermore we guess that ^2 = 1, since this choice seems natural in 
the case of a null fluid. For consistency it is then reasonable to choose ^3 = 1 and ^4 = 1. 
Inserting the above results into (3.32), (3.38), (3.39), (3.40) and (3.41) we get 



2 1 Nt,t.^ 



(3.42) 



{t)j*U ^ (t)r,V 



lu,..,. —(- + -^)-L(t)±j^ 



cNtH Nt 



(3.43) 



-2 Nt,_^ 



(t)ll, 



Tair. = 2c-^a^,f^,, - 2^T.),i, (3.44) 



Nt 



^fitT(t)l± — {^t - 2-^)T(t)ii + K^t)ikTlt) - T( 



N 



(3.45) 



-2 toN' 



I -2- rpi I -2- ''U rp rpi 

+c ani (t)j + c 7 ^9-^ {t)il - ^ {tW 



(3.46) 



Notice that we have separated out the dependence on t in (3.42) and (3.43) since this 
explicit dependence has no kinematical signiflcance. Similarly, when we take into account 
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the explicit dependence on t in the Lie derivatives, we find from equations (3.38), (3.39), 
(3.40), (3.41), (3.45) and (3.46) 



^tT{t)ll - ^ntT^t)ll - + -^)T(t)ll, (3.47) 

^tT(t)jl = ^ritT{t)jl - ^(^ + -^)T(t)jl- (3.48) 

Note that (3.47) and (3.48) tell how the units vary with t; in (3.48) the effects of this on 
the components of a spatial 3-covector is also included (which is why (3.47) and (3.48) 
depend differently on t). 

The relations (3.42), (3.43) and (3.44) (or equivalently (3.45) and (3.46) together with 
(3.47) and (3.48)) represent the local conservation laws for active mass-energy in 77. Note 
that these laws are expressed in atomic units and that the laws take care of the fact that 
gravitational quantities get extra variability measured in such units. Also note that 
equations (3.47) and (3.48) may be violated in the case of local non-conservation (which 
is expected to occur in the early universe and in extremely deep gravitational wells) . On 
the other hand equations (3.44), (3.45) and (3.46) are required to hold always. 

Next, the condition that we must have metric correspondence with Newtonian theory 
for weak stationary fields leads to a postulate, namely {k=^^) 

2^(t)ll = 2(c-^a>|, + c-%^,a> - + U,Kt) = k{T^,)Ii + fl,^,), (3.49) 

where is the Ricci tensor field family corresponding to the metric family (2.18). Note 
that whereas equation (3.49) is postulated rather than derived, it is by no means arbitrary; 
in fact equation (3.49) follows naturally from a geometrical correspondence with Newton- 
Cartan theory. Moreover, it is possible to write equation (3.49) in the form 

N,^% = Nt{R^t)ll + K^t)ikK% - U,K,). (3.50) 

According to the maximum principle applied to closed Riemannian 3-manifolds (see ref- 
erence [16] and references therein), no solutions of equation (3.50) can exist on each FHS 
if the right hand side does not change sign on each FHS. The only exception to this rule 
is the solution Nt =constant for the particular case when the right hand side of (3.50) 
vanishes. Note that (the metric approximations of) equations (3.49), (3.50) also hold in 
GR for projections on an arbitrary spacelike hypersurface since 

G^(t)ll + G(t)ikh\^^ — 2^(i)xi- (3.51) 
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Here is the Einstein tensor family corresponding to the metric family (2.18). 

Equation (3.49) represents one field equation where is partially coupled to space- 
time geometry. This coupling has the particular property that it is does not depend 
exphcitly on the intrinsic curvature of the FHSs. We postulate a second set of field 
equations having this property, namely 

^{t)jl = G(t)jl = ^{t)j\i ~ ^t,j = i^T^^f^jl- (3.52) 

The field equations (3.49) and (3.52) superficially look like just a subset of the Einstein 
equations in ordinary GR. But one crucial difference is that (3.49) and (3.52) hold in 
GR for projections into arbitrary spacelike hypersurfaces, whereas they are required to 
hold only for projections into the FHSs in quasi- metric theory. Anyway the full coupling 
between the Einstein tensor and the stress-energy tensor does not exist in quasi-metric 
theory. One reason for this is that such a coupling would explicitly depend on the intrinsic 
curvature of the FHSs [2], [15], and NKE of the FHSs would no longer be possible. 

We now have a sufficient number of field equations to determine the unknown quan- 
tities in equation (2.18). In order to solve these field equations it is very convenient to 
have general expressions for the spatial curvature of the FHSs. Such expressions can be 
found directly from the spatial part of (2.18). We find 

H(t)ij = c~^(a^|fc - J^-f^)ht)i3 - - c'^aj^iaj^j, (3.53) 

where Hj is the family of spatial Einstein fields intrinsic to the FHSs. Moreover, con- 
tracting (3.53) we can find an expression for the family of Ricci scalar fields Pt intrinsic 
to the FHSs, namely 

= (dlip ~ ^'^"'''^1' + 2c-%^,a>. (3.54) 

We see that given initial conditions on some FHS in JJ", i.e., Nt, Nt,j_, and (N^),±, 
equations (3.52) are immediately determined by these data whereas equation (3.49) is 
not. This means that to determine the time evolution of from equation (3.49) one must 
in addition specify n{T(t)i± + T(^t)i) ^^^^ FHS. Thus equation (3.49) is the dynamical 
equation and equations (3.52) are constraint equations. Note that for illustrative purposes 
it is often useful to study systems where the condition 

N\,^3c-^ariN\ ^ Kt = 3^, (3.55) 

holds, and where the second equation follows from (3.9). The geometrical interpretation 
of the condition (3.55) is that the extrinsic curvature tensor family Kt is traceless. That 
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is, taking the trace of equation (3.10) and requiring that it vanishes yields equation (3.55). 
Note that if the condition (3.55) holds, it implies that it is not necessary to specify an 
initial value {N^),i to be able to solve (3.49) on an initial FHS. 

Equations (3.49) and (3.52) represent one dynamical equation and three constraint 
equations. Compared to e.g. GR, where one has six dynamical equations and four 
constraint equations, we see that the field equation (3.49) yields the dynamical evolution 
of one scalar field rather than of the full unconstrained extrinsic curvature tensor, as 
in GR. Besides, due to the contracted Bianchi identities the local conservation laws 
are identities in GR. In quasi-metric theory, however, the local conservation laws in AT 
represent real constraints. Thus quasi-metric theory is quite different from GR when it 
comes to the way matter is explicitly coupled to the space-time geometry; in this respect 
quasi-metric gravity is more similar to Newtonian theory than to GR. 

In the next section we will see how both KE and NKE are realized in 77 by formulating 
the dynamical theory as an initial value problem. 



3.3 The quasi-metric initial value problem 

Now we want to use the results of the previous section to carry out a further investigation 
of the dynamical theory in J\f, which is not a metric manifold. We do this by aiming 
at a formulation of the theory as an initial value problem for the evolution with time 
of the geometry of FHSs; in contrast to the standard metrical initial value problem our 
non-metric counterpart should involve both KE and NKE of the FHSs. The dynamical 
theory can be extended to be valid in jV by adding extra effects coming from the NKE 
as ilhistrated in section 2.3. 

To begin with we split up the time evolutions of Nt and into kinematical and 
non-kinematical parts. We thus define 

'cN? Nt 

where c~^Xj and c~^yt respectively, define the KE and the NKE of Nf with time. Fur- 
thermore the time evolution of h(t)ij is given from equations (3.11) and (3.12), namely 

On the other hand we define the time evolution of the spatial scale factor to be 

_2- , -lfy_ 1 [1 I Nt,tl Nt,l /Q cro\ 



c '^Xt + c V=^f72 (3-56) 
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where c '^Xt represents the KE and c ^Ht represents the NKE of the spatial scale factor. 
Prom equations (3.56) and (3.58) one easily finds that 

yt^Ht- ^, (3.59) 



and 



cN, 



1 Ni 



= - c-^Ht + 7 + ^ ■ (3-60) 



t Ni 



To have a well-defined initial value problem in Af (for the case when there is no coupling 
to non-gravitational fields) it remains to define the specific form of the NKCE. This will 
be of the type (2.14). To find the right hand side of the NKCE wc first take the particular 
case of a static vacuum in the metric approximation, using a GTCS where the shift vector 
field is required to vanish. (We show in section 4.1 that such a solution can be found.) 
Then, from the metric approximations of equations (3.49) and (3.54) we find 

p 

c'^a'-pUj^i = — . (3.61) 

For this particular case we expect that all the information available for constructing the 
NKCE should be contained in the 4-acceleration field via (3.61). Besides, a second 
particular case is the empty cosmological toy model found in section 4. For this case 
the field equations are vacuous but we can still construct the curvature scalar ^ = 
Thus it may seem natural to construct the right hand side of the NKCE from the scalar 
Y and c. But from equations (3.49) and (3.54) we see that that Pt may take negative 
values inside matter, and since the right hand side of the NKCE should be non-negative 
everywhere we cannot use this scalar explicitly. However, we can try to combine the 
expressions obtained from the two particular cases into one general expression. That is, 
we define the NKCE as the sum of one global term coming from the global curvature of 
space and one term coming from the local gravitational field, i.e.. 



Ht=cFf'jr^,Ft - c-'xt = ^ry- + c- V«Wi, (3.62) 

where Ft — NfCt are the local scale factors. From (3.62) we see that Ht is indeed non- 
negative everywhere. Besides (3.62) becomes vacuous in the Newtonian limit as it should. 

Now we are ready to formulate the quasi-metric initial value problem in JJ" (in absence 
of coupling to non-gravitational fields, see below). This is equivalent to formulating an 
evolution scheme for the geometry ht (or equivalently, for Nt). This may be done as 
follows: On an initial FHS we choose initial data Nt, ^t:l, Nt,t, N\ {N"^),± and the 
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various projections of Tj such that the constraint equations (3.52) hold. Note that the 
chosen equation of state characterizing the matter source on the initial FHS must be 
compatible with the initial data. Also note that Ht and yt are explicitly given from the 
initial data on the initial FHS by respectively equations (3.62) and (3.59), and that c~^Xt 
is given on the initial FHS from the initial data via equation (3.60). 

From the initial data we are able to find Nt on the deformed FHS. Moreover, from the 
local conservation laws (3.45), (3.46), (3.47) and (3.48) we are able to find K>T(f^^f-^i]_ and 
oil the deformed FHS. (However, as mentioned earlier, in special situations one 
expects that the global NKE implies local non-conservation. For such cases equations 
(3.47) and (3.48) do not hold, and how much they are violated must be specified in 
addition.) The remaining projections of on the deformed FHS must be inferred from 
the equation of state. Also the shift vector field iV* on the deformed FHS must be 
compatible with equations (3.52) on the deformed FHS. But in order to predict Nt,j_ 
(and Nt,t) on the deformed FHS, or equivalently Nt on yet another deformed FHS, we 
need to solve the dynamical equation (3.49) on the initial FHS. To do this it is essential 
to supply the quantity 2^(i-)Xl = '^(^(t)ll on the initial FHS. When we have done 

this we are able to advance one time step of the evolution process. We can then repeat 
the above described scheme for each subsequent time step. 

But in general we must also have an evolution scheme in J\f. The reason for this 
is that when there is coupling to non-gravitational fields (e.g., electromagnetism), non- 
gravitational force laws in J\f determine the motion of sources and thus T^, which enters 
into the evolution scheme in JJ". So when there is couphng to non-gravitational fields the 
evolution schemes in jV" and jV are intertwined. 

To have an evolution scheme in jV it is necessary to find yt and hj^ for each step of 
the evolution (see equation (2.23)). Since hjr is intended to represent a generahzation of 
"the local coordinate vector from the center of gravity" which does have meaning for the 
spherically symmetric case, it is natural to seek for an equation which only solution is 
bjF — f-§;: for the spherically symmetric case. Furthermore the wanted equation should 
be linear in hjr to ensure unique solutions, and it should involve a^F since one expects bjF 
to vanish whenever a^F does. By inspection of the spherically symmetric case it turns out 
that it is possible to find an equation which has all the desired properties, namely 

a^lj. + c'^ajr^a'^ 6^ - a^|^ + c'^ajr^a^ - 2a^ = 0. (3.63) 

Since we may calculate b^F from equation (3.63) and then yt for each step of the evolution 
according to the above described evolution scheme, we are also able to construct h^, N 
and for each time step along the lines shown in section 2.3, using equations (2.24), 
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(2.25) and (2.26). Thus with the help of the coupled evolution schemes in U and M we 
are able to find both metric families gj and g^. Prom the relevant equations we see that 
the quasi-metric evolution scheme contains both KE and NKE. The crucial difference 
between the NKE and the KE of the FHSs is that the former can be found on subsequent 
FHSs directly from the spatial metric via the NKCE whereas the latter depends on the 
field R{t)ii_ which must be specified separately for each subsequent FHS. 

Now we must emphasize an important point: since is the active stress-energy 
tensor and thus not directly measurable locally, one must know how it relates to the 
passive stress-energy tensor % in A/". The relationship between Tt and % depends in 
principle on the general nature of the matter source. To illustrate this, consider and 
% for a perfect fluid: 

Tt = (pm + c~^p)ut(g)Ut + pgt, %= {pm + c~'^p)\it®nt + pgt, (3.64) 

where pm is the active mass-energy density in the local rest frame of the fluid and pm is 
the passive mass-energy density measured in the local rest frame of the fluid (p and p 
are the associated pressures). Also, in (3.64) we have used the definitions Pm\/^t=Pm\/ht 
and p\fjit=p^/7h, where ht and ht are the determinants of the spatial metric families 
and ht, respectively. Now the relationship between p^ and pm is given by 

^N^^Pm for a fiuid of material point particles, 

(3.65) 

Pm for a null fiuid, 

and a similar relationship exists between p and p. The reason why the relationship 
between p^ and p^ is different for a null fiuid than for other perfect fiuid sources is that 
spectral shifts of null particles influence their passive mass-energy but not their active 
mass-energy. On the other hand the global NKE does not affect the passive mass-energy 
of material point particles (see section 4.2). Finally, note that more "physically correct" 
local conservation laws than those shown in equations (3.42), (3.43) and (3.44) may be 
found in jV by calculating V-% when gt is known. But these local "physical" conservation 
laws for passive mass-energy do not take any predetermined form. 



3.4 Equations of motion 

In section 2 we derived the equations of motion (2.10) valid for arbitrary test particles. 
However these equations are of no practical value as long as we do not know the rela- 
tionship between the 4-acceleration field at and the degenerate 4- acceleration field a. In 
this section we show that in fact a=at. 
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For starters it is convenient to introduce the difference in coordinate 3-velocity be- 
tween a test particle moving along an arbitrary curve and that of the local FO. We call 
this quantity Wj, and we write its square as w'^. Thus we define (using a GTCS) 

^'>-^Aw- (3.66, 

where Tt is the proper time measured along the curve and 7 is the time dilation factor 
relating the proper time intervals along the curve to the proper time intervals of the local 
FOs. Since Wj is required to be a vector field family tangent to the FHSs, the space-time 
scalar product of Wt and must vanish. That is, if we view Wt as a family of 4- vectors 
we must have 

= 0, w'^t)^(^t)ix = ^(t)M^ft) = 0, =^ W(t)o = jN'w(^t)i- (3.67) 

Note that since the explicit dependence on t of any scalar gravitational quantity is de- 
termined from its dimensionality when measured in atomic units, we must have 

d d 1 

—w^={w(^t)nw'^t))^t^{w(t)iW(^t)kh',})),t^O, Ql^(t)j ^ j'^{t)j- (3-68) 

This means that we must have 

^ — * i ^0 7 _^0 j /o a(\\ 

W{t)j = -W(^to)j=j-^Wj, wl^^ = jwl^^^=jw'. (3.69) 

We then split up the 4-velocity field along the curve into tangential and normal pieces. 
This is done as follows: 

From equation (2.9) we see that to find an expression for a we need to calculate the 
quantity V a. Uf. We will use equation (3.70) for this purpose. By direct calculation using 

dt 

equations (3.67), (3.68) and (3.69) we find 

V^«^{t),.=VA<) =0, =^ VaWj^O, (3.71) 

at ^ at y-' at 

and furthermore we calculate 

Van(i)^=Vanf,)=0, ^ Van, = 0, (3.72) 

confirming the requirement (2.3). But from equations (2.9), (3.70), (3.71) and (3.72) we 
get 



at 



Va^ut = 0, a = at, (3.73) 
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as asserted previously. Using a GTCS we then get the final form of the equations of 
motion, namely 

(Px^ /*,, dt dx^\dx'' /cdTi-\2 , ,, , 

Equations (3.74) are valid for both timelike and null curves. For timelike curves one may 
suitably choose a parameterization in terms of the proper time measured along the curve, 
i.e. one may choose A = cr^. For null curves A must be chosen to be some other affine 
parameter such that ^ = 0. Note that these equations do noi violate LLI. See reference 
[15] for further comments. 

4 Spherically symmetric vacua 

In this section we solve the field equations for two metrically static, spherically symmetric 
vacuum cases. (We use the term "metrically static" to refer to cases where Nf does not 
depend on x^ and Nj vanishes in some GTCS.) 

The first case we analyze is the rather trivial one where space-time is isotropic and 
globally free of matter, i.e. a.jr vanishes identically. The field equations (3.52) are vacuous 
in this case whereas equation (3.49) has the solution Nt =constant. However, from 
equations (3.53) and (3.54) we get (where is the family of Ricci tensors corresponding 
to the spatial metrics hj) 

Ptj 1 T Pt _ 1 

{Ntcty 6 {Ntcty 

or equivalently 



Pf- 

Pm = -^ht)iv (4-2) 

which means that the solution is an Einstein space. It is easy to show that is the 
only compact spherically symmetric solution of (4.2) with Pt > 0. Note that Nt is just 
an arbitrary constant, and that we may set iV^ = 1 by making some appropriate choice 
of scale for the spatial coordinates (and simultaneously scale x° and t). This means that 
each space-time metric of the family gj is of the type S'^xR. Wc can find a relationship 
between these metrics in terms of the epoch t from the NKCE (3.62). This yields Ht = \ 
which means that c^^Xt = yt = from (3.59) and (3.60). Besides this yields Vt = from 
equation (2.23). Thus the FOs in AT do not move with respect to the FOs in J\f, i.e. the 
metric families and g^ are identical. We may then write 

N^Nt^l Ht^^. (4.3) 

V 
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The one-parameter families of S^xR space-time metrics consistent with (4.3) take the 



and we notice that (4.4) is of the form (2.18), as required. Furthermore we see (using 
equation (3.33)) that (4.4) yields that active mass rrit increases with time. We may 
express rrit and by and "if to valid at the arbitrary epoch to, getting 



As usual when rrit is variable this is because by definition, gravitational quantities get 
an extra variability when measured in atomic units. In this particular case it means 
that gravitational length scales measured in atomic units (and vice versa) are different 
at different epochs t. Note that it is just as vahd to interpret (4.4) to mean that atomic 
length units are shrinking as the usual interpretation in terms of an increase of spatial 
dimensions with time. This illustrates that the family of line elements (4.4) expresses the 
gravitational scale (represented by the spatial coordinates in a GTCS) in terms of units 
determined from the atomic scale. Moreover, since the gravitational scale is represented 
by coordinates belonging to a particular class of coordinate systems, this means that no 
particular intrinsic length scales can be associated with gravitation. In contrast to this 
non-gravitational interactions define intrinsic length (and mass) scales. 

It is important to realize that the metrics (4.4) arc vacuum metrics, since we have 
neglected matter creation as mentioned in the previous section. Thus (4.4) is just an 
approximate, not realistic toy cosmological model. However we may insert Ht = j = 
Hq for the value of the Hubble parameter for the present epoch into (4.3) to get a 
realistic value for the age of the Universe. This is possible since in contrast to standard 
cosmology, one should not expect that inclusion of matter into the toy model must 
necessarily result in "braking" of the expansion. In fact one may show that the only 
possible cosmological models allowed with a perfectly isotropic matter distribution are 
models where the source is a null fluid. That is, according to our theory no cosmological 
models with an isotropically distributed material source is possible, see appendix C. This 
indicates that matter density has nothing to do with the time evolution of the cosmic 
expansion. Thus the age Ht^ should be taken as a realistic value and not merely as 
an upper bound. Also one may not neglect time dilation effects in the quantity 
when interpreting redshift measurements. We return to the question of interpretation of 
spectral shifts in section 4.2. 



form (with t — 




(4.4) 



t 



(4.5) 
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4.1 The conformally Minkowski case 

The globally matter-free metric family (4.4) found in the previous section is a natural 
candidate for a cosmological toy model since it is isotropic and thus spherically symmetric 
about every point. On the other hand we need a model for the gravitational field exterior 
to a spherically symmetric body in an otherwise empty Universe. In principle we must 
then find a solution of the field equations which is spherically symmetric about one 
particular point only and where the boundary conditions for this solution should be 
S^xR, i.e. the metric family should be of the form (4.4) at the boundary. But such 
a solution does not exist according to equation (3.50) and the subsequent discussion. 
However, it is possible to find a solution if we relax the boundary conditions, see [14] . To 
avoid the problems at the boundary, in the case where we can neglect the NKE due to a 
global scale factor we shall feel free to use a static metric approximation, i.e., we assume 
that the each metric of the metric family is identical and asymptotically Minkowski. 

Prom (4.4) we see that we get a Minkowski background if we let to^oo and set t — to. 
We may then set Nt{oo) — 1. This means that (using (3.33)) 

^t{r) = ^t{oo)N-\r), (4.6) 

since the acceleration field is purely radial. Equation (4.6) means that reference units 
at infinity are being used for each metric gGgf. However, taking these units as those 
obtained in the limit t^oo in (3.33) we see that they get singular if one compares to the 
units determined by a S^xR background. That is, \E'f vanishes in the limit t^oo. But 
the field equations and equations (3.53), (3.54) still work fine for metric approximations. 

Even if asymptotically Minkowski solutions in principle are the result of using sin- 
gular reference units in the field equations (3.49), (3.52) and in equations (3.53), (3.54), 
the discrepancy between these solutions and those with acceptable boundary conditions 
should be small to a good approximation for many situations. So we feel free to find an 
asymptotically Minkowski, spherically symmetric static solution of (3.49), and apply it 
in appropriate circumstances. Proceeding with this, firstly we notice that the spherically 
symmetric condition makes equations (3.52) vacuous. Secondly, the field equation (3.49) 
together with equation (3.53) give 

H{t)ij + c''^aj:i\j c~^ajFjajrj c'~^ajrka%h{t)ij = 0, a%^^ + c~'^a%aj^k = 0. (4.7) 

When we apply these equations to the general spherically static line element (2.21) we 
get the unique solution 
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where, analogous to the Schwarzschild solution, rs=^^. Some details of the derivation 
of (4.8) are given in an appendix. 

Next we calculate the spatial Ricci scalar P corresponding to (4.8). Using (A. 4) we 
find 

where p is an isotropic radial coordinate defined in (A. 14). Furthermore, from the NKCE 
(3.62) in combination with equation (3.61) valid for this particular case we find 

As can be readily seen, (4.8) contains no coordinate singularity, thus there is no event 
horizon. This is in accordance with the requirement that a global time coordinate must 
exist. Furthermore, there is a curvature singularity at the spatial origin. However this 
singularity is inaccessible due to the fact that the NKE in the tangent space increases 
towards the velocity of light at the spatial origin. Using (2.20) and (4.10) we find 

and from (4.11) we see that v{r)^c in the limit r^O. As can be seen from equation 
(4.13) below, this has the effect of removing the curvature singularity present in the 
intrinsic geometry of the FHSs when one constructs the geometry h from the geometry 
h as described in section 2.3. And although the space-time metric (4.13) is still singular 
at the spatial origin, this is merely a consequence of the fact that the spatial origin 
is inaccessible to material particles and thus that the premise of a gravitational point 
source is unphysical. We interpret this to mean that gravitational collapse cannot, even 
in principle, form physical singularities when the NKE is taken into account. 

Now we must construct the metric g from the solution metric g shown in (4.8). We 
find g by the method discussed in section 2.3. Before we do this we note that equations 
(3.59) and (3.60) yield 

c-^xt = -yt = -Hu c-^xut = Hut = 0. (4.12) 

Thus we find that the construction of the metric g is consistent with the equations of 
evolution. 

We now use equations (2.22), (4.8) and (4.11) to construct the metric g and this is 
sufficient to check if the "classical" solar system tests come in correctly. We find the 
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unique metric 

= -(l-7 + |i + •■■)(<^-V + (l + ^ + ^ + ■■•K + '■''^^i^ (4-13) 

which agrees with the Schwarzschild metric to sufficient order [1]. Thus the metric (4.13) 
is consistent with the four "classical" solar system tests (radar time delay experiments 
included), and our theory seems to have the correct metric correspondence. (This can 
be directly confirmed by applying the equations of motion (3.74) to the metric (4.13).) 
We note that since (4.8) is conformally flat as can be seen from (A. 15), it has vanishing 
Weyl curvature. But since the Schwarzschild metric has vanishing Ricci curvature, it 
is the Weyl curvature of that metric that makes GR agree with the solar system tests. 
This is in contrast to our theory, where the agreement with the solar system tests at the 
post-Newtonian level cannot possibly come about were it not for the fact that the NKE 
has been integrated into the metric structure. Thus the existence of an implicitly coupled 
dynamical degree of freedom is essential in order to make quasi-metric gravity viable. 

4.2 Testing the NKE-paradigm 

To determine whether or not our theory is viable, we must try to confront its predictions 
with data from modern relativistic gravitational experiments. But first we must show 
that it has a consistent Newtonian limit. This limit can be defined in such a way that 
the FHSs are taken as flat and static. Thus all the NKE is ignored in the Newtonian 
limit and we can drop the bar labels where this is convenient. In a static GTCS with 
Cartesian coordinates we have the usual Newtonian form of the metrics 

9{t)oo = g{t)m = -1 + 2c~'^U{x^'), g(^t)io = 9{t)io = 0, h(^t)ij = h(t)ij = Sij, (4.14) 

where —U{x^) is the Newtonian potential. Note that (4.14) is consistent with the gen- 
eral metric family (2.18) since to Newtonian accuracy, we can neglect any contribution 
to h(t)ij of order c~'^U{x^). In appendix B we use (4.14) to construct a well-defined New- 
tonian limit for our theory by calculating the Newtonian approximation to the evolution 
equations (3.44) applied to a perfect fluid. Besides we show that the field equation (3.49) 
reduces to Newton's gravitational equation in the Newtonian limit and that (4.14) is 
consistent with the equations of motion (3.74) with vanishing non-metric terms. 
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We note that the usual Newtonian assumptions of weak gravitational fields and slow 
motions include the condition that the vector fields defined in equation (2.23) yield 
neghgible corrections to the metrics (4.14). Also we must also assume that the system 
is "small" if we wish to have Newtonian theory as a limit of our quasi-metric theory; 
assuming a small system is essential to ensure that effects from the "true" non-kinematical 
expansion coming from a global scale factor are neghgible. The reason for this is that 
in the Newtonian approximation the FOs cannot be defined unambiguously, so we may 
consider any set of observers static with respect to some arbitrary Galilean system as 
FOs. This is a good approximation for small systems and nearly flat FHSs. but for large 
systems one cannot neglect the "motion" of distant FOs with respect to each other due 
to the NKE coming from a global scale factor. This needs not be small even if the change 
of the global scale factor is small everywhere. 

Effects related to that part of the NKE which is integrated into the metric structure of 
gf via Vf we refer to as local effects of the NKE, whereas the remainder of the NKE yields 
the global effects. In essence Newtonian theory should correspond with the metric part 
of our theory, and not with the non-metric part, which has no Newtonian counterpart. 
Thus even if the local effects of the NKE can be neglected at the Newtonian level of 
precision, there is no reason to think that this also applies to the global effects of the 
NKE. 

Most new metric theories of gravity having a satisfactory Newtonian limit can be 
analyzed within the parameterized post-Newtonian (PPN) formalism. This is constructed 
to handle the next higher level of precision beyond the Newtonian approximation for 
systems with weak fields and slow motions [1]. The result is a set of PPN-parameters 
which must take a set of specific values if the metric theory is to be regarded as viable. 
However, quasi-metric theories are not well-suited for a standard PPN-analysis. One 
reason for this is that the global effects of the NKE are missed since any existing PPN- 
formalism does not account for such effects. This means that a PPN-analysis of quasi- 
metric theory is limited to its metric approximations. But even if metric approximations 
to quasi-metric theory can be PPN-analyzed in principle there is another complication; 
namely that the /oca/ effects of the NKE must be taken into account at the post- Newtonian 
level via the construction of g from g. That is, any PPN-analysis applied to quasi- 
metric theory would result in a post-Newtonian metric g and not in its counterpart g. 
This means that any naive attempt to PPN-analyze metric approximations to the field 
equations (3.49), (3.52) and equation (2.18) should not result in an acceptable set of PPN- 
parameters according to metric theory. A particular case confirming this expectation is 
the conformally fiat solution (4.8) which disagrees with the accepted values for two of the 
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PPN-parameters. It is only when the NKE corrections are included that we get a metric 
with the standard values of these parameters, i.e., the metric (4.13). 

Consequently, to analyze our field equations within some sort of weak field formalism 
one might first try to find the PPN-metric g (which does not have an acceptable set 
of PPN-parameters), and then construct a second PPN-metric g from g and Vj in the 
usual way. But for this to have any chance to work one cannot assume (as is done in the 
standard PPN formalism) that all dynamical degrees of freedom are exphcitly coupled 
to T(. Taking into account all the complications, to construct a weak field expansion 
suitable for quasi-metric theory would seem to be a rather messy task but obviously a 
subject for further work. 

Even without a standard formalism within which the theory can be compared to ex- 
periment, we can try to pinpoint some of the "direct" predictions due to the hypothetical 
existence of the global NKE. The fact that the NKE takes the form of an "expansion" 
suggests that one of the obvious effects to look for is a corresponding non-kinematical 
redshift (NKR). As first noted by Synge, spectral shifts in GR (and other metric theories) 
can be analyzed within a unified formalism [11]. This means that gravitational, Doppler 
and standard cosmological spectral shifts are of the same nature; we refer to such spectral 
shifts as kinematical. Kinematical spectral shifts may be blueshifts or redshifts depend- 
ing on the specific circumstances, whereas the NKR is always a redshift. It should be 
simple to extend the unified formalism described in [11] to include the NKR as well by 
including the non-metric piece of the connection. This should be so because the NKR 
essentially arises from the non-metric part of the connection just as kinematical spectral 
shifts arise from its metric part. 

As a straightforward example, an easy calculation of the NKR can be executed for 
the toy cosmological model (4.4). With a change to a convenient radial coordinate x 
(4.4) reads 

dsl = -{dx^'f + c'f (dx^ + sin^xdO') . (4.15) 

We want to find the equation of a null geodesic in the x direction by solving (3.74) for the 
metrics (4.15). To begin with we find that all the relevant ordinary connection coefficients 
vanish, whereas the only relevant non- vanishing degenerate connection coefficient is V}^ = 
t~^. Furthermore the time component of equation (3.74) yields ^ ^constant whereas 
the radial component yields (using t as a time parameter) 
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The solution of (4.16), subject to the condition = (ct) (which follows from 

equations (3.66) and (3.69)), is 

x(t)=Xo + lnf. (4.17) 

Now consider two FOs with radial coordinates Xo and Xi respectively, and with identical 
angular coordinates. At epoch to two pulses of light are emitted in the ^-direction by 
the first FO. The pulses are separated by a short time interval Ato- At epoch ti the light 
pulses are received by the second FO, and the separation time interval is Ati. Using 
(4.17) we find (neglecting terms of higher order in small quantities) 



Ato to' 



(4.18) 



a result identical to that of metric cosmological models, i.e., representing a redshift of 
momentum for decoupled massless particle species. Also, we can find the speed w=Vw^ — 
ct{-^){t) = c'^t{^)(t) with respect to the FOs for some inertial material particle emitted 
in the radial direction. Since {^){t) = ^(^)(to) (which follows from (3.69)) we have 
w =constant, i.e., material particles do not slow down with respect to the FOs. This 
result is different from that of standard cosmology, and illustrates the different nature of 
the NKE compared to the KE. The curves traced out by material inertial test particles 
moving in the radial direction can be found by solving (4.16) subject to the condition 
(^)(t) = c~'^wt~^, and the solution is 

x(t)=Xo + -lnf. (4.19) 
c to 

The calculation leading to (4.18) can be generahzed, i.e., one can in general find the 
equations of null geodesies connecting pairs of FOs and then calculate the frequency shifts 
of fight pulses. The result is a superposition of NKR and kinematical spectral shifts. In 
particular, radiation emitted from a FO in a gravitational well should, as seen from a FO 
far away, show both the gravitational redshift and the NKR. A natural question is how 
this prediction of the theory relates to experiments such as measurements of spectral 
line shifts in the solar atmosphere and the Pound-Rebka-Snider measurements of the 
gravitational spectral shift at the Earth's surface [1]. The answer is that for these cases 
the metric approximation (4.13) represents the gravitational field to sufficient accuracy, 
and (4.13) predicts no NKR. (Even if we replace the Minkowski background present in the 
metric approximation with a x R-background the extra NKR due to this background 
has only negligible infiuence on the above-mentioned experiments, see reference [14].) 
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Another direct consequence of the hypothetical existence of the NKE is that ex- 
tended systems held together by contact forces should experience extra tidal stress due 
to geodesic deviation. An important special case is that of tidal forces experienced by 
objects moving in the gravitational field of the Earth. By using (4.13) we can find the 
corresponding tidal metric field. The difference between this and the analogous tidal 

2 2 

field of the Schwarzschild metric is of the order whereas the dominant contribution 

2 

to each tidal field is of order Thus it does not seem feasible to directly detect the 
NKE by measuring tidal stresses. 

Although it would seem difficult to detect the NKE directly by doing local experi- 
ments, its relevance for explaining astrophysical observations from first principles is quite 
promising. To have any hope of doing better than metric theories in this respect one must 
focus on situations where the metric approximation is not sufficient. In particular, astro- 
physical systems which are large compared to the solar system, but yet non-relativistic in 
the sense of small velocities and weak gravitational fields should be well-suited for testing 
the non-metric aspects of our theory. The most obvious examples of such systems are 
galaxies. Thus one may have some hope that the NKE paradigm may possibly explain 
observed galactic kinematics from first principles, a feat that metric theories seem inca- 
pable of doing [12]. We return to this issue in a followup paper [14], where it is shown 
that the NKE may also have been detected in the solar system. 

More exotic systems which have been modelled successfully by applying GR are the 
binary pulsars [1] . In particular, within GR the decay of orbital periods of such objects 
is naturally explained by assuming that the systems emit gravitational radiation. To 
have any hope of competing, this phenomenon must have a natural explanation within 
our quasi-metric theory also. Of course, to test if this is the case nothing short of a 
quantitative analysis is sufficient. 

The fields where the difference in philosophical basis between our quasi-metric theory 
and metric theories shows up most clearly are cosmology and large-scale astronomy. That 
is, the traditional big bang models (and alternative models also) describe the Universe 
as a dynamic system which evolution can be deduced from field equations given the 
initial conditions. Using this approach it is not surprising that one of the holy grails of 
modern cosmology is to deduce the "correct" big bang model from observational data. 
In contrast to this, according to our theory the NKE due to a global scale factor is the 
dominant phenomenon on cosmological scales, and thus our description of the Universe 
is not primary that of a dynamic system. In fact, a natural beginning of the Universe 
according to our theory is described by the empty Universe toy model (4.4) in the limit 
t^O. Thus even if this model has a curvature singularity at i = it does not represent 
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a physical singularity since no particle world line actually originates at the singularity. 
This removes the vexing question of specifying "natural" initial conditions, a question 
which cannot be avoided in a purely dynamic approach to cosmology. Instead we get 
the problem of constructing a quantitative theory of particle creation based on the NKE 
and vacuum fluctuations of non-gravitational fields. This may perhaps be done along the 
hues described in reference [10]. Thus the possibility exists that one may construct a 
realistic universe model where the Universe has "expanded from a hot dense phase" but 
yet does not actually belong to the class of big bang models. 

It remains to be seen if non-kinematically expanding universe models are capable 
of giving natural explanations of the plethora of astronomical observations which exist 
today; examples are helium abundances, the structure of the microwave background etc. 
In this context one must be aware that interpretations of observational data may be very 
different for a non-kincmatical Universe model than for a big bang model. The main 
reason for this is that the nature of the expansion is different for the two model types. 
In particular this means that in a non-kinematical model one cannot uncritically take 
observed rcdshifts of objects as indicators of distance. That is, any observed spectral shift 
includes not only the NKR due to the expansion on large scales of the Universe, but also 
the NKR due to the intrinsic expansion of gravitational fields. (See [14] for an idealized 
model of an expanding gravitational field.) If the latter is not negligible compared to the 
former, interpretation within standard cosmology may lead to peculiar results. And for a 
sufficiently strong gravitational field the NKR due to the intrinsic expansion will not be 
negligible since the global part of the NKE goes as That is, for an observer far away 

far down into the gravitational well appears larger than the local value of due 
to gravitational time dilation. Thus, if an observer far away registers radiation escaping 
from the gravitational well, this radiation gets an extra NKR mimicking the cosmological 
redshift. Thus it is at least possible that massive expanding objects may systematically 
show excess redshift which should not necessarily be interpreted as cosmological. A well- 
known example of this is the "redshift controversy" regarding the nature of quasi-stellar 
objects [13]. Also measurements of the Hubble parameter may be sensitive to the NKE 
intrinsic to e.g. galaxy clusters. 

5 Conclusion 

In this paper we have introduced the so-called "quasi-metric" framework; this is a non- 
metric space-time framework based on the NKE concept. Furthermore we have con- 
structed a theory of gravity compatible with this framework. The vital question remains; 
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is the theory viable according to basic theoretical and experimental criteria? 

At least our theory is relativistic, and it displays the correct Newtonian limit for 
a large class of situations. However, for "large" systems non-metric effects may not 
be neglected, and for such systems Newtonian theory does not correspond with our 
theory On a higher level of precision our theory corresponds with GR through metric 
approximations; we have illustrated this by constructing one for the static spherically 
symmetric case. This means that our theory should agree with GR for the classical solar 
system tests (to post-Newtonian accuracy). However, due to the fact that our theory is 
not suitable for a standard PPN analysis, a formal weak field comparison to GR is not 
available, see the previous section. 

The theoretical structure of our theory makes it plausible that its equations are well- 
posed and self consistent. This also indicates that its predictions are unique, i.e., that 
they do not depend on the method of calculation. Finally the theory should not be 
regarded as complete, i.e. all aspects of it are not yet fully developed. 

So we do not yet know if our theory is viable or not. However, there should be a 
fair chance that the NKE concept may be relevant for a better understanding from first 
principles of the behaviour of the physical world. This should inspire to further work 
along the direction we have laid out in this paper. 

References 

[I] CM. Will, Theory and Experiment in Gravitational Physics, Cambridge University 
Press (1993). 

[2] K. Kuchaf, Journ. Math. Phys. 17, 792 (1976). 

[3] C. Gerhardt, Commun. Math. Phys. 89, 523 (1983). 

[4] C.W. Misner, K.S. Thorne, J.A. Wheeler, Gravitation, W.H. Freeman & Co. (1973). 
[5] C.H. Brans and R.H. Dicke, Phys. Rev. 124, 925 (1961). 
[6] R.H. Dicke, Phys. Rev. 125, 2163 (1962). 

[7] J. Isenberg and J. Nester, Canonical Gravity, in General Relativity and Gravitation, 

Vol. 7, Editor A. Held, Plenum Press (1980). 
[8] K. Kuchaf, Journ. Math. Phys. 17, 777 (1976). 

[9] S. Weinberg, Gravitation and Cosmology, John Wiley & Sons, Inc. (1972). 

[10] L. Parker, Particle Creation by the Expansion of the Universe, in Handbook of 
Astronomy, Astrophysics and Geophysics, Vol. II Galaxies and Cosmology, 
editors V.M. Canuto and E.G. Elmegreen, Gordon & Breach (1988). 

[II] J.V. Narlikar, Am. J. Phys. 62, 903 (1994). 

[12] V.V. Zhytnikov and J.M. Nester, Phys. Rev. Lett. 73, 2950 (1994). 



93 



[13] H. Arp, Quasars, Redshifts and Controversies, Interstellar Media (1987). 

[14] D. 0stvang, chapter 3 of this thesis (2001). 

[15] D. 0stvang, chapter 1 of this thesis (2001). 

[16] J. Isenberg, Class. Quantum Grav. 12, 2249 (1993). 

A Derivation of a static line element 

In this appendix we give some details of the derivation of the solution (4.8). To do this, 
we insert the general hne element (2.21) into (4.7) to determine the functions A{r) and 
B{r). Connection coefficients for the line element (2.21) can be found in e.g. [9]. The 
components of the acceleration field are (using the notation '=^) 

B' 

a-Tr = c^^' = aj:^ = 0, (A.l) 

and the spatial covariant derivative of the acceleration field has the non-zero components 



c 



B" ,B\o A'B' 



2 



B'r 

sin~^6'a^0l0 = ajr^^g = (A.2) 



B ^B' 2AB 
The spatial covariant divergence of the acceleration field is 

_ 2/ B" B'^ A'B' _i B' \ 
"^l'^ " H2AB ~ 2AB^ ~ 4l^ + ABJ- ^^-^^ 

Next a standard calculation yields the spatial Ricci tensor P for the line element (2.21), 
which we use to find the spatial Einstein tensor H. The result is 



r L A^. 
P 



(A.4) 



- P rA! 
sm~'^eH^^ = Hoe = Pee - -^hee = (^-5) 



Inserting the above relations into equations (4.7) yields 

B" A'B' 



B" B'^ A'B' 2B' _ 
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We can separate the variables A{r) and B{r) by inserting equation (A. 8) into equations 
(A.6) and (A.7). The result is 

(A.9) 
(A. 10) 



A" 

T " 


3A'2 


+ r-^A'{3A- 




= 0, 


B" 




B'^ 
— r — = — h r 
AB^ 


B 


= 0, 




A^ 


1 -r^ + r^ 
t> 


452' 





(A.11) 

One way of solving these equations is to introduce C{r)=r~. One may then show that 
(A. 10) is integrable in this new variable. This yields B[r), and (A. 11) gives A{r) directly 
when B[r) is known. The result is the metric (4.8), which is unique. An alternative way 
of writing (4.8) can be found by using the identity 



That (4.8) is indeed conformally flat, we can show by introducing an isotropic radial 
coordinate p replacing the Schwarzschild coordinate r: 



p^exp(/r-' [f^^ - ^] [l + (^)T'"*) ■ (A-13) 
Evaluating the integral and choosing an appropriate constant of integration, we find 



and (4.8) expressed in isotropic coordinates reads 

di' = (l - ^) - {dxy + dp" + p^dn^^ , (A. 15) 

which is conformally flat. 

B The Newtonian limit 

To find Newtonian approximations to the local conservation laws (3.44) in jV" or equiva- 
lently, to the evolution laws (3.45) and (3.46); in a static GTCS we can use the Newtonian 
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form (4.14) of each metric with Cartesian coordinates, i.e., ^oo~ — 1 + 2c~'^U, hij^Sij, 
where —U is the Newtonian potential. Furthermore, to Newtonian order the active and 
passive stress-energy tensors are equal. For a perfect fluid we then have, to Newtonian 
order 

^°° = p„c^ T'j = p^cwj, rj = pmw'wj+pd'j. (B.i) 

Also, to Newtonian accuracy c~^C/, c~'^w'^, c~'^p/pjn, c~^a^ajFi/ |g|o|^, c~^|a^|/|^| and 
|^^^^|/|g|o I are all small, of second order in the small quantity e, say. To Newtonian order 
one uses the approximations 

this is equivalent to retaining only terms of not higher than first order in e. The approx- 
imations made in (B.2) yield the same Newtonian limit as in the metric case. 

Keeping only lowest order terms and terms not higher than one extra order in e 
in (3.44) yields {t and the ordinary coordinate time are equivalent in the Newtonian 
approximation) 

d 

P:j+Pm{Wj,t+WjnW') = - p^ttj^j + 0{(?) = p^U ,j +0{(?) . (B.3) 

The equations (B.3) are identical to the Euler equations describing non-relativistic hy- 
drodynamics. Here we have used the fact that in the Newtonian limit aj^i — —U,i which 
follows from (4.14). 

Next, use the Newtonian form (4.14) of the metric in combination with (B.I) to get 
the Newtonian limit of the field equation (3.49). Neglecting terms of O(e^) equation 
(3.49) yields (using a Cartesian coordinate system) 

ATiGp^ = = -U\: (B.4) 

which is Newton's gravitational equation as asserted. Note that the field equations (3.52) 
become vacuous in the Newtonian limit. This also applies to the Newtonian limit of 
equations (3.53) and (3.54). 

Finally, to have a well-defined Newtonian limit we must show that our equations of 
motion (3.74) reduce to the Newtonian equations of motion in the Newtonian limit. But 
since equations (3.74) reduce to the usual geodesic equation in the metric approxima- 
tion (since t^oo), they certainly must reduce to the Newtonian equations of motion in 
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the Newtonian limit. With this result we have shown that our theory has a consistent 
Newtonian limit. 

When deciding whether or not it is sufficient to apply the Newtonian approximation 
for any particular system we should always have in mind that in our theory, the Newtonian 
limit is a strictly local limit. That is, even if the gravitational field of the system is weak 
everywhere and the motion of of its constituents is always non-relativistic we must also 
require that the system is "small" so that the global effects of the NKE may be neglected 
if Newtonian theory is to be applied with confidence. This is one of the crucial differences 
between our theory and any metric alternative. 

C Isotropic cosmological models with matter 

Rather than the vacuum toy cosmological model considered in section 4 we may try to 
construct non-vacuum toy models where the source is a perfect fiuid comoving with the 
FOs. That generalization of equation (4.4), each toy model may be described by 

a metric family gt of the form 

d^l = -N^{dx^y + {Ntctfidx^ + sm\dn^). (C.l) 

As a counterpart to the Priedmann models in GR it is reasonable to try to construct 
models where the matter distribution is isotropic. That is, we consider models where Nt 
does not depend on spatial coordinates. We show in this appendix that the only isotropic 
models possible in quasi-metric gravity are models where the source is a null fiuid. 

It is convenient to express the stress-energy tensor for a perfect fiuid comoving with the 
FOs via the coordinate volume density of active mass-energy and the corresponding 
pressure p, i.e. 

T^,, = T^,, = T*,,^(^y-p. (C.2) 

Note that pm does not depend on t if there is no particle creation (i.e., if the local 
conservation laws hold). Moreover, assuming no particle creation, the local conservation 
laws in 77 are sufficiently represented by equations (3.45) and (3.46), which yield 

J^HtPm = --^{pm - 3c~^p) = 0. (C.3) 

One reason why CutPm must vanish is that no other possibilities yield solutions of the 
field equations with a proper vacuum limit (i.e., the empty model found in section 4). 
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This means that a null fluid is the only nontrivial possibility satisfying (C.3). That is, in 
our theory no cosmological model with a perfectly isotropic material fluid source can exist. 
Thus in our theory, except for the vacuum model found in section 4, the only allowable 
cosmological models with isotropic FHSs are models where the source is a null fluid. 

To confirm that cosmological models with an isotropic null fluid source are indeed 
possible in quasi-metric gravity, we must solve the field equations. The field equation 
(3.49) yields, after some straightforward manipulations using equation (3.9) 

or equivalently, for a null fluid 

^lvr>''°^"t^3^"^'' ^ ^ ^ 

where the right hand side of (C.5) is a constant (for t constant). Integrating equation 
(C.5) twice and requiring continuity with the empty model found in section 4 in the 
vacuum limit p^— >0, we find the solution 

Nt = exp[ - -^—(^^PmC \ ■ (C.6) 

We see from (C.6) that Nt is constant in AT. Note that it is also possible to find isotropic 
null fluid models where pm depends on t, i.e., where there is local creation of null particles. 
For such models Nt will depend on t in JJ. Also note that more general cosmological 
models with homogeneous but non-isotropic matter distributions should exist in quasi- 
metric gravity. In particular models where the matter density on each FHS fluctuates 
around some mean density should be possible according to our theory since the FHSs 
cannot be isotropic in such cases. 

The result of this appendix confirms the assertion that in quasi-metric theory, the rate 
of the expansion of the Universe should not depend on dynamical gravitational parameters 
such as matter density. Rather the cosmic expansion must be purely non-kinematical and 
thus a consequence of the global curvature of space. Compared to traditional cosmology, 
this represents a radical difference in philosophy. 



D Symbols and acronyms 

In this appendix we give a table of symbols and acronyms used in the text. Sign con- 
ventions and notations follow reference [4] if otherwise not stated. In particular Greek 
indices may take values — 3 and Latin indices may take values 1 — 3. 
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SYMBOL NAME 



KE/NKE 

FOs/FHSs 

NKR/NKCE 

GTCS 

HOGS 

WEP/EEP 

GWEP/SEP 

LLI/LPI 

SR/GR 

SG 

t = t{x'') 

M/M 
AT/AT 
gt or g(t)f,u 
gt or g(^t)^Lu 
St or g^t)^,u 

ut or n|;) 

or nf,) 
hj or 
ht or 

pa 
pa 

pa 

; ('^'^^ 

pa pa 



ra pa 
t75 ■'-/37 



a* or a^j) 



kinematical/non-kinematical evolution (of the FHSs) 
fundamental observers/fundamental hypersurfaces 
non-kinematical redshift/non-kinematical coupling equation 
global time coordinate system 
hypersurface-orthogonal coordinate system 
weak/Einstein equivalence principle 
gravitational weak/strong equivalence principle 
local Lorentz/positional invariance 
special/general relativity 
Schiff's conjecture 

global time function, conventionally we define t = x^ /c when 
using a GTGS 

metric space-time manifolds obtained by holding t constant 
quasi-metric space-time manifolds (foliated by the FHSs) 
one-parameter family of space-time metric tensors on J\f 
one-parameter family of space-time metric tensors on AT 
one-parameter family of auxiliary space-time metric tensors 
on AT 

unit vector fields normal to the FHSs in the metrics 
unit vector fields normal to the FHSs in the metrics gj 
spatial metric tensors of the FHSs obtained from g^ 
spatial metric tensors of the FHSs obtained from gt 
lapse function in the metrics gt 
lapse function family in the metrics gt 
components of the shift vector field in the metric gtg 
components of the shift vector field in the metric gt^ 
components of the metric connections Vt compatible with gt 
components of the metric connections Vt compatible with gt 
components of the metric connections Vt compatible with gt 
components of the degenerate 5-dimensional connection V 
associated with the family gt 

components of the degenerate 5-dimensional connection V 

associated with the family gt 

family of metric 4-accelerations (any observer) 
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Q-a coordinate expression for metric covariant derivation com- 

patible with a single member of gt or as appropriate 

()*a coordinate expression for degenerate covariant derivation 

compatible with the family 

()|j coordinate expression for spatial covariant derivation com- 

patible with ht or ht as appropriate 

±/l projection symbols (projection on the normal direction to 

the FHSs) 

£x Lie derivative with respect to X in / 

£x Lie derivative with respect to X in Af/Af 

£x. Lie derivative of space objects with respect to X in A^/A^ 

* 

£x Lie derivative of space objects with respect to X in N /N 

ajp/ajr or a^j^/a^j: four-acceleration of the FOs in M./ M. 

hjr or }fj: 3-vector field family representing the local "generalized dis- 
tance vector to the center of gravity" 

efe/e'' or e\/e\ unit vector/covector along bjr 

Wj or w^j) the difference in coordinate 3-velocity between a test particle 

and the local FO 

Uf or -u^^-j 4-velocity of test particle 

or f 3-vector fields giving the motion of FOs in M compared to 

FOs in Jf 

a or a'' degenerate 4- acceleration in J\ (any observer) 

R(t)fiu/R(t)fiiy Ricci tensors in Al/ At 

H(t)ij/ H(t)ij spatial Einstein tensors of the FHSs m. M./M. 

P{t)ij / P{t)ij spatial Ricci tensors of the FHSs m. M./M. 

Pt/Pt spatial Ricci scalars of the FHSs in M/M 

or K{t)i3 family of extrinsic curvature tensors of the FHSs in M. 

K-t or ^(t)ij family of extrinsic curvature tensors of the FHSs in M. 

Ht/vt local measures of the NKE 

Xt local measure of the KE 

dTt/dTt proper time interval (any observer) in Af/JJ" 

drj^/drj: proper time interval for a FO in Af /JI 

Tt or stress-energy tensor as an active source of gravitation 

G/k Newton/ Einstein gravitational constant 

rrit, Mt active masses (scalar fields measured gravitationally) 

Ft=Ntct='^t^ local scale factors of the FHSs in Af 
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The Gravitational Field outside a Spherically Symmetric, 
Isolated Source in a Quasi-Metric Theory of Gravity 



by 

Dag 0stvang 

Institutt for Fysikk, Norges teknisk-naturvitenskapelige universitet, NTNU 

N-7491 Trondheim, Norway 

Abstract 

Working within the "quasi-metric" framework (QMF) described elsewhere [3], 
we find an approximate expression for an isotropic vacuum gravitational field in a 
S^xR-background and set up equations of motion applying to inertial test particles 
moving in this field. It is found that such a gravitational field is not static with 
respect to the cosmic expansion; i.e., distances between circular orbits increase ac- 
cording to the Hubble law. Furthermore it is found that the dynamically measured 
mass of the source increases with cosmic scale; this is a consequence of the fact that 
in the QMF, the cosmic expansion is not a kinematical phenomenon. The relevance 
of this to the observed asymptotically non-Keplerian behaviour of distant orbiters 
in spiral galaxies is discussed. Also it is shown that if this model of an expanding 
gravitational field is taken to represent the gravitational field of the solar system, 
this has no serious consequences for observational aspects of planetary motion. 

1 Introduction 

In metric theory the nature of the cosmological expansion is in principle not different 
from other types of motion; this follows from the most basic postulate in metric theory, 
namely that space-time can be modelled as a semi-Riemannian manifold. Furthermore, 
it is well-known that metric theory predicts that local systems are hardly affected at all 
by the cosmological expansion (its effect is at best totally neghgible, see e.g. [8] and 
references therein). However, when analyzing the influence of the cosmological expansion 
on local systems there should be no reason to expect that predictions made within the 
metric framework should continue to hold in a theory where the structure of space-time 
is non-metric. 

Recently a new type of non-metric space-time framework, the so-called "quasi-metric" 
framework, was presented in [3]. Also presented was an alternative relativistic theory of 
gravity formulated within this framework. This theory correctly predicts the "classic" 
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solar system tests in the case where an asymptotically Minkowski background is invoked 
as an approximation and the cosmological expansion can be neglected. However, for 
reasons explained in [3], the theory predicts a S^xR-background rather than a Minkowski 
background for the gravitational field outside a spherically symmetric, isolated source. 
According to the quasi-metric framework, the nature of the cosmological expansion is 
a direct consequence of the existence of this S^xR-background. Thus, to find how the 
cosmological expansion affects local systems according to the quasi-metric theory, one 
must first calculate the spherically symmetric gravitational field with this background 
and then calculate how test particles move in this field. 

Of particular interest is if the cosmological expansion has any effects on the corre- 
spondence with Newton's theory of gravity in the non-relativistic limit. If this is the case 
one may hope to explain from first principles the asymptotically non-Keplerian rotational 
curves of spiral galaxies as inferred from observations of Doppler shifted 21-cm spectral 
lines of neutral hydrogen. (See reference [1] for detailed information.) 

It has been shown that it is impossible to construct a realistic metric theory of gravity 
(subject to some rather general criteria of structure) such that the solar system tests come 
in and in addition galactic observations can be explained without dark matter [2]. But 
one remaining alternative to the introduction of dark matter may be the possibility of 
constructing a quasi-metric relativistic theory of gravity having the desirable properties. 
In this paper we explore if non-metric aspects of our theory correctly predict galactic 
observations from first principles, without introducing extra ad hoc parameters. 

2 General equations of motion 

We start by briefly summarizing the basics of our quasi-metric theory, see reference [3] 
for details. 

The geometric foundation of the quasi-metric framework consists of a 4-dimensional 
space-time manifold J\f equipped with a one-parameter family of Lorentzian 4- metrics gt 
parameterized by the global time function t. The global time function slices out a family 
St of preferred spatial submanifolds; we call these the fundamental hypersurfaces (FHSs). 
Observers always moving orthogonally to the FHSs are called fundamental observers 
(FOs). Each member of the family gt applies exactly on the fundamental hypersurface 
t ^constant only, but it can be extrapolated at least to a (infinitesimally) small normal 
interval around this hypersurface. (Other spatial hypersurfaces than the FHSs cannot 
be associated with the domains of applicability of single members of the metric family 
gt.) In component notation we formally write gt as g(t)ij,u, where the dependence on t is 
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put in a parenthesis not to confuse it with the coordinate indices. The components of 
each single metric are obtained by treating i as a constant. Since the theory is concerned 
primarily with the evolution of the geometry of the the FHSs along the world hues of 
the FOs we consider only a particular set of coordinate systems; namely the global time 
coordinate systems (GTCSs). By definition a coordinate system is a GTCS if and only 
if the relation ct — x° exists in J\f between the global time function and a global time 
coordinate x'^. 

The evolution of the scale factor of the FHSs may be thought of as consisting of 
a combination of two different contributions; we call these parts the non-kinematical 
evolution (NKE) and the kinematical evolution (KE), respectively. We define these terms 
in the following. 

A peculiarity of the quasi-metric framework is that the kinematical metric family gt 
in general does not represent a solution of field equations. The reason for this is that 
gt may contain some features inherently coming from the NKE. To avoid having these 
features as an explicit part of the dynamics we introduce a second family of metrics 
which by definition does represent a solution of field equations. The first family may 
then be constructed from the second family in a way described in [3]. 

In a GTCS the general form of the family g^ may be represented by the family of line 
elements 



ds] = 



{dx^f + 2-Nidx'dx° + %N^Sijdx'dx^, (1) 

^0 ^0 



where Nt is the lapse function field family of the FOs, ^N^^ is the family of shift 
vector fields of the FOs and where Sijdx^dx^ denotes the metric of the 3-sphere (with a 
radius equal to cto). Note that Ni=N^SikN''. 

The NKE of the FHSs is defined to take the form of an increase with time of the local 
scale factor Ft=Ntct of the FHSs as determined by the FOs using the metric family g^. 
The non-kinematical coupling equation (NKCE) of our theory yields a relation between 
the rate of increase of local scale and the global curvature of the FHSs in addition to a 
term yt coming from the gravitational field; this second term is not "realized" globally. 
That is, it can be argued [3] that the NKCE should read 

1 



where denotes the Lie derivative in the direction normal to the FHSs and h^^^ are 
the components of the inverse spatial metric family h^^. There is also a kinematical 
evolution (KE) of the Ft; this is described by the quantity c~^Xt. The overall evolution 
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of Ft is given by (note that a global term still occurs exphcitly here but that there 
is no term C^yt) 

c-%, + c-5,.-^ + (c*,)-(i + ^), (3) 

where the projection symbol X denotes a scalar product with —fif, which is a family 
of unit normal vector fields to the FHSs. From equation (1) we see that the spatial 
metric family is conformal to a "basic" 3-geometry, namely S^; this implies that the 
FHSs are compact (without boundaries). Furthermore we can find a set of field equations 
whereof gj must be a solution. The detailed form of the field equations can be found 
in [3]. Note the characteristic property of our theory that there are systematic scale 
changes between gravitational and atomic systems; this causes gravitational quantities 
to exhibit an extra formal variation when measured in atomic units (and vice versa). This 
may be thought of formally as if fixed operationally defined atomic units vary in space- 
time. Since c and Planck's constant h by definition are not formally variable the formal 
variation of time units is equal to that of length units and inverse to that of mass units. 
We represent the formal variation of atomic time units by the scalar field ^it=F^^. This 
means that measured in atomic units, the formal variability of gravitational quantities 
with the dimension of time or length goes as ^t^. By definition the gravitational constant 
G does not vary in space-time, but we then have to differentiate between active mass 
rrit as an active source of gravitation and passive mass m (passive gravitational mass or 
inertial mass). Active mass then varies formally as measured in atomic units. 

The formal evolutions of and in the normal direction to the FHSs are defined 

by [3] 

* 

Furthermore, on J\ there applies a linear, symmetric "degenerate" connection V . This 
connection is called degenerate due to the fact that the family may be perceived as 
one single degenerate 5-dimensional metric on a product manifold AlxR whereof H is 
a 4-dimensional submanifold. (Here Al is a Lorentzian space-time manifold and R is 
the real line.) We may then introduce a torsion-free, metric-compatible 5-dimensional 
connection V with the property that 

V^gt = 0, VAnt = 0, VAht = 0, (5) 

dt dt dt 

★ 

on xR and consider the restriction of V to N. It can be shown [3] that in a GTCS the 
components which do not vanish identically of the degenerate connection field are given 
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by (a comma denotes partial derivation, and we use Einstein's summation convention 
throughout) 

* 1 * 1 

The general equations of motion for test particles are identical to the geodesic equation 
obtained from V . In a GTCS they take the form (see [3] for a derivation) 

where drt is the proper time as measured along the curve, A is some general affine pa- 
rameter and at is the 4-acceleration as measured along the curve. 



2.1 Special equations of motion 

In this paper we analyze the equations of motion (7) in the case of a uniformly expanding 
isotropic gravitational field in vacuum exterior to an isolated spherically symmetric source 
in an isotropic, compact spatial background. We also require that the source is at rest 
with respect to some GTCS. It turns out that for an isolated system, also the FOs must 
be at rest with respect to this GTCS and consequently the shift vector field vanishes. 
Furthermore we require that Nf is independent of x^ and t\ i.e., that the only explicit time 
dependence is via t in the spatial scale factor (using the chosen GTCS). We call this a 
"metrically static" case. This scenario may be taken as a generalization of the analogous 
case with a Minkowski background (this case is analyzed in [3]) since the Minkowski 
background is not a part of our theory but rather invoked as an approximation being 
useful in particular cases. 

We start by making a specific ansatz for the form of gf Introducing a spherically 
symmetric GTCS {x°,r, ^,0}, where r is a Schwarzschild radial coordinate, we assume 
that the metric families gt and gt can be written in the form (compatible with equation 

(1)) 

cW' = 5(r)((ix°)2 - {-y(A{r)dr^ + r^dQ^), 

c^dx-BM(..«)^-(i)H^W*VrW). (8) 

where dQ'^=d9'^ + sm^ 9 dcf)'^ and to is some arbitrary reference epoch. Note that the spatial 
coordinate system covers only half of S^, thus the range of the radial coordinate is r < cto 
only. The functions A{r) and B{r) may be calculated from the field equations; we tackle 
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this problem in the next section. The functions A{r) and B{r) may be found from gj 
and yt as shown in [3]. 

We now calculate the metric and the degenerate connection coefficients from the 
metric family gt given in equation (8). A straightforward calculation yields (using the 
notation ' 



A'(r) r 

-ir- _ V ' pr _ ___ pr- _ pr o^n^/: 



r 



(t)oo - ^J) 2A(r)' ^ (t)re - ^ {t)er - i (t)^<^ - -smf^cosy, 



p</' _ p<^ _ 1 p<^ _ p0 _ „„4-/3 po — po — ^ ^^-^ fq^ 

^ (t)r</. - ^ it)<t>r - ^' it)<t>e - ^ it)e4> - ^O^*^' ^ (*)0r " {t)rO - 2B[r) ' ^ 



t = fte=f1, = ]- (10) 

In the following we use the equations of motion (7) to find the paths of inertial test 
particles moving in the metric family g^. Since vanishes for inertial test particles we 
get the relevant equations by inserting the expressions (9) and (10) into (7). This yields 
(making explicit use of the fact that cdt — dx^ in a GTCS) 



(Pr A'{r-) /dr\2 r 
^ 2A{r) ^dX^ ~A{^ 



I (io^2 B'{r) ^dx\2 1 drdx"" _ 

'^U' 2A{r)^d\) ^ ctdXdX ' ^ ' 



d^0 2d0dr ^ ./t^0\2 IdOdx'^ ^ 



c^V 2d6dr dS dO ld6dx^_ 

dX^ + rdXdX + '°^^dA + adXlX ~ ^ 
rf^xO B'{r) fix" dr 



d\^ B{r) dX dX 



0. (14) 



If we restrict the motion to the equatorial plane (12) becomes vacuous, and (13) reduces 
to 

d?(t) 2d(t)dr 1 d(t)dx^ _ 

dx^^rdXdx^adx~dx~ ■ ^ ^ 
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Dividing equation (15) by ^ we find (assuming ^7^0) 



d 

dX 



t 



0. 



We thus have a constant of the motion, namely 

t 2d(f) 
J=—r —. 
to dX 

Dividing equation (14) by ^ yields 



(16) 



(17) 



dX 



'dx\ 



ln(^)+lnS(r) 

^ OA ^ 



= 0. 



(18) 



Equation (18) yields a constant of the motion which we can absorb into the definition of 
A such that a solution of (18) is [4] 

d^ 
dX 



1 



Multiplying equation (11) by 



2^A{r)dr 



Bir) 



and using the expressions (17), (19) we find 



(19) 



d 


rtM(r) 




dX 


L Iq 


^dx) 



B{r) 



+ 



0, 



thus a constant E of the motion is defined by 

1 



t'^A{r) (dr\2 



/2 



/ at \ 



B{r) 



J2 



(20) 



(21) 



The equation (21) may be compared to an analogous expression obtained for the static 
isotropic gravitational field in the metric framework [4]. Inserting the formulae (17), (19) 
and (21) into (8) and using the fact that in a GTCS we can formally write dx^ — cdt 
when traversing the family of metrics, we find 



c'dr^ = EdX\ 



(22) 



Thus our equations of motion (7) force drt/dX to be constant, quite similarly to the case 
when the total connection is metric, as in the metric framework. Prom equation (22) we 
see that we must have £■ = for photons and £^ > for material particles. 

We may ehminate the parameter A from (17), (19), (21) and (22) and alternatively 
use t as a time parameter. This yields 



t 2d<P o/ N T 

— r —- — -D r J, 
to cdt ^ ^ 



(23) 
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B{r) 



1 




(24) 



dT^ = EB\r)dt\ 



(25) 



We may integrate equations (23) and (24) to find the time history (r(t),0(t)) along the 
curve if the functions A{r) and B{r) are known. 

For the static vacuum metric case with no NKE, one can solve the geodesic equation 
for particles orbiting in circles with different radii, and from this find the asymptotically 
Keplerian nature of the corresponding rotational curve [4]. In our case we sec from (23) 
and (24) that we can find circle orbits as solutions; such orbits have the property that 
the orbital speed B~^/'^{r)^r^ is independent of t. 

3 The dynamical problem 

In this section we write down equations which must be solved to find the unknown 
functions A{r) and B{r). (These equations follow directly from the field equations.) 
We then solve these equations approximately and use this to find approximations to the 
functions A{r) and B{r) as well. 

The field equations determining the metric family of the type shown in (8) may be 
set up exactly as for the corresponding case with a Minkowski background treated in [3] . 
The main difference between a Minkowski background case and the case of a compact 
spatial background is the explicit time dependence present for the latter case. However 
this time dependence is only via the scale factor. Now, from equations (1), (2) and (8) 



we get (since Nt 




B{r) is independent of t) 



dh{t)ij _ 2 
dt ~ t 




to cB'{r) 



(26) 




After some tedious but straightforward calculations the field equations yield 




(27) 




(28) 
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Before we try to solve equation (27), it is important to notice that no solution of it 
can exist on a whole FHS (except the trivial solution B ^constant), according to the 
maximum principle applied to a closed Riemannian 3-manifold. The reason for this 
is the particular form of the field equations, see reference [3] and references therein for 
justification. This means that in quasi-metric theory, isolated systems cannot exist except 
as an approximation. 

Even if a non-trivial solution of equation (27) does not exist on a whole FHS we 
may try to find a solution valid in some finite region of a FHS. That is, we want to find 
a solution in the region TZ < r-Ccto, where TZ is the coordinate radius of the source. 
We do not specify any particular boundary conditions at the boundary So=cto since 
the approximation made by assuming an isolated system is physically reliable only if 
^^1. Furthermore it turns out to be difficult to find a useful exact solution of equation 
(27). Thus we may rather try to find an approximative solution in the form of a series 
expansion, i.e., a perturbation around the analogous problem in a Minkowski background. 
That is, in contrast to the analogous case with a Minkowski background there exists the 
extra scale Sq in addition to the Schwarzschild radius r^o determined by the mass of the 
central object. (The extra index used in r^o refers to the time to and is necessary since 
in our theory gravitational quantities get an extra formal variability measured in atomic 
units.) Since we try to model the gravitational field exterior to galactic-sized objects we 
may assume that the typical scales involved are determined by this criterion 

tells how to compare the importance of the different terms of the series expansion. 

One may straightforwardly show that an approximative solution of equation (27) is 

B(r) = + ^ + +... ^ 

Mr) = + ^ + ^ + (29) 

and furthermore this yields, using equations (2) and (26) 



yt = T??[l + ? + 0(%> Ht = yt + ^^. (30) 



t 2r^ 2r ' r^ " - ^ B{r)t 

To construct the family gt as described in [3] we need the quantity {bj^ — r-^) 

2 

^(r)=y,V^(,),fe6^6^ = W^l^r = ^[1 + 0(^)]. (31) 
We note that v{r) does not depend on t. The functions A{r) and B{r) are found from 
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the following relations, valid for the spherically symmetric case [3] 



B{t). 



(32) 



Prom equations (29) and (32) we may then write 




(33) 



This expression represents the wanted metric family as a series expansion. Note in par- 
ticular the fact that all spatial dimensions expand whereas the corresponding Newtonian 
potential —U = —^-^ (to Newtonian order) remains constant for a fixed FO. This means 
that the true radius of any circle orbit (i.e., with r constant) increases but such that the 
orbit velocity remains constant. That is, the mass of the central object as measured by dis- 
tant arbiters increases to exactly balance the effect on circle orbit velocities of expanding 
circle radii. This is not as outrageous as it may seem due to the extra formal variation of 
atomic units built into our theory, rather it is a consequence of the fact that the coupling 
between matter and geometry depends directly on this formal variation [3] . 

What is measured by means of distant orbiters is not the "bare" mass Mt itself but 
rather the combination MtG. We have, however, defined G to be a constant. And as 
might be expected, it turns out that the variation of MtG with t as inferred from equation 
(33) is exactly that found directly from the formal variation of the active mass Mt with t 
by using equation (4). This means that the dynamically measured mass increase should 
not be taken as an indication of actual particle creation but that the general dynamically 
measured mass scale should be taken to change via a linear increase of Mt with t, and 
that this is directly reflected in the gravitational field of the source. That is, measured 
in atomic units, active mass increases linearly with epoch. 

The dynamical measurement of the mass of the central object by means of distant 
orbiters does not represent a local test experiment. Nevertheless the dynamically mea- 
sured mass increase thus found is just as "real" as the expansion in the sense that neither 
should be neglected on extended scales. This must be so since in quasi-metric relativity, 
the global scale increase and the dynamically measured mass increase are two different 
aspects of the same basic phenomenon. 

We finish this article by exploring which kinds of free-fall orbits we get from equation 
(33) and the equations of motion. We also discuss observable consequences thereof. 
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4 Discussion 



To begin with we find the shape of the rotational curve as defined by the orbital velocities 
of the circle orbits; since equation (24) has no time dependence for such orbits we can do 
a standard calculation [4] and the result is 



/ \ ri-1 /2 / \ ^ 

w[r)=B ' [r)—r— — 



^'^'K. (34) 



to dt \ 2B{r) 

However when we apply (34) to the metrics (33) we get a result essentially identical to 
the standard Kcplerian rotational curve; the only effect of the dynamically measured 
mass increase and the non-kinematical expansion is to increase the scale but such that 
the shape of the rotational curve is unaffected. It is true that B{r) as found from (33) 
contains a term linear in r in addition to terms falling off with increasing r; in reference [5] 
it is shown that such a linear term may be successfully used to model the asymptotically 
non-Keplerian rotational curves of spiral galaxies. Unfortunately the numerical value of 
the linear term found from (33) is too small by a factor of order 10^*^ to be able to match 
the data. So at least the simple model considered in this paper is unable to obtain the 
asymptotically non-Keplerian rotational curves of spiral galaxies from first principles. 

Another matter is how the time dependence in the equations of motion will affect 
the time histories and shapes of more general orbits than the circle orbits. Clearly time 
histories will be affected as can be seen directly from (24). However to see if this is vahd 
for shapes as well we may insert equation (23) into equation (24) to obtain r as a function 
of (f). This yields 

^r^^Vl-^ = -^_ (35) 



d(j)^ r2 .PB{r-) J2' 



4 



and this is identical to the equation valid for the case of a single spherically symmetric 
static metric [4]. Thus the shapes of free-fall orbits are unaffected by the non-kinematical 
expansion present in the metric family (33). 



4.1 Expanding space and the solar system 

If one neglects the gravitational effects of the galaxy one may try to apply the metric 
family (33) to the solar system (by using it to describe the gravitational field of the 
Sun). But the solar system is not at rest with respect to the cosmic rest frame; this 
follows from the observed dipole in the cosmic microwave background radiation. Besides, 
if we identify the cosmic rest frame with a GTCS, the FOs will be at rest with respect 
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to this GTCS only in the hmit where the motions of local gravitational sources can be 
neglected. Consequently, if we use a GTCS identified with the cosmic rest frame it is 
necessary to introduce a non- vanishing shift vector field family taking into account the 
motion of the solar system barycenter. And expressed in such a GTCS the metric family 
describing the solar system gravitational field cannot be metrically static since the metric 
coefficients will now depend on x^. However, as a good approximation we may neglect 
the gravitational effects of the galaxy and treat the solar system as an isolated system 
moving with constant velocity with respect to the cosmic rest frame. In this case we are 
always able to find a GTCS in which the solar system barycenter is at rest and where the 
metric coefficients do not depend on by construction (neglecting gravitational fields of 
bodies other than the Sun). Furthermore it turns out that the FOs must be at rest with 
respect to this new GTCS, so the shift vector field vanishes and all calculations we have 
made in this paper are still valid. Thus, if we treat the solar system as approximately 
isolated we can neglect its motion with respect to the cosmic rest frame and use the metric 
family (33) to describe the gravitational field of the Sun. Moreover, the solar system is 
so small that we can neglect any dependence on Sq. The errors made by neglecting terms 
depending on Sq in equation (33) are insignificant since the typical scales involved for 
the solar system are determined by ^ ^ Equation (33) then takes the form 

ds', = -(l - ^ + oC-^)){dx'f + (^)2({1 + ^ + 0C-^)}dr' + r'dn'). (36) 

Equation (35) shows that the shapes of orbits are unaffected by the expansion; this 
means that all the classical solar system tests come in just as for the analogous case of 
a Minkowski background [3]. However, we get at least one extra prediction (irrespective 
of whether the galactic gravitational field can be neglected or not); from equation (36) 
we see that the effective distance between the sun and any planet has been smaller in 
the past. That is, the spatial coordinates are comoving rather than static, thus the solar 
system is not static (measured in atomic units). For example the distance between the 
Sun and the Earth at the time of its formation may have been about 50% less than 
today. But since gravitationally bound bodies made of an ideal gas are predicted to 
expand according to quasi- metric theory [10], a small Earth-Sun distance should not be 
incompatible with palaeo-climatic data since the Sun is expected to have been smaller 
and thus dimmer in the past. Actually, since neither the temperature at the center of 
the Sun (as estimated from the virial theorem) nor the radiation energy gradient times 
the mean free path of a photon depends on t, the cosmic luminosity evolution of the 
Sun should be determined from the cosmic expansion of its surface area as long as the 
ideal gas approximation is sufficient. And this luminosity evolution exactly balances the 
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effects of an increasing Earth-Sun distance on the effective solar radiation received at the 
Earth. 

However, an obvious question is if the predicted effect of the expansion on the time 
histories of non-relativistic orbits is compatible with the observed motions of the planets. 
In order to try to answer this question it is illustrating to calculate how the orbit period 
of any planet depends on t. For simplicity consider a circular orbit r — R ^constant. 
Equation (23) then yields 

f - ^(«'«-^^'^ P^) 

Now integrate (37) one orbit period T « t (i.e., from t to t + T). The result is 

T(*,^.(exp[5-]-l)4wi.5M^...,, Ta.^^, (38) 



where Tgr is the orbit period as predicted from GR. Prom (38) we see that (sidereal) 
orbit periods are predicted to increase linearly with cosmic scale, i.e. 

^/ X t , dT T(to) , , 

m = :rT{to), — = (39) 

and such that any ratio between periods of different orbits remains constant. In particular 
(39) predicts that the (sidereal) year Te should be increasing with about 2.5 ms per year 
and the martian (sidereal) year Tm should be increasing with about 4.7 ms per martian 
year at the present epoch. This should be consistent with observations since the observed 
difference in the synodical periods of Mars and the Earth is accurate to about 5 ms. 

To compare predictions coming from equation (36) against timekeeping data one must 
also take into account the cosmological contribution to the spin-down of the Earth. If 
one assumes that the gravitational source of the exterior field (36) is stable with respect 
to internal collapse (as for a source made of ideal gas), i.e., that possible instabilities 
generated by the expansion can be neglected, one may model this source as a uniformly 
expanding sphere [10]. Due to the cosmic expansion of gravitational fields, the angular 
momenta of test particles moving in the exterior field (36) increase linearly with cosmic 
scale, and this should also apply to the angular momentum Lg of a passive mass source, 
that is 

L,{t) = -Lsih), — i = -L, = (1 + 0{2))HL,, (40) 
to at z 

where the term 0{2) is of post-Newtonian order and where the locally measured Hubble 
parameter H is defined by H=^, or equivalently (tjf is the proper time of the local FO) 



'') = ^(y5w)-;^(f) = (V5M*)"'. (41) 



t drjr to t ' dx^ to 
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Since the moment of inertia IcxMR'^, where M is the passive mass and Rg is the measured 
radius of the sphere, we must have (neglecting terms of post-Newtonian order) 

dt ' dt ' dt ' ^ ^ 

where Us is the spin circle frequency and T, is the spin period of the sphere. (To show (42), 

use the definition Ls=IlJs-) This means that the spin period of a sphere made of ideal 

gas increases linearly with t due to the cosmic expansion. Does this apply to the Earth as 

well? The Earth is not made of ideal gas, so the cosmic expansion may induce instabilities, 

affecting its spin period Tse- However, here we assume that the Earth's mantle is made of 

a material which may be approximately modelled as a perfect fluid obeying an equation 

of state close to linear. Then, if this assumption holds, the effects of instabilities averaged 

over long time spans should be negligible to a good approximation. We may also assume 

that there is no significant tidal friction since given the cosmic contribution, this would 

be inconsistent with the observed mean acceleration of the Moon hm (see equation (46) 

below). We then get 

^^^fM^HT.Eito), nE{t)^^TMto). (43) 

dt Iq 

Prom equation (43) we may estimate a cosmic spin-down of the Earth at the present epoch 

of about 0.68 ms per century (using if~2.5xl0~^^ s"^). To see if this is consistent with 

the assumption that the dominant contribution is due to cosmic effects we may compare 

to historical observations indicating a lengthening of the day of about 1.4 ms per century 

averaged over the last 1000 years [12] and of about 1.7 ms per century averaged over 

the last 2700 years [13]. But the historical observations depend on an assumed value for 

fim obtained from fitting lunar laser ranging (LLR) data to a Newtonian model, making 

these values theory dependent. This theory dependence also affects the comparison of 

equation (43) with palaeo-tidal records obtained from sedimentary tidal rhythmities [11]. 

However, if the quasi-metric prediction of rim (see equation (46) below) rather than the 

LLR value is used as an input to the historical observations, one gets a result close to 

the quasi- metric prediction of the spin-down of the Earth obtained from equation (43). 

Another quantity that can be calculated from equation (43) is the number of days Ny 

in one (sidereal) year Te- This is found to be constant since 

TE^NyT,E^lTE{to), ^ ^ = 0- (44) 

A similar calculation applies to the number of the days in one (sidereal) month T^, 
i.e. 

Tm^NmTsE^r^mih), ^ = 0. (45) 

to dt 
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That Ny and Nm are predicted to be constant is not in agreement with standard (theory 
dependent) interpretations of palaeo-geological data [7], [11]. (The predicted constancy 
of the ratio Ny/N^ agrees well with a standard interpretation of the data, though.) To 
see why, the assumption that Tg is constant is routinely used in the interpretation of 
tidal rhythmities and fossil coral growth data; in particular this applies to [11]. Another 
assumption used in the data interpretation is that active masses do not vary with time. 

The apparent constancy of the sidereal year (as indicated by astronomical observa- 
tions of the Sun and Mercury since about 1680 [9]) represents the observational basis for 
adopting the notion that ephemeris time (i.e., the time scale obtained from the observed 
motions of the Sun and the planets) is equal to atomic time (plus a conventional con- 
stant), but different from so-called universal time (any time scale based on the rotation of 
the Earth). However, from equations (39) and (43) we see that according to quasi- metric 
relativity, ephemeris time should be scaled with a factor ^ compared to atomic time, 
and averaged over long time spans it should be equal to universal time. But within the 
Newtonian framework, any secular changes in the Earth-Moon system are interpreted 
in terms of tidal friction (and external perturbations), so seemingly secular inconsisten- 
cies between different time scales may be blamed on the variable rotation of the Earth. 
In practice this means introducing leap seconds. Given the fact that leap seconds are 
routinely used to adjust the length of the year, the predicted relationship between the 
different time scales should be consistent with observations. In particular, the extra time 
corresponding to an increasing year as predicted from our model may easily be hidden 
into the declining number of days in a year as predicted from standard theory. 

As mentioned above, the value of the so-called mean acceleration of the Moon, 
njn=-^nm, (where rim is the mean geocentric angular velocity of the Moon as observed 
from its motion) is a very important quantity for calculating the evolution of the Earth- 
Moon system. Prom equation (45) we see that quasi-metric theory predicts that 

1T'm{t) = = V^m(^o), =^ nm=4:'>T'm = -HUm, (46) 

at t at 
and inserting the observed value 0.549"/s for rim at the present epoch, we get the corre- 
sponding cosmological contribution to hm, namely about —13.6" per century^. This value 
may be compared to the value —13.68" per century^ obtained from fitting LLR data to 
a lunar model [14]. Note that this second value is the total mean acceleration, wherein 
other modelled (positive) contributions, such as planetary perturbations, are included. 
When the other contributions are removed one deduces a tida/ contribution fitid — —25.8" 
per century^ to hm [14]. Similar results for the tidal contribution to hm inferred from 
LLR data have been found in, e.g., [15] {hud — —25.9" per century^). We see that the 
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quasi-metric result is smaller (in absolute value) than the tidal term. But as long as 
there are no independent measurements of the non-tidal contributions to fim, a model 
containing only the cosmic contribution fits the data well. 

We may also use Bubble's law directly to calculate the secular recession CLqmr of the 
Moon due to the global cosmic expansion; this yields about 3.0 cm per year whereas the 
value htid inferred from LLR is about 3.8 cm per year. However, LLR values are found 
by fitting the LLR data to a Newtonian model where active masses are assumed to be 
constant. Prom Kepler's third law we get the relationship 

_ 2 htid . 

where hqmr is the cosmic contribution. By inserting the numerical values given above we 
see that this is quite consistent with the data. 

To round off this paper: we have seen that the predicted effects of the cosmic expansion 
on the solar system gravitational field have a number of observable consequences, none 
of which is shown to be in conflict with observation so far. That is, it seems that at 
this time no model independent evidence exists that may rule out the possibility that 
planetary orbits expand according to the Hubble law. 
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Abstract 

According to a particular "quasi-metric" theory of gravity developed elsewhere 
[1], the apparently anomalous force acting on the Pioneer 10/11, Galileo and Ulysses 
spacecraft as inferred from radiometric data [2], [9], may be naturally explained as 
resulting from treating comoving coordinates as static ones. More precisely, we 
show that the illusion of an anomalous acceleration of the right size and acting to- 
wards the observer arises as a consequence of the mismodelling of null paths. The 
apparently anomalous acceleration is of order cH [2], [9] (where H is the Hubble 
parameter) as predicted by a simple non-static model of the solar system gravita- 
tional field [3]. 

1 Introduction 

Some time ago a detailed analysis of the observed versus the calculated orbits of the 
Pioneer 10, Pioneer 11, Ulysses and Galileo spacecraft was published [2]. The main 
result was that the observed radiometric data did not agree with calculations based on 
standard theory; rather the data indicated the existence of an "anomalous", constant 
acceleration towards the Sun. 

A short summary of the results presented in [2] is as follows: For the Pioneer space- 
craft the Doppler frequency shift of the radio carrier wave was recorded and analyzed 
to determine the spacecraft's orbits. Two independent analyses of the raw data were 
performed. Both showed an anomalous acceleration towards the Sun, of respectively 
(8.09±0.20)xl0-^ cm/s^ and (8.65±0.03) xlO"^ cm/s^ for Pioneer 10. For Pioneer 11 
only one result is given; an anomalous acceleration of (8.56±0.15) xlO~^ cm/s^ towards 
the Sun. The acceleration did not vary between 40 — 60 astronomical units, within a 
sensitivity of 2xl0~^ cm/s^. 

For the Gahleo and Ulysses spacecraft one also got ranging data in addition to the 
Doppler data. (Ranging data are generated by cross correlating a phase modulated signal 
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with a ground duplicate and measuring the time delay.) For Ulysses one had to model the 
solar radiation pressure in addition to any constant anomalous acceleration. By doing this 
it was found that Ulysses was influenced by an anomalous acceleration of (12±3)xl0~^ 
cm/s^ towards the Sun, consistent with both Doppler and ranging data. For Gahleo the 
corresponding result was an anomalous acceleration of (8±3)xl0~^ cm/s^ towards the 
Sun. 

Very recently a new comprehensive study of the anomalous acceleration was pub- 
lished [9], including a total error budget for the Pioneer 10 data analysis. The new 
result reported in [9] is an "experimental" anomalous acceleration of (7.84±0.01) xlO~^ 
cm/s^ towards the Sun, and including bias and uncertainty terms the final value becomes 
(8.74±0.94)xl0"^ cm/s^. For Pioneer 11 an experimental value of (8.55±0.02) xlO"^ 
cm/s^ was given. Also other new results, such as annular and diurnal variations in the 
anomalous acceleration, were reported in [9]. 

An interpretation of these results according to the standard general relativistic model 
indicates the existence of an anomalous, time- independent force acting on the spacecraft. 
However there are problems with this interpretation since according to the planetary 
ephemeris there is no indication that such a force acts on the orbits of the planets; the 
hypothetical force thus cannot be of gravitational origin without violating the weak prin- 
ciple of equivalence. Thus it is speculated that the effect is due to anisotropic radiation 
of waste heat from the radioactive thermal generators aboard the spacecraft; the design 
of the spacecraft is such that waste heat may possibly be scattered off the back of the 
high gain antennae in directions preferentially away from the Sun [4]. Moreover, besides 
possible anisotropic scattering, an estimate shows that the specific arrangement of waste 
heat radiators on the surface of the spacecraft may perhaps cause sufficient anisotropy in 
the radiative cooling to explain the data [5]. However, it seems that these explanations 
have been effectively refuted [6], [7], [8], [9]. Other possible explanations, such as gas 
leaks, have been proposed [6], [9], but so far it seems that no satisfactory explanation 
based on well-known physics exists. 

However, it is an intriguing fact that the size of the anomalous acceleration is of the 
order cH for all the spacecraft, where H is the Hubble parameter. Since this seems to be 
too much of a coincidence one may suspect that the data indicate the existence of new 
physics rather than a prosaic explanation based on standard theory. This has been duly 
noted by others, see [9] and references therein. But to be acceptable, any non-standard 
explanation should follow naturally from a general theoretical framework. In this paper 
we show that such an explanation can be found, thus the data may indeed be taken as 
evidence for new physics. 
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2 A quasi-metric model 



In reference [1] we defined the so-called "quasi-metric" space-time framework; this frame- 
work is non-metric since it is not based on semi-Riemannian geometry. Moreover, in refer- 
ence [3] we introduced a model of the gravitational field outside a spherically symmetric, 
isolated source as predicted by a particular quasi-metric theory of gravity developed in 
[1]. According to this theory it was found in [3] that at scales similar to the size of the 
solar system, such a gravitational field can be expressed by the one-parameter family 
of Lorentzian 4-metrics 



where r is a comoving radial coordinate and dQ'^=d9'^ + sm'^9d(f)'^ is the squared solid angle 
line element. Furthermore t is the global time function; its value is zero at the beginning 
of the Universe. In equation (1) t is to be viewed as a parameter; the relationship 
between t and the space-time coordinates is given by x° = ct and to represents some 
arbitrary reference epoch setting the scale of the spatial coordinates. The reason that 
one must separate between ct and in (1) is that the afftne connection compatible with 
the family is non-metric. That is, although ct and both can be interpreted as time 
coordinates, the components of the affinc connection containing ct is not equivalent to 
their counterparts containing x^. See [1] for a further discussion. The locally measured 
Hubble parameter H as calculated from (1) reads [3] 



General equations of motion are obtained from the geodesic equation using the non-metric 
connection [1]. But these equations of motion cannot be obtained from the geodesic 
equation using any metric connection. Moreover, as shown in [3] special equations of 
motion for inertial test particles moving in the particular metric family (1) take the 
form (due to the spherical symmetry we can restrict the motion to the equatorial plane 




(1) 





9 = n/2) 




(3) 



t r,d(b 

— r — 
^0 cat 



B{r)J, 



(4) 



dr^ = -c-^ds'^ = EB^{r)dt'^, 



(5) 
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where J and E are constants of the motion. By setting the scale factor ^ = 1 in equations 
(1), (3), (4) and (5) we recover the equations of motion for inertial test particles moving 
in a spherically symmetric, static gravitational field as obtained from General Relativity 
(GR). Note that £■ = for photons and £■ > for material particles, which may readily 
be seen from equation (5). 

The functions A{r) and B{r) may be found as series expansions by solving the field 
equations, this is done approximately in [3]. For our purposes we include terms to post- 
Newtonian order but not higher. Then we have 

B(r) = l-'y + oCf), (6) 

where rso is the Schwarzschild radius at the arbitrary epoch to- 

We now explore some of the differences between the non-static system described by 
equation (1) and the corresponding static system obtained by setting = 1 in (1) and 
using GR. To begin with we notice that the shapes of free fall orbits (expressed e.g. as 
functions of the type r(0)) are identical for the two cases [3]. Moreover, it can be shown 
that the time dependence present in equations (3), (4) and (5) docs not lead to easily 
observable perturbations in the paths of non-relativistic particles compared to the static 
case [3]. However, as we now illustrate, if one considers null paths potential observable 
consequences appear if one treats r as a static coordinate rather than as a comoving 
one. To simplify matters we consider purely radial motion, i.e., J = (one may easily 
generalize to J^O). Since = for photons we get from equation (3) that radial null 
curves are described by the equation 



dr cto 




+ 0(^) = ±i^(l-^ + 0(^)). (7) 



dt t \^ A{r) t \ r t ^ r r 



the choice of sign depending on whether the motion is outwards or inwards. By integrating 
(7) to lowest order we find an extra delay, as compared to standard theory, in the time it 
takes an electromagnetic signal to travel from an object being observed to the observer. 
To lowest order this extra time delay is where R is the radial coordinate distance 
between the object and the observer and H is the Hubble parameter as given from 
equation (2). Also, the fact that the scale factor in equation (1) increases with time 
implies that our model predicts an extra redshift, as compared to standard static models, 
in the Doppler data obtained from any object emitting electromagnetic signals. To lowest 
order this extra redshift corresponds to a "Hubble" redshift 



124 



But the velocity at any given time of an observed object cannot model-independently 
be split up into one "ordinary" piece and one "Hubble" piece. This means that there is 
no direct way to identify the predicted extra redshift in the Doppler data. Similarly, at 
any given time there is no direct way to sort out the predicted extra time delay when 
determining the distance to the object. Rather, to test whether the gravitational field is 
static or not one should do observations over time and compare the observed motion to a 
model. In a model one uses a coordinate system and to calculate coordinate motion one 
needs coordinate accelerations. Accordingly we construct the "properly scaled coordinate 
acceleration" quantity Uc- For photons this is 

„ _ t \[Mr) (fr _ c toTso^ rj^ 
"-to B{r) dt^ ~ ^t^ t r2 ^"^^ r3 ^ 

^^cH + Y-^ + oC-^). (8) 



The point with this is to show that by treating the comoving coordinate system as a 
static one and using GR, an "anomalous" term ^cH will be missed when modelling 
coordinate accelerations of photons. We see that the sign of the anomalous term is such 
that the anomalous acceleration is oriented in the opposite direction to that of the motion 
of the photons. This means that to sufficient accuracy, treating the comoving coordinate 
system as a static one is equivalent to introducing a variable "effective" velocity of light 
Cg// equal to 

ceff = c(l - —)=c{l-HT + 0{{HTf)), (9) 

where T is the light time along the null path. The change of Cg// with T then yields an 
anomalous acceleration 

a„ = ^Ceff = -cH + 0{cH'T), (10) 

along the line of sight of any observed object. That is, if the comoving coordinate system 
is treated as a static one the coordinate motion of any object will be observed to slow 
down by an extra amount if the light time, or equivalently, the distance to the observer 
increases and to speed up by an extra amount if the distance decreases. Hence, judging 
from its coordinate motion it would seem as if the object were influenced by an anomalous 
force directed towards the observer. 

By integrating the anomalous acceleration over the total observation time T-Ct we 
get an "anomalous" speed 

rt+T 

Wa^j ttadt' ^-cHT + 0{cH'^T'^), (11) 
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towards the observer compared to a model where the coordinates are static rather than 
comoving. That is, the coordinate motion of any object observed over time indicates an 
anomalous blueshift compared to a model where the gravitational field is static. 

Now, since any observer is typically located at the Earth, the model in fact predicts 
an anomalous acceleration directed towards the Earth rather than towards the Sun. But 
any directional differences will almost average out over time if the observed object moves 
approximately radially and is located well beyond the Earth's orbit. However, compared 
to a model where the anomalous acceleration is towards the Sun, even if directional dif- 
ferences nearly average out there remains a cumulative difference. We will calculate this 
below. If, on the other hand, the line of sight to the observed object (e.g. a planet) devi- 
ates significantly from the radial direction, observations should not be consistent with an 
anomalous acceleration directed towards the Sun. Rather the direction of the anomalous 
acceleration expressed in Sun-centered coordinates would appear to be a complicated 
function of time. 

To estimate the predicted differences between a model where the anomalous acceler- 
ation is towards the Sun and the result given by equation (10), for the case where the 
observed object moves approximately radially in the ecliptic plane and is located well 
beyond the Earth's orbit it is convenient to define the average anomalous acceleration 
('^c)j„ of a photon away from the Sun during the time of flight T from the object to the 
observer. Thus we define 

1 /■*+^ c r ft+T.dt' _rT^^Ti 

It+Ti ■ 

where Ti is the moment when the photon crosses a plane through the center of the Sun 
normal to a line connecting the object and the Sun. (If the Earth is at the same side of 
this plane as the object, Ti = T.) Similarly we can define the average {ac)out ~ ~^'^c)in 
of a photon towards the Sun during the time of fiight from the observer to the object. 
We are now able to estimate the predicted difference 5a between a model where the 
anomalous acceleration is towards the Earth and one where it is towards the Sun. We 
get 

5a = -cH - (ae)„„, = -2cH{l - ^ + 0{HT)). (13) 

This function has a minimum at solar conjunction and vanishes when the observed object 
and the Earth are at the same side of the Sun. One may show that 5a can be written 
approximately as a truncated sine function where all positive values are replaced by zero. 
The period is equal to one year and the amplitude is approximately 2cif^, where Ro 
and Re are the radial coordinates of the object and of the Earth, respectively. 
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1 /■*+^ C r /"i+^n dt' / Ti \ 



Anderson et al. [9] (see also [6]) found an annual perturbation on top of the anoma- 
lous acceleration op of Pioneer 10. (They claim to see such a perturbation for Pioneer 11 
also.) Interestingly, the perturbation consistently showed minima near solar conjunction 
[9] (i.e., the absolute value of the anomalous acceleration reached maximum near solar 
conjunction). They fitted an annual sine wave to the velocity residuals coming from the 
annular perturbation term, using data from Pioneer 10 when the spacecraft was about 
60 AU from the Sun. Taking the derivative with respect to time they then found the am- 
plitude of the corresponding acceleration; it was found to be aa.t. = (0.215±0.022)xl0~^ 
cm/s^. We may compare this to the corresponding amplitude in Sa calculated above. We 
find an amplitude of about 0.29x10^^ cm/s^, close enough to be roughly consistent with 
the data. Besides the annual term a diurnal term was also found in the velocity residuals 
[9], where the corresponding acceleration amplitude aa.t. is large compared to ap. But 
over one year the contribution from the diurnal term averages out to insignificance (less 
than 0.03x10"^ cm/s^ [9]). 

2 2 

To find the trajectories of non-relativistic particles we may set E=l — where ^ 
is small. Then equation (3) yields 



dr _|_ cto 
'di~ 



(14) 



and the properly scaled coordinate acceleration for non-relativistic particles is 



We see that for non-relativistic particles the effect on coordinate accelerations of treating 
the comoving coordinates as static ones and using GR, is a factor + ^ smaller than 
the corresponding effect for photons. This means that the trajectories of non-relativistic 
particles do not depend crucially on the fact that the gravitational field is non-static. 

On the other hand the paths of photons depend more significantly on whether the 
gravitational field is static or not and this yields the illusion of an anomalous acceleration. 
That is, if one receives electromagnetic signals from some freely falling object located e.g. 
in the outer parts of the solar system, the coordinate acceleration of the object as inferred 
from the signals should not agree with the "real" coordinate acceleration of the object 
if one treats the comoving coordinates as static ones. Rather, from equation (10) we see 
that it would seem as if the object were infiuenced by an attractive anomalous acceleration 
of size cH. The relevance of this is apparent when modehng the orbits of spacecraft and 
comparing to data obtained from radio signals received from the spacecraft; in particular 
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this applies to the analyses performed in [2] , [6] and [9] . An extra bonus for the model 
considered in this paper is that it predicts small deviations during the year if the data are 
compared to a model where the anomalous acceleration is directed towards the Sun rather 
than towards the Earth. And as we have seen, this prediction seems to be consistent with 
the data. 

2.1 Cosmic expansion and the PPN- formalism 

Orbit analysis of objects moving in the solar system must be based on some assumptions 
of the nature of space-time postulated to hold there. The standard framework used for 
this purpose is the parameterized post-Newtonian (PPN) formalism applicable for most 
metric theories of gravity. But the standard PPN-framework does not contain any terms 
representing expanding space via a global scale factor as shown in equation (1) since it is 
inherently assumed that the solar system is decoupled from the cosmic expansion. One 
may try to overcome this by inventing some other sense of "expanding space" where the 
scale factor varies in space rather than in time. But such a model must necessarily be 
different from our quasi-metric model, and we show below that it cannot work. Thus, 
to illustrate the inadequacy of the PPN-framework to model expanding space we now 
consider a specific model where suitable terms are added by hand in the metric. We may 
then compare to the change in light time obtained from our quasi-metric model. 
One may try a post-Newtonian metric of the type 

= -(1 - 7 + ^ + 0(r|))(d.»)^ + (1 + l£ - + oCl))dr^ + r^6a\ (16) 

where Hq is a constant, to describe expanding space within the PPN-framework. (To show 
that this metric yields a spatially variable scale factor, transform to isotropic coordinates.) 
It may be readily shown that the metric (16) yields a constant anomalous acceleration cHq 
towards the origin. But the problem with all metrics of this type is that they represent a 
"real" anomalous acceleration of gravitational origin, and this is observationally excluded 
from observations of planetary orbits [9] . 

Anyway we may calculate the change in light time AT due to the terms containing 
Hq in (16) by integrating a radial null path from to (let Vp > re, say). This yields 

AT = -c-^Ho{rl - rl) + ■ ■ - HoT\ (17) 

where the light time T is equal to c~^{rp — Tg) to first order and where the last approx- 
imation is accurate only if Vp is much larger than rg. Note that AT is negative; this is 
quite counterintuitive for a model representing expanding space. 
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Anderson et al. [9] have considered a phenomenological model representing "expand- 
ing space" by adding a quadratic in time term to the hght time in order to determine 
the coefficient of the quadratic by comparing to data. To sufficient accuracy this model 
may be represented by the transformation 

T^{1 + aguadt)T=T + AT, (18) 

where ag^ad is a "time acceleration" term. This model fits both Doppler and range very 
well [9]. 

If we compare equations (17) and (18) we see that the value ttguad = ~Hqj corresponds 
to the change in light time calculated from the metric (16). This is far too small (of order 
10"^'' s~^ for a light time of a few hours) to be found directly from the tracking data. 
Thus, the fact that aquad was estimated to be zero based on the tracking data alone [9] 
in no way favours a constant acceleration model over a time acceleration model. They 
are in fact equivalent as far as the data are concerned. 

To compare differences in ttguad we may integrate equation (7) to lowest order to find 
the extra delay -^{rp — rg)^ = in the light time compared to the static case. This 
corresponds to a value ag^ad = for the time acceleration term. But a determination 
of o-quad dircctly from the tracking data still reflects the model-dependency explicitly 
present in the orbit determination process. This means that a determination of ag^ad 
directly from the tracking data should in principle be consistent with the model (16) and 
not with our quasi-metric model. However, the quantity aguad is so small that it is not 
feasible to check this. But the fact that a phenomenological model of the type (18) works 
so well should be taken to mean that the explanation of the anomalous acceleration given 
in this paper is sufficient. 

3 Conclusion 

We conclude that a natural explanation of the data is that the gravitational field of the 
solar system is not static with respect to the cosmic expansion. This also explains why 
any orbit analysis program based on the PPN-formalism is insufficient for the task and 
how the largest errors arise from the mismodelling of null paths. (In fact, using the PPN- 
formalism is equivalent to introducing a variable "effective" velocity of light as shown 
in equation (9).) But these explanations, while not involving any ad hoc assumptions, 
are based on the premise that space-time is quasi-metric. That is, rather than being 
described by one single Lorentzian metric, the gravitational field of the solar system 
should be modelled (to a first approximation) by the metric family shown in equation 
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(1). Prom a theoretical point of view this premise is radical; thus it is essential that the 
subject is further investigated to make certain that more mundane explanations may be 
eliminated. However, so far no such explanations based on well-known physics have been 
found. But the facts that the model presented in this paper follows from first principles 
and fits the data very well, together with independent observational evidence in favour 
of the prediction that the Earth-Moon system is not static with respect to the cosmic 
expansion [3] , indicate that explanations based on quasi-metric relativity should be taken 
seriously. 
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Abstract 

It is found that one of the predictions of a so-called "quasi-metric" theory of 
gravity developed elsewhere [1], is that the gravitational field inside a metrically 
static, spherically symmetric, isolated body modelled as a perfect fluid obeying a 
linear equation of state (i.e., pocpm), should expand according to the Hubble law. 
That is, for such a body its radius should increase like the expansion of the Universe; 
this is the counterpart to the similar result found for the corresponding exterior 
gravitational field [2]. On the other hand, if the body consists of a perfect fluid 
obeying some other equation of state (i.e., pgtpm, as for a polytrope for example), 
the expansion will induce instabilities and the body cannot be metrically static; 
such a body will not in general expand. 

1 Introduction 

A new type of non-metric space-time framework, the so-called "quasi-metric" framework, 
was presented in [1] . Also presented was a possibly viable quasi-metric theory of gravity 
compatible with the quasi-metric framework. This theory predicts that the nature of the 
Hubble expansion is different than in metric theory, and as a consequence the Hubble 
expansion should influence gravitationally bound systems, according to the new theory. In 
particular the theory predicts that the gravitational field outside a spherically symmetric, 
isolated source expands according to the Hubble law [2], and this seems to be consistent 
with observations of the orbits of distant spacecraft [3] . 

Now a natural question is what happens to the gravitational fields inside gravitation- 
ally bound bodies; does the quasi-metric theory predict that such gravitational fields may 
expand, too? We show in this paper that this is indeed the case, and that this result 
may support an interpretation of geological data indicating that the Earth is expanding, 
see reference [4] and references cited therein. Moreover an expanding Earth should cause 
changes in its spin rate; we showed in [2] that the main part of the observed secular 
spin-down of the Earth may in fact be of cosmological origin. This also applies to the 
observed recession of the Moon and its observed mean acceleration. 
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2 A brief survey of the quasi-metric theory 



In this section we summarize the main features of the quasi-metric framework and the 
quasi- metric theory. A considerably more extensive discussion can be found in [1]. 

The basic premise of the quasi-metric framework is that the canonical description of 
space-time is taken as fundamental. That is, space-time is constructed as consisting of two 
mutually orthogonal foliations: on the one hand space-time can be sliced up globally into 
a family of 3-dimcnsional spacclikc hypcrsurfaccs (called the fundamental hypcrsurfaccs 
(FHSs)) by the global time function t, on the other hand space-time can be foliated into 
a family of timelike curves everywhere orthogonal to the FHSs. These curves represent 
the world lines of a family of hypothetical observers called the fundamental observers 
(FOs), and the FHSs taken at t =constant represent a preferred notion of space. The 
equations of the theory depend on the geometrical properties of the FHSs. 

The FHSs are associated with two famihes of Lorentzian space-time metric tensors 
gt and gt in such a way that different FHSs correspond to domains of apphcabihty of 
different members of these families in an one to one relationship. The metric family g^ 
represents a solution of field equations, and from gt one can construct the "physical" 
metric family g^ which is used when comparing predictions to experiments. The theory 
is not metric since the affine connection compatible with any metric family is non-metric. 

To be able to compare theory to experiment we have to represent the metric families 
in terms of components with respect to some coordinate system on space-time. It is 
convenient to use a coordinate system {x^} where the relationship between the time 
coordinate and the global time function is given by x° = ct; this ensures that x^ is a 
global time coordinate. A coordinate system with a global time coordinate of this type 
we call a global time coordinate system (GTCS). Note that there exist infinitely many 
GTCSs. 

It can be argued [1] that in a GTCS the most general form allowed for the family g^ 
is represented by the family of line elements (this may be taken as a definition) 

f - f2 _ 

{dxy + 2-Nidx'dx^ + -^N^Sikdx'dx^, (1) 

^0 ^0 



ds] 



where Nf is the lapse function field family of the FOs and where j^N^dx'' is the family of 
shift covector fields of the FOs in g^. Moreover Sikdx^dx^ is the metric of the 3-sphere 
(with radius equal to cto). 

The evolution of the spatial scale factor of the FHSs in the hypersurface-orthogonal 
direction may be split up according to the definition (where a comma denotes partial 
derivation and where the symbol i denotes a scalar product with the unit normal vector 
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field family -fit of the FHSs) 

where c~^Xt represents the kinematical contribution to the evolution and C^Ht represents 
the so-called non-kinematical contribution defined by 

N,. 



We see that the non-kinematical evolution (NKE) of the spatial geometry takes the form 
of an "expansion" . Besides we see from equation (2) that the evolution of Nt with time 
may be written in the form 

^-^ = c x, + c y, (4) 

The spht-ups defined in (3) and (4) are necessary to be able to construct gt from gt [1]- 
A special property of the quasi-metric theory is that operationally defined atomic 
units are defined to vary in space-time. That is, gravitational quantities may get an 
extra formal variability when measured in atomic units (and vice versa). To quantify the 
notion of formally variable units we introduce the scalar field telling how atomic time 
units vary in space-time. Measured in atomic units gravitational quantities then get a 
formal variability as some power of "iff By definition c and h are not formally variable; 
this means that atomic time units vary similarly to atomic length units and inversely to 
atomic mass units. However, for convenience we define the gravitational constant G to be 
a constant measured in atomic units. But if we do this we must separate between active 
mass, which is measured dynamically, and passive mass (passive gravitational mass or 
inertial mass). By dimensional analysis it is found [1] that active mass rrit varies formally 
as ^^7^ measured in atomic units (but passive mass does of course not vary). By definition 
we have ^t=j^, so the formal variation of the units follows directly from the spatial 
scale factor of the FHSs as seen from (1), that is 

mt,t = ( 7 + -TT j"^*' ^ '^t^i ^ ^ a^jm, (5) 

where a^F is the 4-acceleration of the FOs in the family g^. Local conservation laws 
involving projections into and normal to the FHSs of the (active) stress-energy tensor 
take the form (in component notation) 
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where Ca^ denotes Lie derivative of spatial objects in the direction normal to the FHSs 
and 'I' denotes spatial covariant derivation. (A "hat" denotes an object intrinsic to the 
FHSs.) Sec [1] for a derivation of these equations. Also postulated in [1] are the field 
equations (with k=^^) 



P{t)jl — ^lt)j\i ~ ^t,j — i<'T(t)ji, (9) 

where Rt is the Ricci tensor family and is the extrinsic curvature tensor family of the 
FHSs corresponding to the metric family (1). {Kt is the trace of K*.) 

The geometry intrinsic to the FHSs obtained from equation (1) has the curvature 

H{t)ij = c~^(a^|fe - -]^^-^)ht)i3 - c~'^aj^iaj^j - c~'^aj^i\j, (10) 



{Ntcty 



+ 2c- - 4c-'a!rK, (H) 



where is the metric family intrinsic to the FHSs, is the Einstein tensor family 
intrinsic to the FHSs and Pt is the Ricci scalar family intrinsic to the FHSs. 



3 Metrically static, spherically symmetric interiors 

As shown in [2], the gravitational field exterior to an isolated, spherically symmetric 
body can be found from the field equations. In this case the metric families gt and gt are 
spherically symmetric in a S^xR-background. Furthermore they are "metrically static"; 
by definition this means that neither Nf nor (the lapse function field of the FOs in g^) 
depend on (or t) and that the shift vector field appearing in equation (1) vanishes in a 
suitably chosen GTCS. It would seem reasonable to look for interior gravitational fields 
of this form. 

To find the general form of gt and gt for the spherically symmetric, metrically static 
case we start with the general expression (1) and introduce a spherically symmetric GTCS 
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r, 9, (j)}. Then it is straightforward to show that the metric famihes gt and gt can be 
represented by the famihes of hne elements ('=^) 

cWj = 5(r)((ix°)2 - {-f(A{r)dr^ + r^dQ^), 

to ^ ' 

c^drl = B{r){dx''f - {-y(A{r)dr^ + r''dn''), (12) 



^0 

1 — r 



2B(r) 



^ (13) 

where dfl^=d9^ + sin^^d0^ and to represents some arbitrary reference epoch setting the 
scale of the spatial coordinates. Furthermore x'^ = ct. However, x'^ and ct must not be 
interchanged in (12) due to the particular form of the connection coefficients, see [1] for 
a discussion. 

We now seek non-vacuum solutions of the type (12) of the field equations, and where 
the source is modelled as a perfect fiuid. For the vacuum case it was found in [2] that 
putting a line element ds"^ of the form (12) into the field equations results in equations 
separable in A{r) and B{r). It is then possible to find a series solution. And it turns out 
that for non-vacuum cases one still gets separable equations when using line elements of 
the form (12). But we now have the passive pressure p{r, t) as a new variable in addition 
to A{r) and B{r). (Note that the passive mass density pm{r,t) is also a new variable but 
that it is related to p[r, t) via an equation of state.) The explicit dependence on p in the 
equations makes analytical calculations rather impracticable so the equations should be 
solved numerically. In the following we set up the relevant equations. 



3.1 Analytical calculations and a numerical recipe 

In this section we do some analytical calculations in order to write the equations in a 
form appropriate for numerical treatment. Proceeding with this, from the definitions we 
get 

_2_ _ -2- - ^ 

""^^ ~ 2B' " ''^'■l'' ~ 2B " 25^ ~ 4AB' 

c-sin-^a^,l, = ^-'-^^e\e = ^, A = (^)^|(r-(A - 1) + £), (14) 



^t'^2AB 2AB^ AA^B^ rABJ' 

— — — — T A' 

H{t)rr = r-^(l - A), sv^^'^eH^t)u = H^t)ee = (15) 
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Now the active stress-energy tensor varies formally as ^| when measured in atomic 
units. For reasons of convenience we choose to extract this formal variability explicitly 
What is left after separating out the formal variability from the active mass density is by 
definition the coordinate volume density of active mass pm- (The corresponding pressure 
is p.) To be able to have experimental input it is also necessary to define the passive mass 
density Pm and the corresponding passive pressure p. For the case when the perfect fluid 
is comoving with the FOs (i.e., T(^t)lj — 0) we may then write (using (5)) 

where pm is the density of active mass-energy and p is the active pressure as seen in the 
local rest frame of the fluid. Furthermore, using (16) the local conservation laws (6), (7) 
yield C^l) 

B' 

Pm = P = 0> = -c"^a^r(PmC^ - 3p) = -(p^c^ - 3p)^. (17) 

Equations (17) are valid for any metrically static perfect fluid. But to have experimental 
input we need to specify an equation of state p — p{pm)- We see that we need an 
expression relating p^ (or equivalently, pm) to the passive mass-energy density p^. Such 
a relationship is given by [1] 

r^tPm for ^ fluid of material point particles, 

: (18) 

^Nfpm for a null fluid, 
and a similar relationship exists between p and p, namely 

{j-Ntp for a fluid of material point particles, 
I - (19) 
%Nfp for a nuU fluid. 

To find how the active mass mt varies in space-time, note that wc arc free to choose the 
background value of the active mass far from the source to be itlq. We use this to define 
G as the constant measured in a Cavendish experiment far from the source at epoch to- 
Then, using equations (5) and (14) we find 

mt{r,t) ^ B^l\r)^m^. (20) 

Now we can insert the equations (14), (15) and (16) into the field equation (8). Since 
vanishes identically for the metrically static case [1] equations (9) become vacuous and 
equation (8) yields 

B" B'^ A'B' 2B' A, ^ ^ 

+ ^ = «^(Pmc2 + 3p). (21) 



B 2^2 2AB rB B 
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A second equation can be found from the equation (10) for the spatial curvature. The 
radial component just gives equation (13) again but the angular components yield a new 
equation, which can be combined with equation (21) to give 



'A 



'b 



B 



2B^ 



2rA A, 2 



(rfo 



(22) 



We may now eliminate A{r) and A'{r) from equation (21) by inserting equations (13) 
and (22). The result is (with So=c^o) 



(1 



B(r)EV Bir 



B"(r) 



- (1 



2r2 



rB 



Bir)El 



B 



111 



AB^{r) 



B^{r) 



3r^ 



Bir 



rB'ij-] 
2B(r) 



B'{r) 
2B{r)'Ep B{r) 



Pm{ry + 3p(r) . (23) 



To solve equation (23) numerically for B{r) one may proceed as follows. First specify the 
boundary conditions at the center of the body. Prom (13) we easily see that ^4(0) = 1, 
and moreover we must have -B'(O) = p'(0) = 0, so A'{0) — from (22). Furthermore, 
noting that r~^B' must be stationary near the center of the body we have 



B"{Q) = lim. 



-'B'{r) 



B"{0) = ^{pm{0)c' + 3p{0 



(24) 



where the implication follows from equation (21). 

To specify any particular model, choose the boundary condition p(0) at some arbitrary 
time. Also choose some arbitrary value B{0) as an initial value for iteration. It must be 
possible to check how well the chosen -B(O) reproduces the boundary condition for B(TV) 
at the surface r — TZ of the body. That is, to match the exterior solution we must have 
[2] 



sO 



sO 



2El 87^3 



+ 



(25) 



where Tso is a constant which may be identified with the Schwarzschild radius of the body 
at the arbitrary time to. Hence, by definition we set 



.2Mt„G 



Mt=c-'JJJNt 



{t)±± 



{t)i 



* = ^""tX^ [p-^ + Sp/c^Jr^dr, (26) 



where the integration is taken over the whole body. (The particular form of Mt^ follows 
directly from the "metric approximation" [1] of equation (8) applied to the interior of 
a spherically symmetric, static source when extrapolated to the exact exterior solution 
found in [1].) Note that Mt^ is defined to be the total mass-energy of the body as 
measured by distant orbiters at epoch and that depends on the pressure as well. 
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To have a metrically static system it is necessary to specify an equation of state of 
the type pocp^ (e.g. an ideal gas), since this ensures that and p are independent of t. 
Once a suitable equation of state is specified we may use equations (16), (18) and (19) to 
find pm and p and then integrate (23) outwards from r = using (17) until the pressure 
vanishes. The surface of the body is now reached (that \sr — TZ). If the calculated value 
B{TZ) does not match the boundary condition (25), add a constant to B everywhere and 
repeat the calculation with the new value of S(0). Iterate until sufficient accuracy is 
achieved. 

Once we have done the above calculations for an arbitrary time, we know the time 
evolution of the system from equations (16), (17), (18) and (19). That is, a spherical 
gravitationally bound body consisting of a perfect fluid obeying an equation of state 
pccpm will expand according to the Hubble law. But for perfect fluid bodies obeying 
other equations of state the expansion will induce instabilities; mass currents will be set 
up and such systems cannot be metrically static. However, for the metrically static case 
the gravitational field interior to the body will expand along with the fluid; this is similar 
to the expansion of the exterior gravitational field found elsewhere [2]. Wc arc now able 
to calculate the family inside the body. To find the corresponding family gt one uses 
the method described in [1]. 

4 Discussion 

The possibility that cosmic expansion may be relevant for gravitationally bound bodies is 
not a new idea [4]. However, it has not been possible to fit this concept into the standard 
framework of metric gravity. In this paper we have shown that in quasi-metric theory, 
gravitationally bound bodies may be expected to expand with the Universe whenever 
they are made of a perfect fiuid with an equation of state of the type pocpm- (This is 
easily seen from (12) and the fact that the coordinate radius TZ of the body is constant in 
time.) But if the linear equation of state is not a good approximation one expects that 
the expansion should lead to instabilities in order to restore hydrostatic equilibrium. 

The Earth is not made of perfect fluid. But if its equation of state is close to hnear, 
the Earth should be expanding due to cosmic effects. This yields a number of observable 
consequences. One may, e.g., calculate the expected spin-down rate of the Earth due to 
the expansion; this is quite consistent with observations if these are interpreted within 
the quasi-metric framework [2] and provided that tidal friction is negligible. Furthermore 
palaeo-tidal records, obtained e.g. from sedimentary tidal rhythmities, are in general 
in agreement with predictions calculated from quasi-metric theory as long as no extra 
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theory-dependent assumptions, inconsistent with the quasi-metric framework, are made. 

Newtonian stars for which the equation of state takes the form pocp^, 7>|, are called 
polytropes [5]. In quasi- metric theory it is possible to model polytropes by taking Newto- 
nian limits of the relevant equations but such that the t-dependence remains. However, 
polytropes made of degenerate matter will not be metrically static in such a description. 
But if the hydrodynamical effects on the gravitational field coming from instabilities can 
be neglected, we may solve the field equations for each epoch t assuming that the poly- 
trope is at hydrostatic equilibrium. To do this we take the Newtonian limits of equations 
(16), (17) and (23). One may then straightforwardly show that equation (23) reduces to 

^—p' = -47rGrVmo, (27) 

where pmo is the density field pm at epoch to. For polytropes this equation can be 
used only for the fixed epoch Iq. However, the equation is invariant under the scale 

1 3 

transformation r— >(^) 3T^r=£, Pmo^ij^)^^ Pmo=Pm, G^j^^=^t- This means that by 
first scale transforming equation (27), it may be applied to polytropes and solved for any 
fixed epoch t. It then becomes equivalent to its counterpart in Newtonian theory except 
for a variable Gt- Thus the usual Newtonian analysis of polytropes [5] applies, but with 
Gt variable. And as a consequence wc see that the physical radius C{t) = (^) 37-4£(^q) 
of a polytrope made of degenerate matter will actually shrink with epoch (if 7 > |). 

Of particular interest arc polytropes for which 7 = |, since such stars are models 
for Chandrasekhar mass white dwarfs (WDs). From [5] we easily find that the (pas- 
sive) mass rric and radius Cc of such WDs (with identical central mass densities) depend 
on epoch such that mc{t) = (^)^/^mc(to) and Cc{t) = {^y^'^'^c{to)j respectively. Since 
Chandrasekhar mass WDs are believed to be progenitors of type la supcrnovac, one may 
expect that any cosmic evolution of Chandrasekhar mass WDs should imply a system- 
atic luminosity evolution of type la supernovae over cosmic time scales. However, such 
a luminosity evolution would be inconsistent with their use as standard candles when 
determining the cosmological parameters in standard cosmology: luminosity evolution 
could have serious consequences for an interpretation of the supernova data in terms of 
an accelerating cosmic expansion indicating a non-zero cosmological constant [6], [7]. 

Now quasi-metric theory predicts that the cosmic expansion does neither accelerate 
nor decelerate [1]. Moreover, according to quasi-metric theory, the Chandrasekhar mass 
varies with epoch and this means that type la supernovae may be generated from cosmo- 
logically induced collapse of progenitor WDs. The consequences for type la supernova 
peak luminosities due to the predicted evolution of progenitor WDs are not clear. Since 
the luminosity of type la supernovae comes from 7-disintegration of unstable nuclear 
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species (mainly ^^Ni) synthesized in the explosion, this luminosity could depend criti- 
cally on the conditions of the nuclear burn. That is, the detailed nuclear composition 
synthesized in the explosion and its energetics might depend on progenitor mass. Also 
the presence of more massive ejecta during the explosion could have an influence on su- 
pernova luminosities and light curves. But since there is no way to know exactly how 
these observables are affected by the predicted progenitor evolution without doing de- 
tailed numerical simulations, it is not possible to say whether or not the predictions from 
quasi-metric theory are consistent with the data. 

However, what we can easily do is to see if it is possible to construct a simple lumi- 
nosity evolution which, in combination with the cosmological toy model described in [1], 
yields a reasonable fit to the supernova data. That is, we may try a luminosity evolution 
of the form L(t)oc(^)'^, or equivalently L{z)(x{l + z)~'^ and see if the data are well fit 
for some valuc(s) of e. To check this we use the postulated luminosity evolution to plot 
apparent magnitude mqMT versus redshift for type la supernovae. The easiest way to 
compare this to data is to calculate the difference in apparent magnitude Am between 
our model and that predicted using an empty Friedmann model (which is a piece of 
Minkowski space-time). The result is 



where Lqmt/Lmin represents the luminosity evolution of the source in our quasi-metric 
model relative to no luminosity evolution in an empty Friedmann model. 

Comparing the fits obtained from equation (28) to data we find that the luminosity 
evolution Lqmt = (1 + -2) "'^''^-^m/jv fits the data very well. We conclude that quasi-metric 
cosmology combined with a simple luminosity evolution seems to be consistent with the 
data but that a much more detailed model should be constructed to see if the postulated 
luminosity evolution has some basis in the physics of type la supernovae. 
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